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A flow generator is described in which homogeneous axisymmetric turbu- 
lent air flows with varying and fully controllable degrees of anisotropy, includ- 
ing the much studied isotropic case, are generated by the combined agitations 
produced by 32 acoustic mixers focusing at the center of the system. The ax- 
isymmetric turbulence in a central volume of the size of the inertial scale is 
shown to have negligible mean and shear. The Taylor Reynolds number is about 
480. 

The influence of large scale anisotropy on the turbulence is examined from 
three aspects, namely the velocity structure functions, the velocity correlation 
functions, and the integral lengths. The directional dependence of two different 
second order transverse structure functions, in which one of them has separa- 
tions stretched along the axis of symmetry of the turbulence and the other one 
normal to it, is studied. It is shown that the inertial range scaling exponents, 
determined using the extended-self-similarity procedure, and the Kolmogorov 
constants of the two structure functions are unaffected by the direction in which 
the structure functions are measured. 

As an extension, because of its relevance to the study of intermittency, the di- 
rectional dependence of transverse structure functions of the fourth to the sixth 
order is studied. Despite some issues with measurement noise and statistical 
convergence, some indications are found that anisotropy in the velocity field 



intensifies the asymmetry of the probability density of the velocity increments. 
In addition, some evidence is found that the inertial range scaling exponents of 
the fourth, fifth, and sixth order are independent of the anisotropy 

Finally, it is found that, except in the isotropic case, the second order trans- 
verse velocity correlation functions deviate from each other at the large scale 
with increasing anisotropy A self-similarity argument similar to one found in 
the study of critical phenomena is proposed. It is shown that the argument leads 
to a power-law relationship between the large scale velocity fluctuation and the 
correlation length, with an exponent that depends on the inertial range scaling 
exponent of the turbulence. The data collapse predicted by the self-similarity 
hypothesis is verified. It is demonstrated that the value of the power-law expo- 
nent is consistent with the value of the inertial range scaling exponent. 
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4.1 



The inertial range scaling exponents of the longitudinal and 
transverse structure functions. The first two rows are the longi- 
tudinal scaling exponents , , predicted by |Kormogorov| ( |1941a[ ) 
(see equation 1.49) and She and Leveque ( |1994| ) (see equa- 



tion 



1.54[ >. The third and fourth rows show the longitudinal and 



transverse exponents measured by Dhruva et al. ( 1997[ ) using the 



ESS method (see section 1.4 1 in atmospheric turbulence at R\ be- 
tween 10 4 and 1.5 x 10 4 

The equality between longitudinal and transverse scaling expo- 
nents for sheared and unsheared turbulence and their trends 
with Reynolds number. The acronym TH denotes the use of Tay- 
lor's hypothesis in the experiments 

The landscape of anisotropy for geophysical flows. Anisotropy 
refers to the ratio of streamwise to spanwise to vertical, stream- 
wise to vertical, or zonal to meridional RMS velocity fluctua- 



tions, depending on the context. In Sullivan et al. (2000 ) the RMS 



velocity fluctuation data are unavailable. The anisotropy quoted 

is for the mean wind speed 

A bird's eye view of the anisotropy observed in laboratory flows. 
The symbol d refers to distance between grid center lines, cylin- 
der diameter, pipe diameter, jet nozzle diameter, tunnel width, 
jet width, or the boundary layer thickness, depending on the 
context, u'/v' refers to the ratio of streamwise to spanwise, ver- 
tical to horizontal, or radial to axial RMS velocity fluctuations, 
depending on the context. |Shy et al. ( 1997[ ) did not specify the 
Reynolds number, it is estimated from equation 1.7[ 



The chamber geometry and flow statistics from various speaker- 
driven flow chambers. The symbols from top to bottom are: 
energy dissipation rate (e), Kolmogorov length scale {rj), Kol- 
mogorov time scale (r^), Taylor micro-scale Reynolds number 
(R\). For each study, only the case that gave the highest 
Reynolds number is reported 



For forcing anisotropy ratios (A) ranging from 0.33 to 2.5, the 
table shows the velocity fluctuation ratio according to the two 



limits discussed in the text (equations 4.1 and 4.2) 
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4.2 The table shows the turbulence parameters, for loudspeaker 
RMS amplitude ratios (A) ranging from 0.33 to 2.5. K is the tur- 
bulent kinetic energy, \ « 2 + 2 u'? 2 ), and | U | is (U 2 Z + 2 £/ r 2 ) 1 ' 2 . f rira 
and f zr2 are the average normalized Reynolds stresses measured 
along and normal to the axis of symmetry, respectively (see sec- 
tion |46]). e is the energy dissipation rate estimated according to 



the procedure prescribed in section 5.1 rj and r v are the Kol 



mogorov length and time scales, respectively. A is the Taylor 



length (see section 5.3 ) and R\ is the Taylor scale Reynolds num- 
ber. Data were collected at the center of the soccer ball when the 
axis of symmetry of the forcing is horizontal. Similar results (not 
shown) were obtained when the axis of symmetry is vertical. . . 

5.1 The numerical data for radial scaling exponents (Q )/ axial scal- 
ing exponents (C2 ), and the ratio between the two 

5.2 The value of the scaling exponent of transverse second order 
structure functions, £ 2 , measured using the ESS method for vari- 
ous laboratory flows 

5.3 The numerical data for the ratio of the Kolmogorov con- 
stant for transverse structure functions measured in nearly 
orthogonal directions, computed from averaging the ratio 
D zz (ri)/D r2r2 (z) over a range of scales between 7 and 70 mm 
(second column), as well as from taking the ratio between 
the maxima of appropriately compensated structure functions, 
m.ax.(D zz /r\^)/ max(D r2r2 / z 2//3 ) (fourth column). The values in 
the third and the last columns are the corresponding estimated 
error. 

6.1 The numerical integral of the PDF and its moments from the sec- 
ond to the sixth order of an isotropically forced turbulent flow, 



shown in figure 6.4 The isotropy ratio was u' z /u' r = 0.94 and the 



6.2 The numerical integral of the PDF and its moments from the sec- 



ond to the sixth order of figure 6.5 The flow anisotropy was 



IT031 



Reynolds number was R\ = 480 (see table 4.2) 11071 



u'Ju' r = 0.60 at R x = 536 (see tableful) ITT31 

6.3 The values of the fourth order radial and axial scaling exponents, 
£4 and ([ z \ measured in our experiment using ESS. The last col- 
umn shows the ratio between the two 11181 

6.4 The numerical data for the fifth order radial and axial scaling 
exponents, Q r \ Q*\ and the ratio between the two 11181 

6.5 The numerical data for the sixth order radial and axial scaling 
exponents, (q \ Q , and the ratio between the two 11191 



x 



The inertial range transverse scaling exponents, ( p , measured in 
the current experiment using the ESS method and those mea- 
sured by |Shen and Warhaft| ( pBo2| ), denoted as S&W (2002) above, 



in an unsheared wind tunnel homogeneous and isotropic turbu- 
lent flow at two different R\ of 134 and 863. 



The numerical data for axial integral lengths (L z ), radial integral 
lengths (L r ), and the ratio between the two 
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1.1 The figure shows the range of eddy sizes t from the largest (L ) 
to the smallest (77) at very high Reynolds number. The meaning 

of each length scale is explained in the text [7] 

1.2 The second-order longitudinal (D u (r)) and transverse (-D22O") 
and D 33 (r)) velocity structure functions measured in a turbu- 
lent boundary layer with R\ = 1450. When compensated 
appropriately, the structure functions within the inertial range 
of scales are a measure of the value of the Kolmogorov con- 



stant. Reprinted with permission from Saddoughi and Veeravalli 



( [19941 ): /. Fluid Mech., 268:333-372, Copyright (1994) with permis 



sion from Cambridge University Press. 



2.1 The 'soccer ball' chamber. Also shown is the layout of the 
Laser Doppler velocimetry probes in the experiment. A two- 
component probe aligned along the xi-axis measured the com- 
ponents of velocity in the x 2 and 23 directions, and an additional 
one-component probe aligned along the x 2 axis measured veloc- 
ities in the x 3 direction [39] 

2.2 An overview of the experiment showing the 'soccer ball', the 
amplifier unit, and the measurement apparatus. The amplifier 
received signals from the multichannel high-speed voltage out- 
put device on the computer through a BNC connector block, and 
sent the signals to the loudspeaker array on the soccer ball. A 
laser Doppler velocimetry (LDV) system collected velocity sig- 
nals and saved the data on a computer (not shown) [40] 

2.3 A loudspeaker with a conical nozzle |42] 

2.4 (a) A 40-second time series of the voltage applied to the loud- 
speakers, (b) A part of the time series, zoomed in to reveal the 

50 Hz carrier wave |46] 

2.5 The spectrum of the voltage applied to the loudspeakers (a) be- 
fore filtering and (b) after filtering. The solid line in (b) is pro- 
portional to a power law with an exponent of —2 |47] 

2.6 A sketch showing the assignment of RMS amplitude for a given 
loudspeaker. The vector pointing from the center of a sphere, O, 
to the center of a loudspeaker, S ir intersects an inscribed ellipsoid 
at Ei. The distance of Ei from the center of the sphere sets the 
RMS amplitude 6, for the given loudspeaker |49] 

2.7 The figure shows (a) a truncated time series of one component of 
the fluid velocities in response to the forcing and (b) its spectrum 
calculated using the Lomb algorithm for unequally spaced data. [51] 

2.8 The Raveland amplifier. [52] 

2.9 The amplifier unit [53] 
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2.10 A schematic of a laser Doppler velocimeter. 

2.11 The figure illustrates the change in received frequency caused by 
the motion of the particle in a single-beam scattering process. . . 

2.12 The figure illustrates the dual-beam scattering process 

2.13 A schematic of a Bragg cell. A transducer generates an acous- 
tic wave, resulting in a moving periodic refractive index grating 
that diffracts the incident light similar to in Bragg diffraction. 
The arrow points in the direction of travel of the acoustic wave. . 

2.14 An overview of a TSI LDV system. Figure adapted from TSI 
LDV/PDPA system installation manual, page 73 

2.15 The effect of temporal binning on the value of the velocity cor- 
relation, g(r), near zero separation. The values in the legend are 
the separations in millimeters 

2.16 The structure function as a function of temporal bin width. For 
legend, see previous figure 

2.17 The number of samples yielding the velocity correlation as a 



function of temporal bin width. For legend, see figure 2.15 



2.18 The schematic shows the LDV probe volume dimensions. 



3.1 The schematic shows the coordinate system in the laboratory 
frame. s 3 was aligned with the vertical. The axis of symmetry 
of the forcing lay along either x[ or x' 3 , depending on the config- 
uration. Our measurements were made at points lay along x\. . . 

3.2 The figure shows the orientations of the body coordinate system 
(in red) of the forcing (ellipsoidal surface in yellow) with respect 
to the laboratory frame (gray), and the conventions we used: (a) 
when the axis of symmetry of the forcing lay close to the x 3 axis 
of the laboratory frame, and (b) when the axis of symmetry lay 
close to xi axis of the laboratory frame. The coordinate system 
(xi, x 2 ,x 3 ) was fixed in the laboratory frame, and the coordinate 
system (r 1; r 2 , z) was fixed with respect to the symmetry of the 
forcing 

3.3 The figure shows the orientations of the measured velocity com- 
ponents with respect to the forcing coordinate system when the 
separation vector lay (a) along the axis of symmetry, or (b) per- 
pendicular to the axis of symmetry. 
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The velocity fluctuation ratios as a function of the forcing 
anisotropy, A. The circular symbols (O) are a/ (uf (0) ) / (u' 2 (0) ) 
and the triangular symbols (A) are *J (uf (0) ) / (w' r 2 (0) ) . Sym- 
bols in the region shaded in blue have oblate shape asymetry 
whereas those in the region shaded in red have prolate shape 
asymmetry. Spherical symmetry lies at the intersection of the 
two regions. Data are collected with the axis of symmetry lying 
along x' 3 . We obtained similar results (not shown) when the axis 

of symmetry lay along x[ 

For three values of the anisotropy measured at the cen- 
tre of the ball, v / « 2 (0))/« 2 (0)) / we plot (a) the ratio 
of the two radial components of the velocity fluctuations 
y/(u!? 2 (0, 0, z))/« 2 (0, 0, z)) and (b) the ratio of the axial fluctu- 
ations to the radial fluctuations, a/ (u'?(ri, 0, 0))/(u^(ri, 0, 0)) at 
various points moving away from the center of the ball. The val- 
ues in these legends, and in all others in this thesis, are those of 

the data in (b) for r x = 

The comparison of fluctuations at various distances from the 
center of the soccer ball to those in the middle, (a) Radial 
fluctuations along axial separations, a/ (w' r 2 (0, 0, z))/(u'? 2 (0, 0, 0)), 
for anisotropy ratios, (uf (0)) / (u!?(Q)}, ranging from 0.59 
to 1.71. (b) Axial fluctuations along radial separations, 
a/ (uf(r 1 ,0, 0))/(uf(0, 0, 0)), for the same range of anisotropy. 
The anisotropy ratios in the legend are measured at the center 

of the soccer ball 

The mean flow as a fraction of the fluctuating velocity (a) 

(U r2 (0,0,z))/y/«( Q >°> z ))> and ( b ) (^(ri,0,0))/V« 2 (n,0,0)) 
for axial-to-radial velocity fluctuation ratios ranging from 0.59 to 
1.71. The values in the legend are the anisotropy measured at the 

centre of the soccer ball 

The cross correlations of orthogonal velocity components at var- 
ious distances from the center of the soccer ball (a) along the 
axis of symmetry, and (b) along an axis perpendicular to the axis 
of symmetry of the forcing. The values in the color bar are the 
anisotropy measured at the center of the soccer ball for various 

values of the anisotropy of the forcing 

The figure shows the normalized Reynolds stress of orthogonal 
velocity components at various locations (a) along the axis of 
symmetry of the forcing, and (b) perpendicular to it. The val- 
ues in the color bar are the anisotropy measured at the center of 
the soccer ball for various values of the anisotropy of the forcing. 
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5.1 The figure shows the (a) axial structure functions, D r2r2 (z), and 
(b) radial structure functions, D zz (ri). In both cases, the structure 
functions are divided by the inertial range scaling predicted by 
Kolmogorov. The error bars are smaller than the symbols. The 
values in the legends are the anisotropy measured at the center 

of the soccer ball for various values of the anisotropy of the forcing. [90] 

5.2 The ratios between structure functions in different directions but 
in the same flow, D zz (r)/D r . 2r2 (z), approach the isotropic value 
of one when the separation distance, t\ = z, decreases. Results 
for data taken at other values of the large-scale anisotropy are 
consistent with those shown here. The values in the legends are 

the anisotropy measured at the center of the soccer ball [91] 

5.3 The third-order structure functions with separations in the direc- 
tion (a) along the axis of symmetry, and (b) perpendicular to the 
axis of symmetry. The values in the color bar are the anisotropy 
measured at the center of the soccer ball [93] 

5.4 The third-order moments of the absolute values of the velocity 
increments with separations in the direction (a) along the axis 
of symmetry, and (b) perpendicular to the axis of symmetry. The 
dashed lines are proportional to a power law with exponent 1.17. 
The values in the color bar are the anisotropy measured at the 
center of the soccer ball [94] 

5.5 The extended-self-similarity plots for the second-order trans- 
verse structure functions. The dashed lines are proportional to 
the power laws obtained from least squares fits to data. The val- 
ues in the color bar are the anisotropy measured at the center of 
the soccer ball for various values of the anisotropy of the forcing. 

The error bars are smaller than the symbols [96] 

5.6 The scaling exponents for the transverse structure functions 
measured using the extended-self-similarity method. The O 
symbols are the scaling exponents for D r2r . 2 (z) and the □ sym- 
bols are for D zz (r x ) . The last two points in the radial data set have 
large uncertainties due to spatial resolution. They are marked in 
gray to show the general trend of the exponents 197] 

5.7 The figure shows Cl^/Cl / the ratio between the scaling expo- 
nents of the radial structure functions and those of the axial 
structure functions. The last two points with large uncertainties 

due to spatial resolution are marked in gray. [98] 
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5.8 The figure shows the ratio between structure functions in differ- 
ent directions in the same flow, D zz (ri) / D r2r2 (z), for all values 
of the anisotropy of the forcing. The data are the same as in fig- 
ure 5.2 They are plotted separately to display individual trends. 
The two vertical dashed lines demarcate the range of scales be- 
tween 7 and 70 mm where the compensated structure functions 
displayed a reasonable plateau. The values in the color bar are 
the fluctuations anisotropy measured at the center of the soccer 
ball for various values of the anisotropy of the forcing 

5.9 The compensated axial and radial transverse structure functions 
for all values of anisotropy. Only the data in the range of scales 
between 2 and 50 mm are used in the determination of the ratio 
of the Kolmogorov constant. Solid lines are the second-order, 
least-square, linear-log polynomial fits to the structure function 
data. Solid symbols are the maxima of the fitted curves. The 
error bars are smaller than the symbols. The values in the color 
bar are the anisotropy measured at the center of the soccer ball 
for various values of the anisotropy of the forcing 

5.10 The ratio between Kolmogorov constants measured in differ- 
ent directions C^/C^ . The open circular symbols (O) are the 
ratio between the maximum values of the compensated radial 
and axial structure functions. The open square symbols (□) are 
{D zz (ri) / D r2r . 2 (z)} , the ratio of radial to axial structure functions 
averaged over a range of scales 7 < r < 70 mm 



6.1 The extended-self-similarity plots for the fourth-order trans- 
verse structure functions. Dashed lines are proportional to 
power laws obtained from least squares fits to data. The val- 
ues in the color bar are the anisotropy measured at the center of 
the soccer ball 

6.2 The extended-self-similarity plots for the fifth-order transverse 
structure functions. Dashed lines are proportional to power laws 



obtained from least squares fits to data. For legend, see figure 6.1 
6.3 The extended-self-similarity plots for the sixth-order transverse 
structure functions. Dashed lines are proportional to power laws 



obtained from least squares fits to data. For legend, see figure 6.1 
6.4 The probability density of velocity increments, 5u r2 (z), and their 
moments for an inertial range separation z/rj ~ 130 and spheri- 



cally symmetric forcing, u' z /u' r = 0.94, at R\ = 480 (see table 4.2). 

Here x = 5u r . 2 (z)/ ' {\6u r2 (z)\ 2 ) 1/2 

6.5 The probability density of velocity increments, 5u r2 (z), and their 
moments for an inertial range separation z/r] ps 130 and oblate 



shape forcing, u' z /u' r = 0.60, at R\ = 536 (see table 4.2). Here 

x = 5u r2 {z)/(\5u r2 {z)\ 2 ) 112 
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6.6 The kurtoses as a function of separations, (a) The axial kur- 
toses, K A (z) = (|5w r2 (z)| 4 )/(|5M r2 (z)| 2 ) 2 , and (b) the radial kur- 
toses K±(r\) = (|5n z (ri)| 4 )/(|5M 2 (ri)| 2 ) 2 . The dashed lines are 
at the gaussian value of 3. The values in the color bar are the 
anisotropy measured at the center of the soccer ball. Error bars 

(not shown) are smaller than the symbols 11151 

6.7 The hyper-kurtoses as a function of separations, (a) shows 
K 6 (z) = (\5u r2 (z)\ G )/(\5u r2 (z)\ 2 ) 3 and (b) shows K 6 ( n ) = 
{\5u z (ri)\ 6 } /(\5u z (ri)\ 2 ) 3 . The dashed lines are at the gaussian 
value of 15. The values in the color bar are the anisotropy mea- 
sured at the center of the soccer ball. Error bars (not shown) are 
smaller than the symbols 11161 

6.8 The figure shows (4 /Q / the ratio between the fourth-order 
scaling exponents of the radial structure functions and those of 
the axial structure functions. The last two points with large un- 
certainties due to spatial resolution are marked in gray. 11191 

6.9 The figure shows Cs /Cw the ratio between the fifth-order scal- 
ing exponents of the radial structure functions and those of the 
axial structure functions. The last two points with large uncer- 
tainties due to spatial resolution are marked in gray. 11201 

6.10 The circular symbols (O) are our measurements of the inertial 
range scaling exponents, ( p , of the transverse structure functions 
using ESS method. The triangular symbols (A) are the exponents 
measured by Shen and Warhaft ( 2002[ ) in an unsheared wind tun- 



nel homogeneous and isotropic turbulent flow at an R\ of 134, 
and the square symbols (□) are those measured by them at an 
R\ of 863. The dashed line is the Kolmogorov prediction (equa- 
tion 1.49). The solid line is the She-Leveque prediction (equa- 
tion 



1.54) [ml 



7.1 The figure shows the 8 pairs of transverse correlation functions 
obtained in 8 different flows with varying degree of large scale 
anisotropy. Circles (O) are for g r2 r 2 ( z ) an d crosses (x) are g zz {ri). 
The error bars are smaller than the symbols. The values in the 
color bar are the fluctuations anisotropy measured at the center 

of the soccer ball for various values of the anisotropy of the forcing JT301 

7.2 The extrapolation of the partially measured correlation function. 
The blue circular symbols (O) are the experimental data. The 
dashed line is the exponential extrapolated from the data. The 
area in the region shaded in blue is numerically integrated and 
the area in the region shaded in pink is analytically calculated. 

The sum of the two yields the integral length 11331 
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7.3 The sum of squares of residuals in the determination of x max that 
yields the exponential from a least squares fit of an exponential 

tail to the correlation function data 11341 

7.4 The variation in the value of the exponent with the maximum 
separation in the least squares fit of exponential tails to the data. 11351 

7.5 The main figure shows the ratio of integral scales, L z /L r , as a 
function of the anisotropy in the fluctuating velocities, u' z /u' r . 
The solid line is the least squares fit to the data and the dashed 
line is the power law with exponent —3 predicted by the Kol- 
mogorov inertial range scaling argument. The inset shows the 
same data on a logarithmic scale. The straight line in the inset 

has a slope of —2.56 11361 



7.6 The figure shows the same correlation functions from figure 7.1 
plotted against the reduced variables, r 1 u z ~ 2 ^ and z Ur 2 ^' '■ Here, 
C = 0.73. Circles (O) are f or g r2T2 (z) and crosses (x) are for g zz {r\). 
The values in the color bar are the fluctuations anisotropy mea- 
sured at the center of the soccer ball 11381 



A.l Contour for the integral A.5 The pole at zero has been avoided 
by following the small circular quadrant near the origin. The 
integrals over the two circular quadrants tend to zero as R — > oo 
and 5 — > 0. Only the integrals over the line segments contribute 



to the integral [01 IT431 
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CHAPTER 1 
INTRODUCTION 



Homogeneous and axisymmetric turbulence, a simple and unique class of 
anisotropic turbulence, has eluded detailed investigations due to the lack of 
theoretical descriptions and flow generators whose large scale forcing can be 
tuned systematically. The intent of the present experimental study is to expound 
on some ideas for making quantitative measurements of homogeneous and ax- 
isymmetric turbulent flows, in hopes that a better way of characterizing such a 
flow may be found in the future. In his valuable review article on experimen- 
tal methods in turbulence research for the Handbuch der Physik, [Corrsin (1963b) 
states that 

As in other areas of science, the goals of experiment are, loosely, 
of two kinds, a) exploratory and b) to confirm or disprove theories. 
The former puts a premium on imagination, with accuracy often sec- 
ondary; the latter requires accuracies at least sufficient to distinguish 
among competing theories. In exploratory measurements of statisti- 
cal properties of turbulent flows 20% accuracy is sometimes satisfac- 
tory. The more permanent crucial data, especially where signal-to- 
noise ratio is large, are taken with perhaps 2 to 10% uncertainty. 

It is in this spirit that we would like to invite the reader to join in our explo- 
ration of this fascinating field of anisotropic turbulence. As a prelude to the 
experimental study, let us begin with a short theoretical description of turbu- 
lence. 
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1.1 Kolmogorov theory 



The equation of motion of an unforced incompressible fluid, the Navier-Stokes 
equation (in nondimensional form) 



9u _ „ 1 

— + u -Wu = - Vp + — V 2 u . 

at Re 



(1.1) 



is a nonlinear partial differential equation. Here, u(x,t) and p(x,t) are the di- 
mensionless fluid velocity and pressure, respectively, which are functions of 
space, x, and time, t. The dimensionless parameter Re is the Reynolds num- 
ber. Suppose the fluid has a well-defined characteristic length scale, J2?, and a 
characteristic velocity scale, % , the Reynolds number is defined as 



Re 



v 



(1.2) 



where v is the kinematic viscosity of the fluid. At low Reynolds numbers, the 
inertial forces (the terms on the left hand side of the equation of motion) are neg- 
ligible and solutions to the Navier-Stokes equation exist. The flows are smooth 
and laminar, like honey flowing down a plate. At high Reynolds numbers, the 
flows become turbulent and no analytical solution satisfying realistic boundary 
and initial conditions has yet been found. A complete process in solving a phys- 
ical problem usually involves finding the right governing equations, followed 
by finding the solutions to these equations. In turbulence, we have only the 
equations and the solutions are still missing. Nelkin (1992), Warhaft ( 2002[ ), and 
chapter 41 of Feynman et al. (1963 ) give concise introductions to the turbulence 
problem. For a slightly idiosyncratic review of the milestones of turbulence, the 



reader could refer to Lumley and Yaglom (2001 ) 



For simplicity, research in turbulence has focused on idealized turbulent 
flows detached from the influence of externally imposed boundary. As a re- 



suit, there are no apparent velocity and length scales with which one could use 
to define the Reynolds number. For this reason, the Reynolds number is con- 
structed from the statistics of the turbulence. The characteristic velocity scale is 
taken as the root-mean-square (RMS) of the velocity fluctuations 



u 



1 3 



a, 



1/2 



(1.3) 



and the characteristic length scale is taken as the integral length, L, which is the 
integral of the velocity autocorrelation function (see section 7.1) . This yields a 
Reynolds number 

Re L = — , (1.4) 



intrinsic to the turbulence and independent of the external geometry of the flow. 



More frequently, the Taylor scale Reynolds number (Taylor| 1935) 



Ry 



u'X 



(1.5) 



based on a different intrinsic length scale is chosen to describe the turbulence 
level. The Taylor scale, A, is the separation at which the parabolic expansion of 
the velocity autocorrelation function near the origin becomes zero, and, in the 
case of isotropic turbulence (to be described shortly), is given by ( Taylor[ 1935} 



A 



15 vu 



12 



(1.6) 



where e is the energy dissipation rate per unit mass of the fluid. The Taylor scale 
Reynolds number is related to Re^ by ( Taylor} 1935) 



R x ^ (15Re L ) 1/2 



(1.7) 



It is routinely employed in the literature and is useful for comparison to the 
plethora of theoretical, numerical, and experimental results. 
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Taking the divergence of equation 1.1 and using the incompressibility con- 
dition, V • u = 0, yields a Poisson equation for the pressure field 

V 2 p= -V-(u- V)u. (1.8) 



Equation 1.8 shows two features of the pressure field. Firstly, the pressure term 



in the Navier-Stokes equation is a nonlinear term of the same order as the ad- 
vective term, (it • V) u. Secondly, the pressure at a given point is determined by 



the velocity field everywhere in space at the same time. Equation 1.8 thus high- 
lights the crux of the turbulence problem - the Navier-Stokes equation is nonlin- 



ear, spatially nonlocal, and contains no small parameter in the theory (Warhaft 



2002[ ). Any small perturbative expansion around the linear part of the equation 
will be amplified by the nonlinear terms and eventually diverge. 

To make progress, we turn away from the equation of motion and construct 
phenomenological models based on physical assumptions. A recurring theme 
in many branches of physics is the search for universality. Physical systems 
exhibiting universality display macroscopic phenomena that are independent 
of the microscopic details of the systems. The emergence of simple continuum 
laws, like the Navier-Stokes equation, from complex underlying microscopic 
interactions would not have been possible had the fluid motions depended in 
great detail on the atomic interactions in the fluid molecules. In turbulence, 
there is an additional layer to the universality of the physical phenomena we 
are seeking. The universality must not only be unaffected by the quantum- 
mechanical properties of the fluid, but it must also be impervious to the influ- 
ence of the initial and boundary conditions of the flow. This is the vision of 
the forerunners of turbulence theorists (see Taylor| 1935} Kolmogorov 1941a 



Prandtl} [1945} |Heisenberg| [1948| |Weizsacker| [1948) |Onsager[ |1949) . These re 



lated works, built upon Richardson's idea of the energy cascade (Richardson 
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1922) ), consider localized turbulent motions (eddies) as the putative mediators 



of energy. Turbulence is viewed as a hierarchy of eddies spanning over a wide 
range of length scales. Each eddy of size t has a characteristic velocity u(£) and 
timescale r(£). Larger eddies may contain smaller eddies. The largest eddies, of 
typical size L , are vigorously created by the forcing mechanisms of turbulence. 
Feeding on the energy they received from the forcing, they become unstable and 
break up into successively smaller unstable eddies. As the cascade proceeds, en- 
ergy is being passed from the larger eddies to the smaller ones, until the effect 
of viscosity on eddies of a particular size becomes significant, and energy is 
dissipated as heat. 



Kolmogorov ( 1941a[ ), building on the foundation of statistically homoge- 



neous isotropic turbulence laid down by Taylor ( 1935 ), is credited for quanti- 
fying Richardson's energy cascade and for introducing the concept of univer- 
sality in turbulence. The tenets of homogeneous isotropic turbulence put forth 



by Kolmogorov ( 1941a[ ) can be summarized as follows. The presentation has 



been abridged from Pope| ( 2000) ) 



(a) For very high Reynolds number turbulent flows, the small-scald^]turbulent 
motions are statistically isotropic. 

This hypothesis is usually referred to as the postulate of local isotropy. 
A turbulent flow is locally isotropic if it is locally homogeneous and if it is 
invariant with respect to rotations and reflections of the coordinate axes. 

(b) The statistical properties of the small scales (i <C L ) is governed univer- 
sally by the energy dissipation rate per unit mass, e, and the fluid viscosity, 



v. 



1 In the language of turbulence, the term 'scale' is suggestive of a characteristic length scale 
for the size of an eddy. 
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The range of scales i <C L is known as the universal equilibrium 
range, see figure The motion of eddies in the energy-containing range, 



£ ~ L , is shaped by the initial and boundary conditions, and is therefore 
anisotropic and non-universal. 

Kolmogorov gave a prescription for calculating the length, time, and 
velocity scales of the smallest eddies 

77 = I/ 3/4 e -i/4 j (19) 
7 - jj = I/ i/a e -i/2 > (L10 ) 

^ = ^/4 e l/4_ (ln) 

They are appropriately referred to as the Kolmogorov length, time, and 
velocity scales. The Reynolds number formed by these scales u^ri/u is 
unity, the smallness of which indicates that the motion of eddies at scales 
t < r] is laminar. 

(c) The statistical properties of those scales in the inertial range is determined 
by e only. 

The universal equilibrium range is split into the inertial range {q <C 



t <C L Q ) and the dissipation range (£ < rj), see figure 1.1 Eddy motion in 
the inertial range is neither affected by molecular viscosity nor the forcing 
mechanism. 

Recasting this hypothesis in the cascade picture, this means that the 
rate of energy input (e) received by the largest eddies is equal to the rate 
of energy transferred in succession from the largest eddies to the smallest 
eddies. This rate, in turn, is equal to the energy dissipated by the smallest 
eddies. An implicit assumption in this premise is that e is independent 
of v. Experimental evidence for this assumption has been collected by 



Universal iequilibrium range 
Dissipation range i Inertial range 




^ L 

Figure 1.1: The figure shows the range of eddy sizes £ from the largest (L ) to the 
smallest (77) at very high Reynolds number. The meaning of each length scale is 
explained in the text. 



Sreenivasan (1984) and, more recently by Pearson et al. (2002). 



The validity of Kolmogorov's hypotheses is constantly under debate (e.g. 
Frisch[ 1995 ). Its limitation to describing only theoretically ideal homogeneous 
and isotropic turbulence has not only very little relevance to anisotropic and 
non-infinite Reynolds number turbulence occurring in geophysical environ- 
ment and in the laboratory but the practical aspect of these hypotheses is also 
in question. For example, it is not clear to which statistical quantities they can 
be applied (see Nelkin| 1994[ ). There is, however, considerable predictive power 
in them and a myriad of the existing work on turbulence rests on them. As an 
illustration, consider the second order velocity structure function 



D i:i (r) = (6ui{r,x) 5u j (r,x)) 



8ui(r, x) is the velocity increment 

5v,i(r, x) = u'Ar + x) — u'Ax) 



(1.12) 



(1.13) 



where v! i is the fluctuation in the z-th velocity component, r is the spatial in- 
crement and x is an arbitrary reference point in space. The symbol (■ ■ • ) is an 



ensemble average in theoretical analysis, and a time average in experiment. The 
commonly studied one-dimensional structure functions are the projections onto 
Cartesian r-axes of Z)y(r). These functions are important because they are ex- 
perimentally measurable. For simplicity, the coordinate system is chosen such 
that the separation r is in the x\ direction. The two other mutually orthogo- 
nal directions x 2 and x 3 are then normal to r. For homogeneous and isotropic 
turbulence, -Dy(r) is related to the longitudinal structure function 

D n (r) = (K(r, 0, 0) - ^(0, 0, 0)) 2 ) , (1.14) 

and the transverse structure functions 

D 22 (r) = ((u' 2 (r, 0, 0) - u' 2 (0, 0, 0)) 2 ) , (1.15) 

D 3 3(r) = (K(r, 0, 0) - u' 3 (0, 0, 0)) 2 ) , (1.16) 

by the equation 

Ai(r) = D 22 {r) + [D u {r) - D 22 {r)\ ^ , (1.17) 



where D 22 (r) = D 33 (r) by symmetry. See e.g. Pope (2000) for a derivation. 



For an incompressible fluid, the transverse structure functions are related to 
the longitudinal structure function 

D 22 {r) = D 33 (r) = (l + T - ^\d 1x {t) , (1.18) 

so that only one scalar function, either -Dn(r) or D 22 {r), is needed for describing 
structure functions with distances measured along arbitrary direction. 

In the inertial range, 77 «C r «C L , an application of dimensional analysis 
and the universality hypothesis leads to the scaling for the structure functions 
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at sufficiently high Reynolds numbers (Kolmogorov 1941a[ ) 



jDll (r)=C 2 (er) 2 / 3 , 
D 22 {r) =D m (r) = ^C 2 (e rf\ 



(1.19) 
(1.20) 



where C 2 is presumed to be a universal scaling constant, known as the Kol- 
mogorov constant. These power laws are the first nontrivial predictions of Kol- 
mogorov's theory^} (See, however, Frisch ( 1995 ), chapter 6, for a modern deriva- 
tion assuming only the symmetries of the Navier-Stokes equation, i.e. time and 
space invariances, rotational invariance, Galilean invariance, and scale invari- 
ance.) 



The power-law behavior (equations 1.19 and 1.20 1 has its Fourier space coun 



terpart. This can be seen by applying the Wiener-Khinchin theorem (e.g. equa- 
tion 4.52 of Frisch ( 1995[ )), which states that the correlation function 



Rij(r) = (ui(r + x)u' j (x)) , 
and the three-dimensional velocity spectrum tensor 



(1.21) 



(1.22) 



are Fourier transforms of each other 



R tj (r) 



dk §ij(k) e 



ik • r 



(1.23) 
(1.24) 



Here, the Fourier representation of itj(fc) was introduced to turbulence theory 
by |Heisenberg| ( 1948 ) and it is written as 



uAk) 



(2 7T) 



dru'(r)e- lk ' r . 



(1.25) 



2 An experimental discovery of the law of pair dispersion by Richardson ( 1926 ) predates Kol- 
mogorov's theoretical discovery. 



The asterisk in u*(k) denotes complex conjugate and reality of u'^r) implies 
u*{k) = Ui(-k). 

Just as with the structure function D^r), is usually projected onto 

Cartesian fc-axes to yield the experimentally accessible one-dimensional spectra 



/oo roo 
dk 2 / dfoQijik) 
-oo J — oo 



(1.26) 



Using the Wiener-Khinchin theorem, |Taylor (1938) showed that the one- 
dimensional spectra and the correlation functions form a Fourier transform pair 



Ru(r) 



71 



Rij(r) e~ ikir dr 



1 f°° 

^ J -oo 



(1.27) 
(1.28) 



In the above, r is taken asre^ The normalization factor is chosen so that when 



r = 0, equation 1.28 gives, for example when i = j = 1 

POO 

Ru(0) = K 2 ) = / E 11 {h)dk 1 . (1.29) 
Jo 

Here, E n (k 1 ) is extended to negative frequencies by E X1 {— kx) = Eu(ki). If the 
turbulence is locally isotropic]^], the correlation function Rij{r) is then linked to 
the structure function D^r) by 



D ij {r)=2{R ij {0)-R ij {r)), 



(1.30) 



which, when equation 1.28 is substituted into the above, yields 

/oo 
(l-e^^Od*!. 
-oo 



(1.31) 



In two short notes that predated his prediction for the structure function, 
Kolmogorov presented relevant mathematical description for turbulence energy 



3 If the turbulence is locally isotropic, then (u-(r) u'j(r)) = (u-(0) u'j(0)) = Rij(0), and also 
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spectra and their connection to structure functions (see Kolmogorov 1940a|b[ ). 
He considered the case when the spectra follow a power law, taking the longi- 
tudinal case, i — j = 1, as an example 



E u (h) = B \h\ 



-(7+1) 



(1.32) 



where B is a positive constant. A substitution of the above into equation 1.31 
with a change of variables, x = k\r , yields the following power law for the 
longitudinal structure function 



Dn(r) = BA^\r\ 



(1.33) 



with the dimensionless integral 

— 



l-e l 



x 



dx. 



(1.34) 



Convergence of the structure function limits the value of 7 to < 7 < 2 (see 
e.g. page 90 of |Monin and Ya glom ( 1975[ >), in which case the integral can be 



evaluated analytically to give (for a derivation, see appendix A I 

Av = —, t^—, nrr , (0 < 7 < 2) . 



(1.35) 



7 r( 7 + 1) sin(7 7r/2) 
Here, is the gamma function. The conclusion to be drawn from this analysis 
is that if the spectrum exhibits a power-law behavior E\\{k\) = £>|&;i|^ 7+1 ), then 
there is a corresponding power-law behavior -Dn(r) = B A y |r| 7 for the structure 
function. 



Kolmogorov ( |1941a ) and Obukhov ( 1941a|b[ ) made the first attempt at deter- 



mining the spectral power-law exponent by assuming that the energy spectrum 
En(ki) has a universal form uniquely determined by the energy dissipation rate 
e. A straightforward dimensional analysis yields 



(1.36) 
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where C k is a dimensionless universal constant. It is related to the Kolmogorov 
constant, C 2 , through (equation 21.25 of Monin and Yaglom ( 1975[ )) 



9i 

Ck 



A 2/3 = fr(i)« 



(1.37) 



Just as in the case with the transverse structure functions, the transverse spectra 
are also related to E n (ki) through 

E 22 (h) = E 33 (h) = l - (l - h E u (h) = ^C k e 2 / 3 A;" 573 . (1.38) 

Equations 1.36 and 1.38 are the famous one-dimensional Kolmogorov — | spec- 
tra. 



1.2 Verification of Kolmogorov theory 

As a preliminary to the investigation of their angular dependence in chapter |5j 
we review, in the following, efforts made in establishing the exact value of the 
scaling exponents, the Kolmogorov constants, C 2 and C k , and their Reynolds 
number trends. 



Work on verification of the two-thirds law, equations 1.19 and 1.20 as well 



as estimation of C 2 , up to the early 1970s has been reviewed by Monin and Ya 



glom ( 1975[ ) (see page 461). Kolmogorov made the first attempt to compare his 



theoretical prediction for the structure functions to the data obtained by Dryden 
et al.| ( 1937| ) with hot-wire anemometry, in which Taylor's frozen- turbulence hy- 
pothesis ( Taylor[ 1938) has been invoked to surrogate temporal velocity fluc- 
tuations for spatial velocity fluctuations (see the sidebar Taylor's Hypothe- 
sis). The measurements were obtained from a wind tunnel grid turbulence, 
whose Reynolds numbers (R\) we estimated to be in the range from 15 to 140. 
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TAYLOR'S HYPOTHESIS: When studying a multiple-scale problem like tur- 
bulence, it is difficult to obtain an instantaneous snapshot of all the eddies 
formed within the turbulence. It is easier, and cheaper, to make measure- 
ments at one point over a long period of time. To study the turbulence from 
a continuous record of measurements from a single point, we need to as- 
sume that the advective velocity, (U), is uniform and the turbulence inten- 
sity, u'/(U),is low. In other words, the turbulence is frozen, or evolves on a 
time scale much larger than the advective one, so that the conversion from 
time to space, x = (U) t, can be made. For a reference on Taylor's hypothesis, 
see 



Frisch(1995), page 58. 



Kolmogorov suggested that C7 2 = 3/2 (Kolmogorov 1941b[ ). (See Kolmogorov 



( |1991 ) for an English translation of this work.) Shortly thereafter, a study car- 
ried out in the atmospheric surface layer by means of hot-wire anemometry 
was undertaken by Obukhov ( 1942[ ), but the study was halted by World War 
II. Obukhov later resumed the study using the hot wire anemometer data of 
Godecke ( Obukhov] 1949a Godecke} |1935 ). Around the same time, many at- 
mospheric surface layer studies employing hot wire anemometry were initiated 



in different parts of the world; for example Townsend ( 1948 ) in the UK; Mac- 



Cready| (1953) and |Cramer| (1959) in the USA; |Shiotani| (1955) in Japan; and [Tay- 



lor | ( 1955) in Australia; after Kolmogorov's work received significant exposure 



due to |Batchelor | ( |1946b[ [1947) . The 1970s saw the emergence of precise hot-wire 



measurements. Measurements of Van Atta and Chen (1970) in the atmospheric 
boundary layer over the ocean showed that both r 2 / 3 and r° 722 power laws fit 
well to their second order structure functions data, and that the value of C7 2 is 
2.3 (error not reported). Comte-Bellot and Corrsin ( 1971) measured longitudinal 
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and transverse velocity autocorrelations but did not report on measurements of 
the structure functions. Anselmet et al. (1984]) investigated structure functions 
up to the 18th order on the axis of turbulent jets, at R\ ranging from 536 to 
852, and in a turbulent duct flow, at R\ = 515. They found that C 2 = 2.2 and 
Dn(r) ~ r ° 71 in the inertial range. An experimental study by Saddoughi and 



VeeravaflT| ( 1994 ), which later became a benchmark experiment in the field, was 



performed on the test-section ceiling of the 80 x 120 foot wind tunnel at NASA 
Ames Research Center. The flows are boundary layer flows at R\ ranging from 
500 to 1450. The Reynolds numbers achieved in this experiment are superior 
to those in many laboratory flows. It has only been surpassed by a few flows 
generated in high-Reynolds-number facilities, including the ONERA wind tun- 
nel with R\ = 3374 (e.g. Gledzer et al. 1996[ ), the Princeton University turbulent 
pipe (e.g. Zagarola and Smits 1998| ) with a R\ of 2.3 x 10 4 (estimated from equa- 
tion 1.7 where Re L is based on the average velocity and the pipe diameter), 



and the wind tunnel of the Central Aerohydrodynamic Institute in Moscow, 
Russia, with a R\ of 3200 ( Praskovsky and Oncley[ 1994[ ). Saddoughi and Veer- 



avalli (1994) observed approximately one decade of r 2 / 3 scaling (reproduced in 



figure 1.2). The value of the Kolmogorov constant reported is 2.0 ± 0.1. Sub- 
sequent measurements of the Kolmogorov constant by Anfossi et al. ( 2000[ ) in 
the atmospheric surface layer (C 2 = 1.7), |Degrazia et alT ( 2008[ ) and Welter et al 



( |2009 ) in the wind tunnel and atmospheric boundary layers (C 2 = 2.47 ± 0.30) 
at R\ in the range between 373 and 2 x 10 4 are marginally consistent with the 
value reported by jSaddoughi and Veeravalli ( 1994[ ). The discrepancy between 
the values of the Kolmogorov constant measured in the laboratory and in the at- 
mosphere may arise due to the uncertainty in atmospheric flow conditions, for 
example variability of wind speed and direction. We note that measurements 
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from the controlled laboratory studies in Degrazia et al. ( 2008[ ) and Welter et al. 
(2009) give C2 = 2.34 ± 0.26, in closer agreement with the value reported by 
Saddoughi and Veeravalli ( 1994[ ). 



Many of the past studies however, have been devoted to verifying the spec- 



tra (equations 1.36 and 1.38). Experimental studies prior to the early 1970s have 
been reviewed by Monin and Yaglom ( 1975[ ) (see page 467), and subsequent 



studies up to the early 1990s have been extensively reviewed by Sreenivasan 
( |1995 ). Data from numerical simulations have been collected by Yeung and 
Zhou ( |1997[ ). We review in the following a few milestone experiments and sup- 
plement the list with several recent ones. 

In a pioneering tidal-channel experiment at R\ ranging from 3000 to 18000, 
Grant et al. (1962^ observed more than two decades of —5/3 scaling and mea- 
sured a value of 0.47 ± 0.02 for C& from data obtained with hot-film anemom- 
etry, which uses the same principle as hot-wire anemometry but the probe is 
less likely to be contaminated by dirt in water The hot-wire measurements of 



Saddoughi and Veeravalli (1994) show slightly less than a decade of —5/3 scal- 
ing, and in that range, Ck = 0.49. The reduced scaling range in this study might 
be an effect of the Reynolds number; we note that their Reynolds number is an 



order of magnitude smaller than that of Grant et al. (1962). The data of over 



100 spectra, mostly collected with hot-wire technique, from various geophys- 
ical and laboratory flows collected in Sreenivasan ( 1995[ ) support a universal 
value of C k , approximately 0.53 ± 0.055, beyond R\ of about 50, as can also be 
inferred from recent experiments of Welter et al. ( 2009[ ). Similar conclusions can 



4 These estimates for the Reynolds numbers were obtained by Sreenivasan (1995 1. See the 
foo tnote therein for d etails . 

^Kraichnan ( 1966 1 has analyzed the same data and observed that the estimate would be 



higher by more than 10% if one looked, instead, for flat region in plots of e- 2 /3 fc5/ 3 Enih) 

VS fci. 



15 



I 
I 




\ t i h i h i| — ^ 



I 



K 

rr. 



■inn i n) — ri < mm — i i i mill — i i i 1 nnj r-ri i mi 



C 2 =2.0 



C' 2 =(4/3) C 2 




] i i inm | — t i ii ni i | — i i i mn | — i i n i nr | — y 



— 

(c) : 




Figure 1.2: The second-order longitudinal (Dn(r)) and transverse (D 2 2(r) and 
-D33( r )) velocity structure functions measured in a turbulent boundary layer 
with R\ = 1450. When compensated appropriately, the structure functions 
within the inertial range of scales are a measure of the value of the Kolmogorov 
constant. Reprinted with permission from Saddoughi and Veeravalli ( 1994[ ): 
/. Fluid Mech., 268:333-372, Copyright (1994) with permission from Cambridge 
University Press. 



16 



be drawn from numerical simulation studies (e.g. Yeung and Zhou 1997 Go- 
toh et al. 2002 Kaneda et al. 2003[ ). On the other hand, these findings seem 
to be at odds with the outcome uncovered in a few experiments, which sug- 
gest Ck follows a power law of the form Ck oc R^- The value of the exponent, 
/i, is 0.094 ± 0.004 in the atmospheric surface layer ( jPraskovsky and Oncley| 
1994), 2/3 in a wind tunnel with specially designed active grid (Mydlarski and 
Warhaft 1998[ ), and 0.16 in a hydrodynamic simulation forced by Taylor-Green 
vortex (Mininni et al. 2008); whereas theory predicts Ck oc (In RxY 1 (Barenblatt 
and GoldenfiId][T995l ). 



These conflicts seem to have been resolved in a recent study by pDonzis and 
Sreenivasan ( |2010 ). The apparent Reynolds-number dependence is ascribable 
to inconsistent procedures followed in the estimation of Ck- A typical proce- 
dure for estimating Ck is to plot the compensated spectrum, e -2 / 3 k^ 3 Eu(ki), 
and seek a plateau, whose height in the inertial range is taken as Ck- This proce- 
dure is, however, complicated by the absence of a distinct plateau in the inertial 
range, and the existence of a spectral bump at the high-wavenumber end (the 
near-dissipation range), known as the bottleneck, whose effect is most promi- 
nent in the three-dimensional spectrum. The occurrence of the spectral bump 
was first captured in an analysis of |Qian ( 1984[ ). The physical mechanism lead- 



ing to its appearance has been attempted by various workers (e.g. Yakhot and 



Zakharov, 1993, Falkovich, 1994; Martinez et al- 


1997, 1 


<urienetal. 2004,Verma 




and Donzis , 2007, Bershadskii 


2008;Frischet al. 


2008) 


Donzis and Sreenivasan 



(2010) showed that the height of this bump in the three-dimensional spectrum 
varies as R^° M , while the height and the location of the local minimum pre- 
ceding the bump remain practically constant. Donzis and Sreenivasan ( 2010[ ) 
suggest evaluating Ck from the height of the local minimum before the bump, 
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as this leads to a Kolmogorov constant independent of R\. It remains for future 
investigations to elucidate the underlying mechanism of this operation. 

Since the spectra are uniquely related to the structure functions, one may 
expect that the bottleneck effect in wavenumber space will have its physical 
space correspondence near the transition between the dissipation and the iner- 
tial ranges. Dobler et al. ( 2003[ ) argued that localized features in wavenumber 



space will become sufficiently non-local in physical space that they are practi- 
cally undetectable. Lohse and Miiller-Groeling ([1995), on the other hand, ar- 



gued that a rapid transition between these ranges in the physical space is the 



manifestation of the spectral bump. The findings of Donzis and Sreenivasan 



(2010) confirm this belief. To assess the effect of the Reynolds number on the 



transition region, they studied the Batchelor interpolation formula (Batchelor 
1951) 



D ( r ) _ K ( r A?) fl39 v 

with K = 1/15 for the longitudinal and 2/15 for the transverse structure func- 
tions, u v is the Kolmogorov velocity (see equation 1.11[ ), and ( 2 = 0.67 is the iner- 



tial range scaling exponent. c B and q are fit parameters. They found c B = 0.076 
and 0.102 for longitudinal and transverse structure functions, independent of 
R\, and the power laws q oc R^ 1 06 and q oc R^ 92 for the longitudinal and 
transverse structure functions, respectively. The decrease of q with R\ confirms 
the belief that the transition at high Reynolds numbers is smoother, which is 
the manifestation of a flatter spectral bump. Such a subtle Reynolds number 
dependence is currently beyond the resolution limit of experiments and may 
have escaped definitive detection. Notwithstanding this subtle effect, at least 
the contamination by a 'physical' bump in the determination of C 2 can be obvi- 
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ated in the physical space^j Between the approach in spectral space and that in 
physical space, the approach in physical space remains the better one. 



1.3 Relation to other inertial-range constants 

The following discussion is intended to illustrate the interrelations between 
the Kolmogorov constant and the scaling constants of the Lagrangian structure 
functions, the Richards on-Obukhov constant for turbulent relative dispersion, 
and the Obukhov-Corrsin constant for passive scalar structure function. The 
scalar-valuedness of the Kolmogorov constant to be established in chapter [5] 
may be taken as a constraint on the angular dependence of the various con- 
stants. 

It is often the case that analyzing the same physical problem from a differ- 
ent perspective may reveal unexpected results and advance our understanding 
of the problem. The structure functions introduced in the previous section are 
functions of separations fixed in space. Analyzing fluid motion at fixed spatial 
locations is known as the Eulerian description. To inspect the turbulence prob- 
lem in a different way, researchers have devised the Lagrangian description, 



where an observer moves along with the flow (e.g. Toschi and Bodenschatz 



2009 j ). The basic tools in the Lagrangian description of turbulence are the La- 



grangian structure functions, defined as 

Di(r) = (8u i (T)5u j (r)), (1.40) 
where the velocity increments, 5ui(r) = u[{t + r) — are taken along the 



6 The typical procedure in physical space is to plot, for example, (e r) 2 / 3 Dn(r) vs r and seek 
a plateau whose height is taken as the value of C2 ■ 
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trajectory of an individual fluid particle. 



Obukhov and Landau applied Kolmogorov's theory to derive scaling 
forms for the Lagrangian velocity structure functions valid in the inertial range 

(S Ui (r) 5 Uj {r)) = C er% , (1.41) 

where i is any Cartesian component x, y, or z; e is the energy dissipation rate 



per unit mass; and Sij is the Kronecker delta. See e.g. Monin and Yaglom 



( 11975 ), page 358 for a derivation. The constant C , presumed to be univer- 



sal, is known as the Lagrangian velocity structure function constant. It is an 
important parameter in stochastic models of turbulent transport and disper- 



Rodean 


1991 


Sawford 


1991 



measurements of C are scarce and its value is very uncertain, partly because 
Lagrangian experiments, where the trajectories of fluid particles are followed 



in both time and space (e.g. La Porta et al. 2000), have historically been very 
difficult. Reported values of C range from 1.0 to 7.0 ( [Poggi et aT| |2008| ). [In 
( 1951[ ) applied the Kolmogorov theory to the Lagrangian velocity spectrum 



oue 



(f)fj(u>), defined as the Fourier transform of the Lagrangian velocity autocorrela- 
tion tensor, and obtained <j>%j{oj) = B eu~ 2 6ij, where u is the angular frequency. 



B and C are simply related by C /B = A\ = n, see e.g. Monin and Yaglom 
( |1975 ), page 361. Under the assumption that the velocity fluctuations integrated 
from Eulerian and Lagrangian spectra are equal ( Lumleyj 1957, 1962), Corrsin 
( |1963a ) approximated the Eulerian and Lagrangian spectra by their inertial 
range forms between appropriate wave number and frequency limits and ob 



tained B/C\ 



3/2 



constant, where B is the scaling constant of the Lagrangian 



velocity spectrum. The same power law has recently been derived by Franzese 



7 The scaling form for the Lagrangian velocity structure function was independently obtained 
by Obukhov and Landau shortly after the publication of Kolmogorov's 1941 paper. See|Monin 



and Yaglom i 1975|, page 359, for a historical account. 
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and Cassiani| ( |2007j ) based on a statistical diffusion theory of relative dispersion 
of fluid particles. 

The Lagrangian velocity structure function constant is, quite astonishingly, 
also related to the Richards on-Obukhov constant governing the turbulent rela- 
tive dispersion from a point source in isotropic turbulence. Obukhov ( 1941b[ ), 
refining the empirical finding of Richardson ( 1926[ ), applied Kolmogorov's di- 
mensional reasoning to the mean square separation distance between two par- 
ticles (r 2 (t)) and derived an expression valid in the inertial range of scales 

(r 2 (t))=get\ (1.42) 

where e is the energy dissipation rate. The scaling constant g is known as the 



Richardson-Obukhov constant. See e.g. Ouellette et al. (2006) for a modern 
derivation. Experimental and numerical studies seem to be converging on a 
value in the range 0.5 < g ^ 0.6 ( |Salazar and Collins} |2009~| ). The interest in 
the law of two-particle dispersion is immense because of its utility in describing 
transport and mixing processes in environmental and engineering problems. 



Lin 


1960; 


Novikov 


1963 


Ivanov and Stratonovich 



1963[ ) exploited the assumption that the accelerations of fluid particles are un- 
correlated at high Reynolds number to connect two-particle statistics with La- 



grangian one-particle statistics. All their analyses yield the result g = 2 Co. Bor 
gas and Sawford ( 1991[ ) showed that taking into account the two-particle accel- 
eration correlation provides a value of g smaller than 2 C . The recent work of 
Franzese and Cassiani ( 2007[ ) based on a statistical diffusion theory of relative 
dispersion provides an estimate g w Co/11. 

The relative motion between particles is, remarkably, connected to the con- 
centration profile of a substance, the so-called passive scalar. The substance 
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could be a contaminant, like smoke dispersing in air; it could also be odor, like 
the smell of perfume traveling in air; or heat, as in the case when a weakly 
heated object is cooled in the flow. In his seminal paper on turbulent relative 
dispersion, Richardson ( 1926[ ) introduced the distance-neighbor function q(r, t), 
which is the equal-time probability density function for the separation r be- 
tween fluid particle pairs randomly chosen from a scalar field 6'(r, t) passively 
advected by the turbulence. They are related by (e.g. Ott and Mann[ 2000) 



q(r,t)= (e'(x,t)9 , (x + r,t))dx. 



(1.43) 



Assuming homogeneity, isotropy, and stationarity, we also have by definition 

POO 

(r 2 (t)) = / A7rs A q(s,t)ds. (1.44) 
Jo 



Equations 1.43 and 1.44 provide a formal link between the mean square sep- 
aration distance between particle pairs and the autocorrelation of the passive 



scalar fluctuations. Assuming the classical inertial range scaling for q, Thorn 



son 



( 1996) showed that the connection between the scalar dissipation rate and 
relative dispersion is given by 

d(9' 2 ) 3 /fl/2 . d\n(r 2 ) 



> 



(1.45) 



dt 2 x ' dt 

Obukhov (1949b]) and Corrsin (1951) independently extended Kolmogorov's 
theory to the passive scalar and derived for the second order structure function 
of scalar increments 



{(e'(r,t)-6'(0,t)) 2 ) = C e eee- l / 3 r 2 / 3 



(1.46) 



where e is the energy dissipation rate, e e = 2 k ((W) 2 ) is the rate of dissi- 
pation of concentration fluctuations, and k is the diffusivity. Ce is known as 
the Obukhov-Corrsin constant. Experimental values of Ce lie mostly in the 



22 



range 1.21 < C e < 2.01 (Sreenivasan 1996[ ). Thomson ( 19961 ), applying the 



two-particle statistical theory of Batchelor ( 1952 ), obtained a result relating the 
Obukhov-Corrsin constant to the Richards on-Obukhov constant 



g — 243 a l o e 



(1.47) 



for some constant «i that shapes the spreading of the Richardson ( 1926[ ) 
distance-neighbor function. 

While many of the inertial range constants still elude precise theoretical 
modeling and experimental quantification, we believe that a complete theory 
of turbulence should and must unveil the interrelations between them. Here, 
we make no attempt at seeking such a fundamental relation, but we verify in 
chapter [5] the scalar invariance of C%, which is inherently implied in all theoreti- 
cal investigations that invoked the local isotropy hypothesis, as well as in many 
single-point measurements that assumed Taylor's hypothesis. 



1.4 Higher-order structure functions 

In preparation for the examination of the influence of anisotropy on the scaling 
exponents of structure functions of order greater than 2 in chapter [6j we intro- 
duce in the following the higher-order structure functions and their role in the 
study of intermittency. 



According to Landau and Lifshitz ( 1959[ ), the arguments leading to Kol- 



mogorov's prediction presented in the preceding sections do not take into ac- 
count the possible 'bursty' nature of the energy dissipation rate. A distinctive 
feature of turbulence is its abrupt and extremely intense fluctuations in the ob- 
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servables. Landau's main objection is that the e in Kolmogorov's theory is a 
mean taken over time. He argued that since e fluctuates in time, (e n ) ^ (e) n , 
except when n = 1, and therefore the mean rate of energy (e) is insufficient to 



describe turbulence. This phenomenon is known as intermittency. (See Frisch 



( 1995 ), 6.4, for a historical account of Landau's objection.) 



It has now become common in studies of intermittency to look at higher- 
order velocity structure functions. Kolmogorov's arguments can be naturally 
extended to moments of velocity increments, Su(r) = u'(x + r) — u\x), to order 
higher than two. For example, the pth-order longitudinal structure function, 
Dp(r) = ((SuL(r))' p ), where the subscript L indicates that the direction of veloc- 
ity is along the separation vector, scales in the inertial range as 

D L p (r)=C p (erf^ (1.48) 

with universal scaling exponents given by the Kolmogorov prediction 

Ci = p/3 • (1-49) 

The values of Of given by the above equation for order 1 to 6 are listed in ta- 



ble 1.1 The constants C p s are presumed to be universal. To date, there ex- 
ist no measurements that display unambiguous inertial range scaling over an 
enormously wide scaling range ( jSreenivasan and Dhruva 1998[ ). Thus, several 



self-consistent procedures have been devised for better estimation of the scal- 



ing exponents. The Batchelor interpolation formula (equation 1.39), based on 
matched asymptotics, is useful but involves a very elaborate nonlinear fitting 
procedure (see e.g. Kurien and Sreenivasan 2000[ ). Another simple and intrigu- 



ing procedure has been proposed by Benzi et al. ( 1993[ ). They observed that 



structure functions plotted against any other structure function of a given or- 
der, rather than against the separation, generally yields a measurably improved 
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inertial scaling range at low Reynolds numbers. This procedure is known as 
the extended-self-similarity (ESS) method. The reason for its success in Navier- 
Stokes turbulence is unclear, but recent progress in Burgers turbulence explains 
this transformation of variables from physical to structure function space as a 
method to deplete the subdominant inertial and dissipative range contributions 
that mask power-law scaling ( Chakraborty et aE| 2010[ ). 



In practice, the third-order longitudinal structure function, ((SuL(r)) 3 ), is 
chosen as the surrogate for separation. The foundation for this choice lies in 



an exact result obtained by Kolmogorov ( 1941b ) from the Navier-Stokes equa- 



tion for freely decaying isotropic turbulence, and it remains the only exact re- 
sult derived from the equations of motion. Starting with the Karman-Harwarth 



equation (Karman and Howarth 1938), Kolmogorov rigorously showed that 



er + 61/ 



(1.50) 



5 dr 

In the inertial range, this reduces to the celebrated Kolmogorov's four-fifths law 
( jKolmogorov 1941b[ ) 



D L ,{r) 



er . 



Thus, relating the above to equation 1.48 the results Q 



(1.51) 

-4/5 and C 3 L = 1 



are both exact and nontrivial. To be precise, ESS can only be used to determine 
relative scaling exponents. Since (3 is exactly 1, the relative scaling exponents 
may be directly related to the true scaling exponents. 

Another variant of ESS, which gives even better scaling, modifies structure 
functions of order p by taking the moments of the absolute values of the ve- 



locity increments, (\5u L \' p ), the so-called generalized structure functions (Vain- 



shtein et al. 1994| ). Strictly speaking, the four-fifths law only applies to ((Sul) 3 }, 



with no absolute value. Sreenivasan et al. (1996) have found small differences 
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in the scaling of ((Su) p ) and (\Su\ p ). The difference is small, however, and is 
very difficult to observe experimentally. One then replaces the scaling laws in 



equation 1.48 by 



(\6u L \*) oc (\6u L \ 3 )$ . (1.52) 

Plotting the structure functions relative to (\Sul\ 3 ) shows cleaner inertial range 
scaling behavior than plotting them against the separations. Because of its util- 
ity in the longitudinal structure functions, the ESS technique has been extended 
to the transverse case. The transverse ESS scaling laws 



(\8u T \ p ) ex (\5u T \ 3 }^ , 



(1.53) 



where the subscript T denotes velocity direction perpendicular to the separa- 
tion vector, and Cj = p/3 is the inertial range scaling exponent given by the 
Kolmogorov prediction, Z?J(r) = C p {er) c i '. This version of ESS is less justi- 



fied, but is widely in use (e.g. Herweijer and Van de Water 1995 : Camussi et al. 



1996 


Kahalerras et al. 


1996J Noullez et al. 


1997 


Van de Water and Herweijer 




1999, 


Zhou and Antonia 2000; Pearson and Antonia 2001 Zhou et al.| 2001), 



although defining the transverse structure functions in this way always yields 
values for the odd moments, even when they do not exist for homogeneous 
isotropic turbulence (Shen and Warhaft 2002[ ). We will see in chapter [5] that at 



moderate Reynolds numbers, the inertial range is too narrow to fit a power law 
in r, and we will make use of ESS to determine ( p . 



1.5 Anomalous transverse scaling exponents 



As we have seen in the previous section, in order to verify the Kolmogorov 
prediction for the higher-order scaling exponents (equation 1.49), we naturally 
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(p Source 



Order p 



Dhruva et al. 



(1997) 



1 


2 


3 


4 


5 


6 


0.333 


0.667 


1 


1.333 


1.667 


2 


0.364 


0.696 


1 


1.280 


1.538 


1.778 


0.366 


0.700 


1 


1.266 


1.493 


1.692 



Dhruva et al. 



(1997) 



0.359 0.680 0.960 1.200 1.402 1.567 



Table 1.1: The inertial range scaling exponents of the longitudinal and trans- 
verse structure functions. The first two rows are the longitudinal sc aling ex- 
ponents, C^, pr edicted by Kolmogorov ( 1941a[ ) (see equation 1.49) and She and 
Leveque ( 1994b (see equation 1 .54 ( . The third and fourth rows show the longi- 
tudinal and transverse exponents measured by Dhruva et al. ( 1997[ ) using the 
ESS method (see section 1~4~| ) in atmospheric turbulence at R\ between 10 4 and 
1.5 x 10 4 . 



need to measure the higher-order structure functions. Measurements of the Ion 
gitudinal structure functions up to the fourth order were first obtained by |Van 
Atta and Chen ( 19701 ), and subsequently, up to the 18th order, by Anselmet et al 



(1984). The longitudinal scaling exponents, (r, from both measurements show 



departure from the expression 1.49 (the difference is less than 5% at fourth order, 
and increases monotonically up to 30% at 12th order), a phenomenon known as 



anomalous scaling, but it is well predicted by the She-Leveque model (She and 



Leveque) , 1994[ ) using a hierarchical model for the energy dissipation rate, which 
gives for the scaling exponent 



Sp 



P 



+ 2 



1 - 



P/3l 



(1.54) 



The values of Of given by the above equation for p up to the sixth order are 
listed in table ll.ll 



It was thought that the transverse scaling exponent, (E (defined in equa- 
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tion 



1.53), is equal to (£. This is true only for p = 2 because of the incompress- 



ibility condition (see equation 1.18[ >. For p ^ 2, the transverse and longitudinal 
scaling exponents are not related by any constraint. In fact, measurements of 



Herweijer and Van de Water (1995) in shear flows at R\ = 340 — 810 did sug- 



gest that the difference Cj — p/3 is larger than Cp — p/3 for p ^ 4, implying that 
Cp < Cp f° r V ^ 4. The difference persists even at R\ up to 10 4 , as shown by 
the hot-wire measurements of Dhruva et al. ( 1997[ ) in the atmospheric surface 



layer, reproduced in table 1.1 The same result is found in other experiments 



( |Camussi and Benzi[ 1997 Antonia and Pearson 1999 Zhou and Antonia 2000 
Romano and Antonia| |2001f |Hao et aL| [2008) and numerical simulations ( |Chen 



eTaTl[T997t[Boratav and PeTz|[l997][Grossmann et"aT|[T997]|Gotoh et al4|2002"l ). A 



number of experiments, however, suggest that the two exponents are equal (Ca 
mussi et al4[l"996t|Noullez et al4|l997}|Kahalerras et al.||l998}|He et al4|l999l ). A 



few authors (Pearson and Antonia 2001 Zhou et al. 2001 : Antonia et al. 2002 



Shen and Warhaft 



2002 



Zhou et al. , 2005 1 report both Cj ~ Cp and Cj < Cp in 



their experiments. 



Table 1.2 is an assessment of the issue. We would like to caution the reader 



of a too simplistic interpretation of table 1.2 The flow conditions, measure- 



ment techniques, and the definitions of Cp and Cj vary from experiment to ex- 
periment. We note that in most experiments, Taylor's hypothesis has been in- 
voked to surrogate temporal series of measurements for spatial ones (see the 
sidebar Taylor's Hypothesis in section 1.2). On the other hand, Lin ( 1953[ ) has 
shown that the hypothesis breaks down in the case of flows containing high 
shear. Correction to Taylor's hypothesis based on spectral information gathered 
in only one direction and on a local derivative in the remaining direction has 



been proposed by del Alamo and Jimenez (2009). We note that in some cases 
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Flow TH R\ Scaling exponent Source 



Shearless 


No 


100 - 300 


C T < C L 

Sp pj Sp 


Shearless 


No 


863 


C T ~ C L 

Sp Sp 


Shearless 


Yes 


100 - 300 


Sp Sp 


Shearless 


Yes 


863 


Sp Sp 


Shearless 


Yes 


10 4 - 1.5 x 10 4 


C T < C L 

Sp r^j Sp 


Sheared 


No 


254 


C T < C L 

Sp ^ Sp 


Sheared 


No 


875 


C T ~ C L 

Sp Sp 


Sheared 


Yes 


254 


ESS not applicable 


Sheared 


Yes 


875 


Sp Sp 



Zhou et al. 


(20051 


Shen and Warhaf t 


2002|) 


Zhou et al. 


(20051 


Shen and Warhaf t 


2002|) 


Dhruva et al. (1997 





Shen and Warhaft 



Shen and Warhaft 



Shen and Warhaft 



( |2002| ) 
( |2002| ) 
( |2002| ) 



Shen and Warhaf t| < |2002[ > 



Table 1.2: The equality between longitudinal and transverse scaling exponents 
for sheared and unsheared turbulence and their trends with Reynolds number. 
The acronym TH denotes the use of Taylor's hypothesis in the experiments. 



when the scaling exponents do agree, Taylor's hypothesis has not been invoked 
(e.g. |Noullez et al.Hl997||Shen and Warhaf tH2002)|Zhou et ai]|2005) or a modified 
version based on the instantaneous velocity u(t), instead of the mean velocity 



{U), has been used to obtain the spatial separation Ar = u(t) At (e.g. Kahalerras 



et al. 1998[ ). Because many authors report only structure functions constructed 
from single-point measurements of the velocity signals, we have deliberately 



limited the scope of comparison in table 1.2 to those that report two-point mea- 
surements that truly scan the physical space. 

Our second observation is that there is no consensus among researchers on 
the choice of the third order structure function for the ESS method. While all 
researchers agree on using the ESS scaling (\6u L \ p ) oc (|5n L | 3 )^ for (£, the ESS 
scalings (\5u T \ p ) oc (I^mlI 3 )^ and (\5u T \ p ) oc (|<5wt| 3 )^ have been proposed for 
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Cp. We observe in the data of Pearson and Antonia ( 2001| ); Zhou et al. (2001 



2005) that ESS scaling based on (\5ul\ 3 ) tends to produce a wider scatter and a 



value for (J that is smaller than the one obtained with (\5ut\ 3 )- We think that 
the uncertainty in the reported values of Cj partly stems from the inconsistency 
in the procedure used in obtaining the transverse scaling exponent and our at- 
tempt here is to test the directional dependence of {\Su T \ p ) oc (\5u T \ 3 )^ ■ 

The difference between the longitudinal and transverse scaling exponents 
remains one of the unresolved issues of turbulence. The reasons for this in- 



equality have been proposed by various authors. For example, Romano and 



Antonia ( 2001[ ) and Zhou et al. ( 2001[ ) list the following as plausible explana- 



tions: a) the anisotropy of the large scale, b) the Reynolds number effect, c) the 
differences in initial and boundary conditions, d) different intermittency. e) the 
effect of Taylor's hypothesis. Among these possible causes, which may be re- 
lated, [RomanoandAntonia ( 2001[ ) have singled out the anisotropy of the flow 



as the most likely source for the observed inequality. 

A complementary approach has been developed to extract the anisotropic 
contributions to the scaling of structure functions (Biferale and Procaccia 2005). 



This approach uses the irreducible representation of the SO(3) symmetry group. 
Instead of projecting the structure functions onto the longitudinal and trans- 
verse components in the traditional way, we expand the structure functions in 
terms of spherical harmonics^J, Y™, sidestepping the questions about the dif- 
ferences between longitudinal and transverse scaling exponents. Because the 
isotropic (£ = 0) and anisotropic parts (£ = 1, 2, . . .) of the structure functions are 
naturally decoupled in this representation, one could ask how does the value 



8 Following the usual convention, the orbital angular momentum is represented by I. Its 
projection onto the z axis is represented by m. 
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of the anisotropic scaling exponent compare to that of the isotropic one? The 
first-order anisotropic scaling exponent has been calculated using perturbation 
theory ( |Grossmann et al" 1994 Falkovich and L'vov 1995] ) and is found to be 
higher than the isotropic one. Borrowing results from passive scalar turbulence 
( jFairhall et al. 1996 ), it is then conjectured that all higher order anisotropic scal- 
ing exponents are positive, greater than the isotropic scaling exponent, and in- 
crease with increasing order ( LVov and Procaccia| 1996[ ). If this is true, then 
these nondecreasing anisotropic exponents would neatly explain the diminish- 
ing anisotropic scaling contributions with decreasing scale. Initial tests have 
been conducted in atmospheric boundary layer flows ( Arad et al. 1998} Kurien 



et al. 2000} Kurien and Sreenivasan[ |2000 ) and the results suggest a hierarchy 



of increasingly larger anisotropic scaling exponents with increasing order. Suc- 



ceeding measurements in laboratory homogeneous shear flows (Warhaft and 



Shen 2002[ ), however, have cast doubt on the conclusion reached by the above 



authors and more refined experimental analyses have found the results to de- 



pend on the geometrical configuration of the measurement probes (Staicu et al. 



2003). The SO(3) decomposition is in its infancy. The technique clearly awaits a 



more impartial analysis and more refined experimental techniques^ 



1.6 Anisotropic turbulence 



As discussed in section 1.1 Kolmogorov's theory presumes that, at sufficiently 
high Reynolds numbers, the small scales of the turbulence are statistically ho- 
mogeneous, isotropic, and free of externally imposed boundary. Few geophys- 



9 The well-crafted argument that the method might improve at higher Reynolds numbers is 
problematic because it is precisely at moderate Reynolds numbers that the method makes many 
of its claims (e.g. Arad et al. 1999 Biferale and Toschi 2001} . 
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Flow 



Anisotropy Direction 



Source 



Urban surface layer 
Rural surface layer 
Marine surface layer 
Mixing layer 
Canopy 

Sea surface current 
Martian boundary layer 



1.9 : 1.5 : 1 streamwise/ spanwise/ vertical 

2 : 1.5 : 1 streamwise/ spanwise/ vertical 

1.83 : 2.6 : 1 streamwise/ spanwise/ vertical 

1 . 3 streamwise / vertical 

1 . 7 streamwise / vertical 

0.44 - 2.18 zonal /meridional 

1.06 zonal / meridional 



Roth 



(20001 



Counihan (19751 



Frieheet al. (1991 1 



Finnigan 



(2000) 



Finnigan (2000) 



Ducet et al. 



(2000) 



Sullivan et al. (2000) 



Table 1.3: The landscape of anisotropy for geophysical flows. Anisotropy refers 
to the ratio of streamwise to spanwise to vertical, streamwise to vertical, or zonal 



to meridional RMS velocity fluctuations, depending on the context. In Sullivan 



et al. ( 2000[ ) the RMS velocity fluctuation data are unavailable. The anisotropy 
quoted is for the mean wind speed. 



ical and laboratory flows fulfill these requirements. Field measurements in the 
atmospheric boundary layer show that atmospheric flows attain the highest 
ever Reynolds number (R x ) achievable on Earth, typically of the order 10 4 (e.g. 
Dhruva et al.|[l997| ). Complex terrestrial terrain and planetary dynamics, how- 
ever, may influence the energy-containing scales of turbulent motion. Conse- 



quently, one would expect motion at these scales to be anisotropic. Table 1.3 
offers a few glimpses of the degree of anisotropy in nature. And how does 



man-made turbulence compare with nature? Our reaction to table 1.4 is one of 



amazement. Despite the moderate Reynolds number, laboratory flows created 
under different conditions possess approximately the same degree of anisotropy 
as nature. 



What causes this anisotropy in the energy-containing scales of turbulence? 
Is it a reflection of the asymmetry in the mechanism that agitates the fluid? 
There is ample evidence to support a causal relationship between asymmetry 
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and anisotropy, but the relationship has not been tested experimentally. This 
evidence includes flows where asymmetry and anisotropy are present at the 
same time, such as in turbulent pipe flows (Pearson and Antonia 20011 ), turbu- 
lent jets (e.g. Romano and Antonia 2001[ ), counter-rotating von Karman flows 



(e.g. Voth et al. 2002 ), and wind tunnels with specially designed shear gen- 
erators (e.g. Shen and Warhaft 2000 Isaza et aTj 2009). Turbulence produced 



by computer simulation can also be forced in an asymmetric way, and Yeung 
and Brasseur ( 1991[ ) studied the influence of this asymmetry. One exception is 
the turbulence produced by a grid in a wind tunnel with a specially designed 
contraction, where the axis of the tunnel introduces a clear asymmetry, yet the 
turbulence produced is nearly locally isotropic, but decaying ( |Comte-Bellot and 



Corrsin 1966). 



There is a pattern in the relationship between asymmetry of the forcing and 
anisotropy of the turbulence. Machines with a single axis of symmetry, such 
as the von Karman flow generators (e.g. Voth et al., 2002) and wind tunnels 
(e.g. Comte-Bellot and Corrsin 1966[ ), produced axisymmetric turbulence. Ma- 
chines with more axes, such as the one developed by Hwang and Eaton (2004), 
produced isotropic turbulence. A careful study by Zimmermann et al. ( 2010| ) 
showed that six axes were sufficient to produce turbulence without a preferred 
direction. The machine to be described in this thesis has 16 axes. 
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Batchelor ( 1946a[ ) initiated the study of anisotropic turbulence by consider- 



ing turbulence with reduced symmetry, namely symmetry about an axis, draw- 
ing motivation from wind-tunnel and pipe flows. Applying the invariant the- 
ory of Robertson ( 1940[ >, Batchelor (1946a]) showed that the correlation tensor, 



(w-(r + x) u'Ax)), can be described by a single scalar function in sys- 



tems of high symmetry (isotropic), but in axisymmetric turbulence additional 



functions are required. Chandrasekhar ( 1950 ) showed that it can be written in 
terms of only two independent scalar functions, which are not directly measur- 
able. For two special forms of this tensor, which represent possibly the simplest 
kinds of non-isotropic turbulence, Sreenivasan and Narasimha| ( 1978 ) showed 
that it is sufficient to give two functions of a single variable or even one function 
of a single variable with an additional parameter. An experimentally-oriented 



formulation developed by Lindborg ( 1995 ) re-expressed the correlation tensor 
in terms of 6 measurable two-parameter scalar functions, but because of reflec- 
tional symmetry only two of them are independent. 

As the rank of the tensor increases, so does the number of independent func- 
tions we need to describe the tensor. Simple symmetry arguments, however, 
may reduce the number of terms and greatly simplify the expression. Let us 
consider Rij{xk) where i, j, k are cyclic permutations of the indices 1, 2, 3. Re- 
flectional symmetries of the Navier-Stokes equation about the x — y, y — z, and 
x — z planes yield (e.g. Lindborg 1995[ > 



Rij{xk) = o 



(1.55) 



for all permutations of (i, j, k) = (1, 2, 3). This is nothing but a consequence of 
the conservation of parity. 

In principle, to fully describe the axisymmetric turbulence, all 6 functions 
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would need to be measured. Here, we ingeniously designed our flow appara- 
tus to allow rotations of the anisotropic large scale forcing, instead of the mea- 
surement apparatus, in order to measure parts of two out of these six functions 
with separations that extended radially from the center of the turbulent region. 
In the same spirit as the previous sections, we shall examine in chapter [7] the 
directional dependence of these correlations in anisotropic turbulence. 

1.7 Overview 

The central idea that interweaves all the chapters in this thesis is anisotropic 
turbulence. From this idea, we develop three interrelated sub-themes. First, we 
discuss the generation of anisotropic turbulence in the laboratory. In chapter |2j 
we introduce the turbulence generation apparatus and describe the accompany- 
ing flow measurement technique. Following this, in chapter |3j we explain the 
coordinate system we have used to present our measurements, and describe the 
measurement protocol and methods. In chapter]?} we examine the performance 
of the turbulence generator. We discuss three aspects of the large scale velocity 
fluctuations: the uniformity in the anisotropy, axisymmetry, and homogeneity. 
Combining these and the auxiliary measurements of the mean flow, the veloc- 
ity fluctuations mixed correlation, and the Reynolds stress, we evaluate to what 
degree the turbulence in the fluid is homogeneous and axisymmetric. 

The second idea concerns the influence of large scale anisotropy on the in- 
ertial scales of turbulence. This influence is traditionally studied through the 
structure functions. In chapter |5j we review the existing works on anisotropy 
in structure functions. We examine an aspect that has been overlooked in the 
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literature, namely the behavior of second-order transverse structure functions 
in different directions of the flow. We found that the two structure functions be- 
come equal at the inertial scales. This contrasts with existing results obtained in 
shear flows and suggests that anisotropy produced in the fluctuations and that 
produced by shear are likely to be fundamentally different. Our measurements 
suggest that the Kolmogorov constant is independent of the direction in which 
it is measured. This we exploit to determine the parameters of the turbulence. 
Following this, in chapter [6j we examine structure functions of order four, five, 
and six. We illustrate the difficulty in obtaining precise measurements of these 
structure functions with measurements of their probability density functions. 
We demonstrate the possible influence anisotropy might have on the shape of 
the probability density functions. 

The third idea covers the influence of large scale anisotropy on the integral 
scales of turbulence. This is studied through the correlation functions. In chap- 
ter |7j we devote ourselves to the analysis of transverse correlation functions 
measured in different directions. We observed that they deviate from each other 
at the large scale. We propose a scaling form similar to one that is studied in 
problems dealing with critical phenomena, and combine it with Kolmogorov- 
type inertial range scaling to show that the data collapse when appropriately 
scaled. We then derive a power-law relationship between the integral length 
and the large scale velocity fluctuation and demonstrate that the power-law ex- 
ponent is indeed linked to the inertial range scaling exponent. 

Finally, in chapter |8j we summarize the results presented in this thesis and 
discuss several shortcomings of our arguments, as well as future measurements 
that may complement and extend the results of this thesis. 
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CHAPTER 2 

THE FLOW APPARATUS AND THE MEASUREMENT TECHNIQUE 



The goal of the present experimental study is to generate a closer approxi- 
mation to tunable and anisotropic, but unsheared, turbulent flows than those 
obtained by |Shen and Warhaft ( 20021 ), which we accomplished through the con- 



struction of a flow apparatus described hereunder. The apparatus makes sys- 
tematic studies of anisotropic turbulence feasible, and a direct comparison of 
anisotropic turbulence with isotropic turbulence in a single apparatus possible. 
The flows generated by this apparatus are characterized using laser Doppler ve- 
locimetry (LDV). We review the LDV measurement technique and discuss some 
of the essential features, only to a degree sufficient for the understanding of this 
thesis. 

2.1 The flow apparatus 

We now describe our turbulence chamber, the loudspeaker turbulence genera- 
tor and the loudspeaker amplifier unit. We review in some detail the algorithm 
we used for driving the turbulence generators. 



2.1.1 The chamber 

The geometry of our flow chamber was inspired by the apparatus developed by 



Zimmermann et al. ( 2010[ ) at the Max Planck Institute for Dynamics and Self- 



Organization in Gottingen, Germany. As shown in figure |2.1[ our turbulence 
chamber had the shape of a truncated icosahedron and it had a diameter of 
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1m 




Figure 2.1: The 'soccer ball' chamber. Also shown is the layout of the Laser 
Doppler velocimetry probes in the experiment. A two-component probe 
aligned along the a^-axis measured the components of velocity in the x 2 and 
x 3 directions, and an additional one-component probe aligned along the x 2 axis 
measured velocities in the £3 direction. 
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Figure 2.2: An overview of the experiment showing the 'soccer ball', the ampli- 
fier unit, and the measurement apparatus. The amplifier received signals from 
the multichannel high-speed voltage output device on the computer through a 
BNC connector block, and sent the signals to the loudspeaker array on the soc- 
cer ball. A laser Doppler velocimetry (LDV) system collected velocity signals 
and saved the data on a computer (not shown). 
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99 cm. It was made of wood with twelve regular pentagonal and twenty regular 
hexagonal faces joined together with nylon straps and glue. The inner and outer 
surfaces were coated with several layers of lacquer to make the chamber water 
repellant. In the center of each face, a circular hole was cut for the jet generator. 
In addition, further circular holes adjacent to the jet generators were cut out for 
optical access. As the shape resembles that of a soccer ball, we also refer to the 
flow apparatus as the 'soccer ball' chamber. 



2.1.2 The loudspeaker 



As depicted in figure 2.3 the turbulence mixers were jets, or jet-like winds that 
are produced by powerful acoustic sources in air, developed for flow control 
applications ( jGlezer and Amitay[ 2002[ ). A 150 W loudspeaker with a diame- 



ter of 16.5 cm and a uniform frequency response between 50 and 3000 Hz, was 
mounted on each face of the soccer ball and pointed towards the center of the 
chamber. The loudspeakers pushed air in and out through conical nozzles of 
opening angle 30°, length 4.3 cm and orifice diameter 5 cm. Each nozzle was 
held between a face of the soccer ball and a square wooden plate. By using 
a sinusoidal driving of the loudspeaker, we generated a pulsating turbulent 
jet, a phenomena known as acoustic streaming ( LighthilTj |1978 ). When back- 



ing away from the orifice, the loudspeaker diaphragm ingests air from all di- 
rections, while it blows air out in a single direction when moving toward the 
orifice. As the jet travels downstream and reaches the center of the soccer ball, it 
interacts with other jets created by other loudspeakers. These interactions pro- 
duce a turbulent region that is more intense than the turbulence due to a single 
turbulent jet. The sound level inside the soccer ball was typically 135 dB, which 
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Figure 2.3: A loudspeaker with a conical nozzle. 



corresponds to 0.2% of the measured turbulent kinetic energy. Thus, we expect 
the sound to have a negligible effect on the turbulence and the measurements 
of velocity. 



The flow chamber was similar to the ones described by Hwang and Eaton 
( |2004l ); |Webster et aL] ( [2004] ); |Warnaars etaE] ( [2006] ); |Lu et al.| ( |2008| ); |Goepfert 



et al. (2010). The chamber by Hwang and Eaton was a 410 x 410 x 410 mm 3 



cubical Plexiglas box, from which the corners were cut off to make the inter- 
nal volume of the chamber closer to spherical. A loudspeaker-driven jet was 
mounted at each vertex, and aimed toward the center of the chamber. The 
loudspeakers were driven with sine waves, each with a random frequency and 
phase, to discourage the formation of standing waves or periodic structure 
inside the chamber; the power spectra were uniform in the range from 90 Hz 
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to 110 Hz. The approach by Hwang and Eaton had inspired the design of flow 
apparatuses by various workers interested in the generation of homogeneous 
and isotropic turbulence; see Webster et al. ( 2004[ ); Warnaars et al. (2006); Lu 



et al. (2008[);[Goepfert etaT (2010[). We compare the flow parameters in table 2.1 



2.1.3 The driving algorithm 

Here, we drove the loudspeakers with signals chosen to generate a flow with 
the desired anisotropy and low mean velocity. We achieved this by modulat- 
ing the amplitude of a sine wave with independent noise ( Fox et aTj 1988[ ) of 



correlation time 0.1 seconds, which was approximately equal to the large-scale 
eddy turn over time, Jz? /u', Jz? being a characteristic length scale describing the 



large-scale motions of the flow (see section 7.1 for a definition), and u' being 
the RMS velocity fluctuations. This condition ensured that fluctuations in the 
energy input rate to the turbulence occurred on time scales that were equal to 
or faster than the turbulence decay time, so that the turbulence was in a steady 
state. In any case, we found that the statistical properties of the turbulence were 
insensitive to the correlation time. 

At any given time t, the voltage for the i-th speaker was an amplitude- 
modulated sine signal 

a i = 5 i sm(2 7ift), (2.1) 

where 5i is the modulating noise and / is the base frequency of the sine wave. 
We drove the loudspeakers in phase and 50 Hz produced the strongest jet. 
The exponentially correlated noise 5i was calculated according to the algorithm 
given by Fox et al. ( |1988[ ). We first set the step size At, and the parameters D 
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and A. At each time step, we picked two random numbers p and q, uniformly 
distributed on [0, 1], and computed the following 

E = exp(— A At) , (2.2) 
h — [—2DX(1 — E 2 ) ln(p)] 1/2 cos(2 tt q) , (2.3) 
6i(t + At) = E5i(t) + h. (2.4) 



It follows that 5i is an exponentially correlated noise with the properties (Fox 

eTaTl|l988l ) 



<*(*)> = 0, (2.5) 



5i(s)) = D A exp(-A|t - s|) . (2.6) 
It can be seen froml2.6lthat A 1 is the correlation time for the noise. 



Figure 2.4[ a) shows a typical time series of the voltages generated at a sam- 



pling frequency 1/ At = 3000 Hz. A MATLAB® script to generate this is given 
in appendix [C] We had set A = 10 s _1 , D = 0.1 s, and the initial condition 
5i(t = 0) = 0. The sequence was about 40 s long and the RMS voltage was 1 V. 



Zooming into the sequence, figure 2.4 'b) shows the fine structure within 1 s of 
the signal. Because the correlation time of the noise is five times larger than the 
period of the carrier wave, the voltage amplitude varies slowly with time. 



Figure 2.5 'a) shows the unfiltered power spectrum of the voltage signal. The 



trend being masked by noise, in figure 2.5 'b) we filter the spectrum with a run- 



ning average to show the intended data. The spectrum was uniform up to about 
30 Hz, at which point the spectrum increased steadily with increasing frequency 
and reached a peak at the carrier frequency of 50 Hz, after which the spectrum 
fell off approximately inversely proportional to the square of the frequency. 
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10 20 30 40 

Time (s) 

(a) 




Figure 2.4: (a) A 40-second time series of the voltage applied to the loudspeak- 
ers, (b) A part of the time series, zoomed in to reveal the 50 Hz carrier wave. 
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Frequency (Hz) 

(a) 




Frequency (Hz) 

(b) 

Figure 2.5: The spectrum of the voltage applied to the loudspeakers (a) before 
filtering and (b) after filtering. The solid line in (b) is proportional to a power 
law with an exponent of —2. 
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Finally, to avoid mean flows, the amplitudes of each driving signal are con- 
strained so that they sum to zero. This is ensured by sending 

1 32 

ai ~ 32 a i ' ( Z7 ) 

to the i-th loudspeaker. If the response of the loudspeakers to the driving signal 
is linear, this adjustment maintains a constant volume of air in the flow appara- 
tus. It also reduces the amplitude of the sound generated by the loudspeakers, 
and minimizes the amount of air exchanged between the inside of the apparatus 
and the room. 

The code for driving the loudspeakers was programmed in a Lab View envi- 
ronment. The 32 independent voltage signals were generated with a National 
Instruments 32-channel DAQ card (NI PCT6723) on a computer and amplified 
by an amplifier unit before being sent out to the loudspeakers. 



2.1.4 The distribution of amplitudes 



The distribution of loudspeaker amplitudes is shown in figure 2.6 By calibrat- 
ing the RMS amplitude of each loudspeaker along the surface of the soccer ball, 
we were able to select the desired anisotropy. As we restricted ourselves to 
cylindrically symmetric forcing, the asymmetry can be characterized by the ra- 
tio of the axial amplitude, & ax i a i, to the radial amplitude, 6 ra diai : 

A = ^ . (2.8) 

"radial 

Here, the RMS amplitudes are taken as 6 axial = (a axia i) 1/2 and Radial = K 2 ad i a i) 1/2 / 
where (a axial ) and (ar adial ) refer to the RMS amplitude of the polar and equato- 
rial loudspeakers, respectively. A forcing is then described as exhibiting oblate 
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b axial 




Figure 2.6: A sketch showing the assignment of RMS amplitude for a given 
loudspeaker. The vector pointing from the center of a sphere, O, to the center of 
a loudspeaker, S if intersects an inscribed ellipsoid at Ei. The distance of E; L from 
the center of the sphere sets the RMS amplitude 6j for the given loudspeaker. 



(pancake) asymmetry when < A < 1, spherical symmetry when A — 1, and 
prolate (cigar) asymmetry when A > 1. 

As we have 32 loudspeakers distributed in an icosahedral symmetry, the 
RMS amplitude for a given loudspeaker was set according to the following 



scheme. As shown in figure 2.6 a vector was drawn from the center of the 
sphere to the center of the loudspeaker. We calculated the intersection between 
the vector and the surface of an inscribed ellipsoid, which had the same center 
as the sphere. The distance between this intersection and the center of the soccer 
ball set the relative RMS amplitude for the given loudspeaker. Given an asym- 
metry A and any one of the principal radii, 6 axia i and & ra iai/ the geometrical shape 
of the ellipsoid is uniquely specified. Therefore, the distances may be precalcu- 
lated once and used for setting the RMS amplitudes for all the loudspeakers. 



Figure 2.7 'a) shows a typical fluid velocity time series when the 32 loud- 
speakers were driven according to the prescription described above. Fig- 
ure |2.7| b) shows the velocity spectrum calculated using Lomb algorithm for 
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unequally spaced data[j](see e.g. Press et al. (2007)). We did not detect the 50 Hz 
oscillation in the velocity spectrum. We think that any traces of such periodic 
excitation may have been erased as a result of the fluid mechanical interactions 
between multiple turbulent jets. 



2.1.5 The amplifier 

Sixteen Raveland XCA 700 two-channel amplifiers, each being capable of de- 
livering up to 700 Watts of power, amplified the digital signals sent through a 
PCI-to-BNC interface. The design and building of the amplifier unit was done 
by Ortwin Kurre and his technical assistants from the electronic shop at the Max 
Planck Institute for Dynamics and Self-Organization in Goettingen, Germany. 



2.2 Laser Doppler Velocimetry 

Much of the information about the flows we generated inside the soccer ball was 
collected using laser Doppler velocimetry (LDV). It is a technique for extracting 
the components of the velocity of individual particles from the light scattered 



by them (e.g. Albrecht et al.| 2003 ). A laser light source is an essential part of the 



technique. As shown in figure 2.10[ two beams of collimated, monochromatic, 



and coherent laser light are focused by a converging lens and made to cross at 
their respective focal point, where they interfere and generate a set of straight 
fringes. Focusing is not essential to the explanation that will be made, but is 
only required for making the intersection volume small so that only one par- 



1 The MATLAB® code by C. Saragiotis for calculating spectrum with Lomb algorithm is 
greatly acknowledged. 
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Figure 2.7: The figure shows (a) a truncated time series of one component of the 
fluid velocities in response to the forcing and (b) its spectrum calculated using 
the Lomb algorithm for unequally spaced data. 
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Figure 2.8: The Raveland amplifier. 

ticle can be inside this volume at any given time. As individual particles pass 
through the fringes in the measurement volume, defined by the crossing of the 
two beams and is approximately ellipsoidal in shape, they reflect light from the 
regions of constructive interference into a photo detector. The component of ve- 
locity perpendicular to the fringes can be calculated from the frequency of the 
signal received at the detector. There exist, in the literature, two explanations 
for the frequency dependency of the particle velocity (e.g. Albrecht et al.| 2003[ ). 



Both yield the same result but the explanations differ slightly in their starting 
point and emphasis. We review them in the following sections and discuss the 
limitations of LDV for turbulence measurements in air. 



52 



* II 




n 


i 


■ 


■ 


■ 
■ 


■ 


■ 


■ : 


■ \ 


i 


Si 




§ 


HI ■ 

9i ■ 




Amplifier la - Id 
Amplifier 2a - 2d 

Amplifier 3a - 3d 
Amplifier 4a - 4d 

Amplifier 5a - 5d 

Amplifier 6a - 6d 

Amplifier 7a - 7d 

Amplifier 8a - 8d 

BNC connectors 
32 Analog 
4 Digital 

Fans 



Figure 2.9: The amplifier unit. 
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Figure 2.10: A schematic of a laser Doppler velocimeter. 




Laser 



Receiver 

Figure 2.11: The figure illustrates the change in received frequency caused by 
the motion of the particle in a single-beam scattering process. 



2.2.1 The Doppler model 



The first explanation is based on the Doppler effect. Consider the scenario in 



figure 2.11 in which a stationary monochromatic light source with wavelength 
X s and frequency u s impinges on a particle moving with velocity v p in the flow. 
Denoting the Doppler frequency perceived by the particle as v v , we have 



v P = v, I 1 



e s - v p 



(2.9) 



where e s is the unit wave vector of the incident beam and c is the speed of 
light. Invoking the Doppler effect a second time, the frequency received by a 
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stationary detector is 

1 - e s -v p /c 



(2.10) 



1 — e r -v p /c \1— e r -v p /c / 

where e r is the unit wave vector of the scattered beam at the receiver's end. For 
typical turbulent flows, the particle speed is much less than the speed of light. 



Thus, we have for the expression in equation 2.10 the Taylor expansion 



v s + Vp ' [e : Gs) (hg« C)C =A s ^). (2.ii) 

As 



For typical flow velocities (1~100 m/s), the Doppler shift frequency is of 
the order 1 MHz to 100 MHz, and is contained in the light frequency, which is 
approximately 10 14 Hz. To resolve this small fraction of Doppler shift from the 



light source is technically challenging, if not at all impossible; see Paul and Jack 



son| ( |1971[ ). In order to make the Doppler frequency lie in a more manageable 



range, two coherent laser beams of equal intensity and wavelength are used 



instead of one. Figure 2.12 shows the dual-beam scattering process. When a 
particle moves across the measurement volume formed at the intersection of 
the two beams, it scatters light from both beams, each with a frequency 

v p - (e r - ei) Vp - (e r - e 2 ) . . 

v\ = v s -\ , i/ 2 = u s H - . (2.12) 
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The light is focused and mixed on the detector to yield the beat frequency 

v p -(e 1 -e 2 ) 

vb = v 2 - vi = . (2.13) 

As 

If the angle formed by the intersection of the two beams is denoted by 9, and the 
component of particle velocity perpendicular to the bisector of the beams (the 
z-axis) is denoted by v px , then the beat frequency can be expressed as 

vb = M . (2.14) 

Thus, the frequency difference is linearly proportional to the velocity compo- 
nent perpendicular to the bisector of the two beams. 



2.2.2 The fringe model 



The expression in equation 2.14 is derived largely based on kinematic consider- 
ation and it does not explain the origin of the beat frequency. The fringe model 
fills this gap of our understanding. We illustrate it using the configuration in 



figure 2.12 If two incident beams of equal amplitude E and frequency tu s , and 
linearly polarized perpendicular to the plane in which the two beams lie, then 
the electric field of each beam at an arbitrary point r in space is described by 

Ei = e y E exp(i[u s t - k x -r + 4>i\) , (2.15) 
E 2 = e y Eq exp(i[u s t - k 2 • r + fa]) . (2.16) 

ki and k 2 are the wave vectors of the two beams, given by 

k 1 = A; s ei , k 2 = k s e 2 , k s = 2ir/\ s . (2.17) 

(pi and 2 are the initial phases of the waves at the origin, respectively. At the 
crossing of the two beams, the two electric fields add according to the superpo- 
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sition principle to give an intensity / proportional to 



I oc |Ei + E 2 | 2 . 



(2.18) 



By substituting the expressions for the electric fields from equations 2.15 2.16 



and 2.17 we have for the intensity 



/ oc 2 E 2 



l + cos(27r( 2sin A (g/2) ) x - {fa - fa) 



(2.19) 



It can be seen that the electromagnetic wave intensity varies periodically in x, 
with a fringe spacing given by 



Ax 



(2.20) 



2 sin(0/2) ' 

A particle of diameter d p <C Ax passing through the intersection volume sam- 
ples the local intensity of the interference pattern and produces a burst of flick- 
ering light. If At is the particle's fringe crossing time, then the frequency of 
flickering is given by 

1 2 sin(0/2) Ax 2 sin(0/2) 



^detector 



At 



A, 



At 



A, 



' px I ) 



(2.21) 



which is precisely equation 2.14 The fringe model offers a very physical and 
intuitive approach for understanding the laser Doppler technique. This model 
is only valid when the particle diameter is less than the wavelength of light. For 
particles larger than the wavelength of light, the arguments have to be modified 
in order to take into account the additional phase shift due to wave propagation 
inside the particles. This additional phase shift only causes the fringe pattern 



to be shifted in phase and equation 2.21 remains valid. The reader is referred 
to chapter 2 of Albrecht et al. ( 2003[ ) for a full account of the particle finite size 
effect. 
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Figure 2.13: A schematic of a Bragg cell. A transducer generates an acoustic 
wave, resulting in a moving periodic refractive index grating that diffracts the 
incident light similar to in Bragg diffraction. The arrow points in the direction 
of travel of the acoustic wave. 



2.2.3 The Bragg cell 

The reader may note that particles with equal, but opposite velocities, moving 
across the measurement volume formed in the dual beam arrangement will pro- 
duce the same frequency of flickering. To resolve the particle velocity direction, 
a Bragg cell is placed in one of the beam paths to modify the Doppler frequency. 
It consists of a mechanical oscillator, such as a piezoelectric transducer, attached 
to an optically transparent material, such as quartz. An oscillating electric field 
drives the transducer to vibrate, emanating a traveling acoustic wave. Under 
the mechanical pressure of the wave, a periodic refractive index grating results. 
This traveling periodic grating diffracts the incident light much like the atomic 
planes of a crystal in Bragg diffraction. 



If a Bragg cell is placed in the path of beam 1 in figure 2.12 the frequency of 
the incident beam is shifted by an amount to v\ + z/ s h when the acoustic wave 
moves toward the beam, and V\ — z/ sh when the acoustic wave moves away from 
the beam. In the fringe model, the first case corresponds to a movement of the 
fringes in the —x direction, whereas the second case corresponds to a movement 
of the fringes in the +x direction. The light from beam 1 and 2 is mixed on the 
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detector to yield a beat frequency 

2 sin(0/2) c . . , . 

^sh H t e x . • v p fringes move in the — x direction , 

UB= < 9 ■ \,n ' (Z22) 

2 sm(6 , /2) . . , ,. . 

z/ sh • v p fringes move m the + x direction . 

X s 

It is easily seen that a stationary particle will yield a signal with frequency u Bh , 
a movement in the direction of the fringes a lower frequency, and a movement 
against the fringes a higher frequency 



Figure 2.14 shows an overview of the TSI LDV system used in our inves- 
tigation. An argon ion laser by Coherent (Innova 90C-5) provided a multi- 
colored beam. A multicolor beam separator (TSI Model 450200), converted 
the beam into three beam pairs of different wavelengths: green, blue, and vi- 
olet. A Bragg cell in the beam separator shifted one in each beam pair by 
40 MHz. Two fiberoptic probes focused the light from the beams to form a 
measurement volume. One probe measured two components, and the other 
one only one component, of the velocity of individual particles that passed 
through the approximately ellipsoidal measurement volume, which measured 
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100 \xm in diameter and 2 mm in length. The effect of a finite measurement vol- 
ume on the measurement of particle velocity is discussed in greater detail in sec- 



tion 2.2.5 This resolution, however, was sufficient to resolve scales larger than 
the dissipation scale of the turbulence in this experiment. A receiving fiber in the 
probe collected light scattered by particles in the volume and sent the light to a 
photodetector (PDM 1000). The photodetector converted the scattered light into 
electrical signal which was then sent to a signal processor (FSA 4000). The sig- 
nal processor extracted Doppler frequency information from the electrical sig- 
nals and sent the digitized signal to be analyzed by software (FLOWSlZER™) on 
a computer for data acquisition, yielding a single velocity measurement. Data 
consisted of a stream of velocity samples with irregular time intervals between 
them. The mean data rate was between 300 and 3000 samples per second, de- 
pending on the probe. For statistical convergence, we typically collected 2 x 10 6 
data points per aspect ratio of the forcing and per spatial position of the probes. 



2.2.4 Data processing 

In LDV measurement, the flow velocity is sampled in time by random passages 
of particles that go through the measurement volume. This has two conse- 
quences. First, the velocity samples are irregularly spaced in time. Second, 
the measurement volume is sampled by high-speed particles more often than 
low-speed ones. As a consequence, even if the particle spatial distribution is 
independent of the velocity field, the rate of particle arrival at the measurement 
volume is correlated with the velocity field ( |McLaughlin and Tiederman 1973[ >. 



Without careful attention to this latter fact, statistical average of flow quantities 
can be biased. Thus, velocity samples must be weighted by appropriate statis- 
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tical factors to yield unbiased statistical averages (McLaughlin and Tiederman 
1973[ ). The statistical factor we used for correcting single-point statistics is the ar- 
rival time weighting factor (see e.g. Adrian and Yao 1987[ ). Other methods may 
also work (see e.g. |Buchhave| \1975\ |Buchhave et al.f |1979f |McDougaIIj 119801 ), but 
the arrival time method has the advantage of being computationally efficient 
and is well-suited for non-homogeneously seeded flow fields. 

The method works by weighting the observables with the arrival times be- 
tween particles. For example, the estimator for the mean and the RMS fluctua- 
tions are 



EjliVi (ti-*i-i) 



(v) 



i-l 



T^ =x {vi-{v)) 2 {ti-t 



i-l) 



T N t -t 



(2.23) 
(2.24) 



i-l 



where N is the total number of velocity samples. As noted by Albrecht et al. 



(2003), this method is suitable for moment estimation of single point statistics 



but fails for estimation of correlation functions. Therefore, we did not apply 
this method to our two-point correlation measurements, and we did not take 
any measures to correct for possible bias in our two-point correlations. 

Because different probes record slightly different particle arrival times, in or- 
der to calculate two-point statistics, one must look for coincidences of velocity 
samples coming from these probes. We found these coincidences by discretiz- 
ing the velocity time series using a slotting technique ( Mayo et al.| 1974) and 
searched for samples coming from each of the two probes with time separations 
falling within a given temporal bin, At, of duration 0.5 ms. This duration was 
chosen because it maximized both the number of coincidences and the value of 



the velocity autocorrelations (see section 7.1 ) near zero separation while keeping 



the value of the structure functions near zero separation small (see figures 2.15 
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Figure 2.15: The effect of temporal binning on the value of the velocity correla- 
tion, g(r), near zero separation. The values in the legend are the separations in 
millimeters. 
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Figure 2.16: The structure function as a function of temporal bin width. For 
legend, see previous figure. 



2.16 and 2.17). We typically collected about 10 4 — 10 5 velocity samples. 



For each of the two probes, an average was taken of velocity signal coming 
from the same probe falling within At. The product of all averaged velocity 
pairs was then added to a sum. After processing the entire temporal sequence, 
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At (s) 

Figure 2.17: The number of samples yielding the velocity correlation as a func- 



tion of temporal bin width. For legend, see figure 2.15 



the total sum was divided by the number of accumulated products, giving us 
the value of the velocity autocorrelation for that given spatial separation. We 
found that so long as At was smaller than 15 ms, the difference in the values of 



the correlation was less than 1%, as shown in figure 2.15 



By definition, the velocity correlation is unity at zero separation. Noise, how- 
ever, may contaminate the data and mask the true value of the correlation at 
zero separation. In order for the correlation to comply with this definition, we 
normalized the velocity correlation with its value at zero separation. The mea- 
sured structure functions could have been corrected to yield better accuracy at 
the smaller separations, but since the measured structure functions had leveled 
off at the low end (see figures 5.1 and |5.3[ >, and correction of any data is never a 
satisfying proposition, no matter how carefully executed, because in some sense 
it presumes the answer, therefore no attempt was made to apply corrections to 
these measurements, but they may be required in some future investigation. It 
can certainly be said that the uncorrected data were internally consistent and the 
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inertial scale measurements were done with sufficient accuracy that no special 
correction needed to be developed. 



2.2.5 Limitations of the LDV technique 



In addition to the bias discussed in 2.2.4 other physical parameters may limit 



the accuracy of LDV. We review these limitations in the following by dividing 
them into those arising from the flow and those arising from measurement in- 
struments. 



Tracer particles 

The motion of an individual spherical particle in a turbulent flow has been ex- 
tensively studied (see e.g. Corrsin and Lumley| 1956} Manton 1977| Maxey and 
Riliyj pgfy and the response of a tracer particle to isotropic turbulence has 
been examined by |Mei ( 1996 ). In making accurate velocity measurements of a 
fluid with LDV, we require the velocity of the tracer particle be close to the lo- 
cal, instantaneous fluid velocity. The particle inertial effects are described by the 
Stokes number 

(2.25) 



7V, 



where t p is the particle response time defined as 

1 fg P \ ( D 2 



TV, 



(2.26) 



where g p , g /, D, and v are the particle density, the density of the fluid, the par- 
ticle diameter, and the fluid viscosity, respectively. t v is the Kolmogorov time 



defined in equation 1.10 Buchhave et al. (1979) considered the particle-fluid 
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relative velocity, using a model by Lumley ( 19761 ), and estimated an upper limit 
of the particle Stokes number of St ^ 1/74 with the criterion that the error in 
the velocity be less than 1%. We used oil droplets (density 0.9 g cm 3 ) as tracer 
particles in air (density 1.2 kg m 3 ), generated by a Palas AGF 10.0 aerosol gen- 
erator. The droplets had a most probable diameter of about 3 \xm. They settled 
in air at 200 \xm s _1 . The Stokes number based on the most probably diameter 
was about 0.02, the smallness of which qualifies these oil droplets as passive 



tracers (see e.g. Bewley et al. 2008). 



A second problem connected to seeding the fluid with particles is turbulence 
modification due to particle-fluid interaction. This is an area of active research 
(see e.g. He stroni| 1989| Elghobashi and Truesdell 1993). The addition of parti- 
cles at a volume fraction of as low as 10~ 5 to a turbulent flow modifies its mo- 
tion (Elghobashi 1994[ ). For low values of particle volume fraction (^ 10 6 ) the 
particles have negligible effect on the flow (Elg hobashi} 1994[ ). The oil droplets 
used in our experiment were present at a volume fraction of less than 10~ 7 . This 
particle concentration was too small to alter the turbulence. 



The effect of refractive index fluctuations 

Since the laser beams of the LDV probes must propagate through a medium 
having random index of refraction fluctuations, the beam paths will no longer 
be straight, but will take on random fluctuating passages through the medium, 
and the phase and amplitude of the receiving signal will fluctuate randomly. 
The fluctuations in the index of refraction are primarily due to pressure fluc- 
tuations and temperature fluctuations in the medium. Buchhave et al. ( 1979| ) 



estimated that, by considering an isentropic variations in the pressure in a com- 
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pressible flow and assuming a quasi-Gaussian relationship between pressure 
and velocity fluctuations, a root-mean-square (RMS) velocity fluctuations of the 
order of 70 m s 1 are necessary to produce significant changes in the index of re- 
fraction. For temperature fluctuations, they obtained an estimate for the rate of 
decrease of index of refraction with temperature of the order of 1 x 10 6 K 1 for 
air at 273 K. Since the typical velocity fluctuations in our experiment was lms" 1 
(<C 70 m s -1 ), and the rate of change of index of refraction with temperature is 
small, we can neglect the effect of index-of -refraction fluctuations. 



The probe volume size 



The probe volume is defined by the surface on which the light intensity of the 
fringes is 1/e 2 of the maximum intensity in the beam intersection region. With 
Gaussian beams the probe volume is an ellipsoid, and we have the expression 
for the probe volume dimensions in terms of the principal diameters of the el- 



lipsoid in the coordinate system defined in figure 2.18 



4 = — , (2.27) 

COS K 

d y = d , (2.28) 

4 = — , (2.29) 
sin k 

where d , the width of the waist of the focused beam, is 

4A f 

rfo = — , (2.30) 

7r a 

for A, the wavelength of the light, /, the focal length of the transmitting lens, and 
a, the diameter of the beam. For our LDV system, the probe volume formed by 
the two beams of wavelength A = 514.5 nm had the largest dimension. Each 
beam had a diameter of 3.5 mm, and for our optics, / = 512 mm, n = 2.8°, 
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the probe volume dimensions were d x « d y ~ 96 |J.m, and = 1.97 mm. In 
order to resolve the smallest length scales in the flow, we required that any of 
the three probe volume linear dimension be smaller than the Kolmogorov scale, 



rj (see equation 1.9). Our probe volume dimensions were d x rs 0.6r),d y m 0.6 rj, 



and d z ~ 12 77 (see table for values of 77 in the experiment). Thus, our optical 
setup had not met this condition but it had sufficient dynamic range to resolve 
the intermediate and the very large scale motion. 
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CHAPTER 3 
THE EXPERIMENTAL PROCEDURE 



Here we describe the coordinate system we used to present our results. We 
describe the measurement protocol that exploited the symmetry of the flow gen- 
eration system to rotate the turbulence relative to our LDV measurement appa- 
ratus. We discuss the quantity we fixed when comparing flows with different 
values of the forcing anisotropy. 



3.1 The coordinate system 



The turbulence coordinate system was not fixed to the laboratory frame, 



(xi, X2, Xs), as is shown in figure 3.1 but was fixed to the turbulence symme- 
try axis. Because the turbulence was axisymmetric, it had two principal axes; 
one being the axis of symmetry of the forcing, and the other one, a radial axis, 



perpendicular to it. However, as described in 3.2 our measurement apparatus 
measured particle velocities only at points that lie along a straight line fixed in 
the laboratory frame. In order to sample the statistics of turbulence along both 
axes, we rotated the axis of symmetry of the turbulence by taking advantage of 
the symmetries of the flow apparatus. 

The two orientations of the symmetry axis in the laboratory frame are shown 



as x\ and x' 3 in figure 3.1 each of which passes through the centers of opposite 
hexagonal faces of the flow apparatus. The x\ axis was nearly parallel to the 
line along which we collected velocity samples, x\, and it lay along 9 = 107° 
and <p — — 4°. The x' 3 axis lay along 9 = 4° and <p — 180°, and was close to 
the £3 axis. These two axes, x[ and x' 3 , were also nearly perpendicular and 
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Figure 3.1: The schematic shows the coordinate system in the laboratory frame. 
X3 was aligned with the vertical. The axis of symmetry of the forcing lay along 
either x[ or x' 3 , depending on the configuration. Our measurements were made 
at points lay along x±. 



parallel to the vertical axis, respectively. We assume that the two orientations 
are equivalent because the acceleration of gravity is negligible relative to particle 



accelerations in the turbulent field (see e.g. Voth et al. 2002 ). Because the primed 



and unprimed coordinate systems are close to each other, we do not distinguish 
between them in the rest of this thesis. 



Figures 3.2 'a) and (b) show the coordinate system, (r 1 ,r 2 ,z), which was 
aligned with the symmetry axis of the forcing, z. The two orientations of the 
forcing coordinate system with respect to the laboratory frame correspond to 
the two cases described above. Hereafter, 'axial' refers to both the direction of 
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(a) (b) 

Figure 3.2: The figure shows the orientations of the body coordinate system (in 
red) of the forcing (ellipsoidal surface in yellow) with respect to the laboratory 
frame (gray), and the conventions we used: (a) when the axis of symmetry of 
the forcing lay close to the x 3 axis of the laboratory frame, and (b) when the axis 
of symmetry lay close to x x axis of the laboratory frame. The coordinate system 
(xi, x 2 , x 3 ) was fixed in the laboratory frame, and the coordinate system (r 1: r 2 , z) 
was fixed with respect to the symmetry of the forcing. 

the particle velocity component measured along the axis of symmetry, and, in 
discussing two-point statistics, separations along the axis of symmetry. Like- 
wise, 'radial' refers to both the direction of the particle velocity components 
normal to the axis of symmetry, and to separations normal to the axis of sym- 
metry. Additionally, we refer to two-point quantities whose separation vector 
lies along the axis of symmetry as axial quantities, and those with radial sepa- 
rations as radial quantities. For example, a two-point quantity f(xi) is denoted 
as f(z) when the axis of symmetry of the forcing lies along x\, and is called 
an axial quantity. In keeping with the geometry described above, our 'axial' 
measurements were in reality about 15° away from the axis of symmetry, and 
the 'radial' measurements were about 4° from its normal. Our results do not 
depend on the angle between the 'axial' and 'radial' measurements. 
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z 




(a) (b) 



Figure 3.3: The figure shows the orientations of the measured velocity compo- 
nents with respect to the forcing coordinate system when the separation vector 
lay (a) along the axis of symmetry, or (b) perpendicular to the axis of symmetry 



3.2 The measurement protocol 



We used the LDV probes in two configurations. In one configuration, all three 
orthogonal components of the fluid velocities were measured at a single point. 
For these measurements, the probes observed a fixed volume (aligned to within 
10 |a.m with a pin hole) close to the center of the soccer ball, which is the point 



(0, 0, 0) in the laboratory frame shown in figure 3.1 One probe measured a 
single component of the fluid velocities, namely u xi (0, 0, 0). The second probe 
measured two orthogonal components of the fluid velocities, namely u X2 (0, 0, 0) 
and ^(O, 0, 0). In the second configuration, we aligned the probes in a simi- 
lar way, except that we positioned the two-component probe at different points 
along the x\ axis, using a programmable linear traverse. This probe now mea- 
sured u X2 (xi, 0, 0) and u X3 {x\, 0, 0). In addition, the single-component probe was 
rotated by 90° to measure u X3 (0, 0, 0), which was coincident with one of the com- 
ponents measured by the two-component probe when x\ equaled zero. These 
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velocity components, expressed in the {r l) r 2l z) coordinate system, are shown 
schematically in figure 3.3 The two orientations of the forcing coordinate sys- 
tem relative to the laboratory frame are described in figure 3.1 



3.3 Methods 



We collected two data sets for each of a series of forcing anisotropy, A, given by 



equation 2.8 For one data set, the axis of symmetry of the forcing was aligned 



close to £3, and for the other, the symmetry axis was aligned close to x\. For 
each value of the forcing anisotropy, and each orientation of the forcing sym- 
metry axis, data were collected with the two-component LDV probe positioned 
at various locations along the xi-axis. We varied the forcing anisotropy while 
fixing the turbulent kinetic energy in the center of the soccer ball 



K 



< (0,0,0)) + « (0,0,0)) + « (0,0,0)) 



(3.1) 



where u' x . is the RMS fluctuations of the fluid velocity in the %i direction, and 
(• • ■ ) denotes temporal averaging. 
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CHAPTER 4 
THE FLOWS 



In this chapter, we evaluate the quality of the flows generated by the tur- 
bulence apparatus introduced in chapter [2] Three-dimensional velocity mea- 
surements made at single points in the center region of the flow apparatus 
show that the anisotropy of the turbulent velocity RMS fluctuations followed 
the anisotropy of the forcing signal. Two-point velocity measurements made at 
points that lay along a single line in the laboratory frame show that the flows 
were approximately homogeneous and axisymmetric. For the case of spheri- 



cally symmetric forcing, with A = 1 (see equation 2.8 1, we expected and indeed 
found that the turbulence was isotropic. In all cases, the mean velocities were 
less than 0.24 ms _1 , or less than 23% of the fluctuations and the shear stresses 



were less than 7% of the turbulent kinetic energy, as shown in table 4.2 
4.1 Anisotropy of the fluctuations 



Figure 4.1 shows the ratio of axial to radial velocity fluctuations as a function 
of the forcing anisotropy, A. It is evident that the anisotropy of the turbulence 
followed the variation of A. In addition, the anisotropy was nearly the same 
whether measured in the r% or r 2 directions. Therefore, the turbulence was close 
to cylindrical symmetric, though this was less so at extreme values of A. The 
loss of cylindrical symmetry might be explained by the decrease in the number 
of loudspeakers driving the turbulence as the value of A moved away from 



one. This is the nature of the forcing algorithm described in section 2.1.3 We 
think that, as the number of loudspeakers effectively decreased, the turbulence 
became more sensitive to mechanical differences between the loudspeakers, and 

74 



2.5 

g 2 

4 1 



0.5 





A 

A O 

6 

prolate 


■I 

oblate ^ 

I 


spherical 



0.5 



1.5 

A 



2.5 



Figure 4.1: The velocity fluctuation ratios as a function of the forcing anisotropy, 
A. The circular symbols (O) are ■>/ (uf (0) ) / (u® (0) ) and the triangular symbols 
(A) are ^/ (uf (0)) / (ug(0)) . Symbols in the region shaded in blue have oblate 
shape asymetry, whereas those in the region shaded in red have prolate shape 
asymmetry Spherical symmetry lies at the intersection of the two regions. Data 
are collected with the axis of symmetry lying along x' 3 . We obtained similar 
results (not shown) when the axis of symmetry lay along x\ . 



to misalignments of the nozzles. 



4.2 Axisymmetry of the fluctuations 



Here, we characterize the degree to which the axisymmetry was spatially uni- 



form. In figure 4.2 'a), we examine the ratio of the two radial fluctuating veloci- 
ties, a/ (w' r 2 2 (0, 0, z))l 1 Jvfl (0, 0, z)), as a function of the distance from the center of 
the soccer ball along the axial direction for three values of the forcing anisotropy. 
Within 50 mm from the middle of the flow chamber, the values of this ratio de- 
viated by less than 10% from the value 1, indicating that the turbulence was 
close to being cylindrically symmetric. To gauge the extent of spatial unifor- 



mity of the anisotropy, we plot in figure 4.2 b) the ratio of the axial fluctuating 
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Figure 4.2: For three values of the anisotropy measured at the centre of the ball, 
a/ (w /2 (0)) I (u'?(0)), we plot (a) the ratio of the two radial components of the ve- 
locity fluctuations (w' r 2 2 (0, 0, z))/(m^(0, 0, z)) and (b) the ratio of the axial fluc- 
tuations to the radial fluctuations, J Jvffir\, 0, 0))/(ii' r 2 2 (ri, 0, 0)) at various points 
moving away from the center of the ball. The values in these legends, and in all 
others in this thesis, are those of the data in (b) for r\ = 0. 
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velocity to one of the radial fluctuating velocities, a/ (uf(ri, 0, 0))/ (u'? 2 (ri, 0, 0)), 
as a function of the distance from the center of the chamber. Here, the measure- 
ment points lie on a line nearly normal to the axis of symmetry, and the values 



of this ratio at r± = correspond to the values shown in figure 4.1 Again, the 
velocity fluctuation ratio is approximately constant within 50 mm of the center 
of the chamber. Thus, we inferred that the degree of anisotropy was approxi- 
mately uniform within a central spherical region with radius 50 mm. In both 
plots, the error bars were the standard errors of the measurements obtained by 
truncating the original data sets into 20 equal parts. They indicate the sampling 
accuracy and not the measurement accuracy. These error bars were always less 
than ±8.5%. We concluded that, within the error bars, the ratios of the fluc- 
tuating velocities at different locations from the center of the soccer ball were 
indistinguishable from each other. 

Since we take axial and radial measurements through the turbulence re- 
gion, we obtained twice the information on the radial velocity. Thus, the ratio 
a/ (m' z 2 (0))/(m^ 2 (0)) in our experiment may be calculated in two different ways. 
The first one is the ratio of the two orthogonal velocity fluctuations in the limit 

as n — > 

</< = lira J(uf(r u 0,0))/(ug( ri ,0,0)) , (4.1) 

7*1 — HI V 

and the second one is 

</< = lim ^(uf( ri ,0,0))/(ug(0,0,z)) . (4.2) 

z-»0 
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We calculated the ratio in each case by averaging the RMS velocity fluctua- 
tions at locations within 10 mm from the center of the chamber, and then taking 
the ratio between the axial and radial fluctuations. Their values are summa- 



rized in table 4.1 The variation in the fluctuations in this range was less than 
5% and we used this variation for estimating the error in the velocity fluctua- 
tions. Our single-point statistics do not depend on this choice of the definition 
for the anisotropy ratio, and we discussed the results in terms of u' z /u' r2 defined 
4.1[ The two-point correlation functions, as we will show in chapter [7| do 



m 



depend on the anisotropy of the fluctuations and we discuss the results in terms 



of u' z /u' r2 defined in 4.2 since these are the velocity fluctuations associated with 
the radial and axial correlation functions. 



4.3 Homogeneity of the fluctuations 

To evaluate the homogeneity of the turbulent fluctuations, we present mea- 
surements of y/(u%(0,0,z))/(ug(0,0,0)) and y/(uf{r 1} 0, 0))/« 2 (0, 0, 0)) in fig- 
ures 4.3 These ratios compare the amplitude of the velocity fluctuations at loca- 
tions away from the center to those in the center of the chamber. The error bars 
were the standard errors of the measurements obtained by dividing the original 
data sets into 20 equal parts. Within a radius of 50 mm, the difference between 
y/(ug(0, 0, z))/(ug(0, 0, 0)) and its value at the origin is within 5% of the value 
at the origin. The inequality also holds for a/ {vfj(rx, 0, 0))/(tt' 2 2 (0, 0, 0)). Thus, 
we conclude that the RMS fluctuations are homogeneous. 
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Figure 4.3: The comparison of fluctuations at various distances from the center 
of the soccer ball to those in the middle, (a) Radial fluctuations along axial sep- 
arations, v/«(0, 0, z))/« 2 2 (0, 0, 0)), for anisotropy ratios, \/(uf(0))/(u?(0)), 
ranging from 0.59 to 1.71. (b) Axial fluctuations along radial separations, 
a/ (uf(ri, 0, 0))/(m' z 2 (0, 0, 0)), for the same range of anisotropy. The anisotropy 
ratios in the legend are measured at the center of the soccer ball. 
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4.4 Anisotropy of the mean flow 



Because anisotropy can manifest itself not only in the fluctuations, but also 



through spatial variation of the mean velocity, we show in figure 4.4 'a) 
and (b) the ratio of the mean flow velocity to the velocity fluctuations, 
(U n {W,z))l(u%{W,z)) x l 2 and (U z {n, 0, 0))/« 2 (r 1; 0, 0)) 1 / 2 . It is difficult to 
predict a trend in the data, but may indicate that the actual center of the tur- 
bulence was displaced from the center of our coordinate system. We concluded 
that the variation in the mean velocity was negligible because it was typically 
less than 10% of the fluctuations. Only at the extreme values of the anisotropy 
was the mean velocity as much as 15% of the fluctuations. This may have been 
due to the sensitivity of the turbulence to small differences between the loud- 



speakers in this range of anisotropics, as discussed in section 4.1 



4.5 Reflectional symmetry 



As described in section 1.6 reflectional symmetry in a plane containing the axis 



of symmetry implies that (e.g. Lindborg 1995 ) 

«(0,0,*)< (0,0,0)) 



Xr 



(4.3) 



v /«(0,0,z))(«£(0,0,0)) 
In addition, reflectional symmetry in a plane normal to the the axis of symmetry 
implies that 



X 



zr 2 



Kfo, o,ok 2 (o,o,o)) 

V<««(r 1 ,0,0))(««(0,0,0)) 



(4.4) 



The measurements of these correlations are shown in figure 4.5 The shape of 



the curves are difficult to interpret, but may suggest that the symmetry planes 
of the turbulence were slightly tilted from the plane of measurement. In both 
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Figure 4.4: The mean flow as a fraction of the fluctuating velocity (a) 

(U n (p,0,z))/y/(u%(0,0,z)), and (b) (r/,(r a ,0,0))/V«(ri,0,0)) for axial-to- 
radial velocity fluctuation ratios ranging from 0.59 to 1.71. The values in the 
legend are the anisotropy measured at the centre of the soccer ball. 
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Figure 4.5: The cross correlations of orthogonal velocity components at various 
distances from the center of the soccer ball (a) along the axis of symmetry, and 
(b) along an axis perpendicular to the axis of symmetry of the forcing. The 
values in the color bar are the anisotropy measured at the center of the soccer 
ball for various values of the anisotropy of the forcing. 
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plots, the error bars were the standard errors of the measurements obtained by 
truncating the original data sets into 20 equal parts. These errors were always 
less than ±3%. We do not believe that the trend and the scatter in the data were 
significant enough to be the basis of a theory Because the variation in the cross 
correlations was always less than 10% (the mean was no more than 5%), we 
concluded that the turbulence possessed remarkable reflectional symmetries 



4.6 Reynolds shear stresses 



Two observations of turbulent shear flows may be useful. First, results obtained 



in the strongly sheared part of the boundary layer by Saddoughi and Veeravalli 



(1994) indicated that the width of the Kolmogorov scaling range of the second 
order transverse structure function is significantly shorter than that of the longi- 



tudinal function (see figure 1.2). In an overview of turbulence spectra obtained 



in various flows, Sreenivasan ( 1995| ) noted that shear-flow turbulence exhibits 
inertial range scaling in the transverse spectrum only at values of R\ larger 
than 10 3 . A later study by |Noullez et al.| ( |1997| > lowered this threshold to 500. 
Second, there is evidence from experiments (e.g. |Stolovitzky and Sreenivasan 



1993 Shen and Warhaft 2002[ ) and numerical simulations (e.g. Benzi et aTj 1996[ > 
that the presence of a strong shear at low Reynolds numbers may destroy the 
extended self-similarity method (ESS) introduced by Benzi et aT| (1993) (see sec- 



tion 1.4). We surmise from these observations that a delayed emergence of a 
scaling region can be avoided by eliminating shear from the turbulence. The 
process of elimination entails creating a uniform mean profile, (U), and mini- 
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mizing the Reynolds shear stresses, (uv). The total shear stress being^] 

T(y)=pv £> ^-p(uv), (4.5) 

by having a uniform mean profile, the gradient term, d(U)/dy, can be kept 
small; and if the Reynolds shear stress, (u v), is small too, the total shear stress 
is then negligible. 

It was hoped that randomizing the forcing would suppress the growth of the 
Reynolds shear stress term, (uv). We gauged the Reynolds stresses by measur- 
ing the normalized correlations between orthogonal velocity components, 

«(o,o^X(o,o,z)) 

and 

«(r 1 ,0,0) li ' r2 (r ll 1 0)> 



zr2 v^olMMf ' 

at various locations, z and r lr along the axial and radial axes. Figure |4.5[ a) and 
(b) show tests at different levels of anisotropy. The shear stresses were always 
closer to zero than ±0.1 over the entire range of distances from the center of the 
soccer ball; |Shen and Warhaft| fl2000| ) and |Shen and Warhaftj fl2002| ) have mea- 



sured the shear stress in sheared turbulence and found that the value is about 
—0.4 and is constant throughout the flow. We concluded that the turbulence 
is essentially shearless because the mean shear stresses were less than 7% of 
the corresponding kinetic energies, u' Tl u' T2 and u' z u' r2 , respectively. The present 
forcing randomization scheme is a simple method that gave remarkably satis- 
factory shearless turbulence. We believe the shear can be reduced further, for 
example, by covering the nozzle with specially designed passive grid with so- 



It is customary to use U, V, and W to denote the mean velocities in the x, y, and z directions, 
respectively; and u, v, and w their corresponding RMS of the velocity fluctuations in the x, y, 
and z directions. 
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Figure 4.6: The figure shows the normalized Reynolds stress of orthogonal ve- 
locity components at various locations (a) along the axis of symmetry of the forc- 
ing, and (b) perpendicular to it. The values in the color bar are the anisotropy 
measured at the center of the soccer ball for various values of the anisotropy of 
the forcing. 
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lidit^cr = 0.34. For extreme values of solidity, o > 0.4, the resulting turbulence 
behind the grid is unstable and non-uniform^] (Corrsin 1963a[ ). Useful design 
guidelines regarding pressure losses across the grid, the uniformity of the tur- 
bulence behind the grid, and its stability are collected in|Roach|([T987|). 



4.7 Conclusions 

Our apparatus is unique among other loudspeaker-driven flow chambers be- 
cause it permits a systematic exploration of anisotropic turbulence. We have 
shown that we are able to control the large-scale anisotropy of the turbulence, 
as summarized by the ratio between axial and radial components of the root- 
mean-square velocity, to any value between 0.6 and 2.3, with 1 being the value 



for isotropic turbulence. Table 4.2 summarizes the values of various parame- 
ters for the flows. Because the anisotropy, the mean flow, and the strength of 
the fluctuations are approximately constant within 50 mm of the center of the 
soccer ball, we call this the region of homogeneity. 



2 The solidity of a grid, a, is defined as the projected solid area per unit total area. For a 
biplane grid with square mesh and round rods, the solidity is a = A (2 — 4r), where d is the 
rod diameter and M is the distance between rod centerlines, or mesh size. The popular grid 
solidity a = 0.34 is a compromise between low turbulence intensity v! / (U), and high flow 
stability. 

3 The critical solidity a = 0.4 limits the turbulence intensity behind the grid, since the tur- 
bulence intensity is proportional to the drag coefficient, and that, in turn, is proportional to the 
solidity. 
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.4 


0.33 


0.4 


0.5 


0.8 


1.0 


1.25 


2.0 


2.5 


u' z (m s" 1 ) 


0.74 


0.75 


0.79 


0.91 


1.08 


1.21 


1.42 


1.47 


K 2 ( m s_1 ) 


1.26 


1.24 


1.27 


1.19 


1.15 


1.03 


0.83 


0.75 


u'Ju' r2 


0.59 


0.60 


0.63 


0.77 


0.94 


1.17 


1.71 


1.97 


K (m 2 s- 2 ) 


1.87 


1.83 


1.92 


1.82 


1.91 


1.79 


1.69 


1.65 


\u z /u r2 \ 


0.31 


0.16 


0.13 


0.15 


4.04 


2.57 


7.58 


23.9 


\U\ (m s- 1 ) 


0.17 


0.18 


0.13 


0.10 


0.05 


0.07 


0.14 


0.24 


T rir2 


-0.07 


-0.06 


-0.06 


-0.07 


-0.07 


-0.05 


0.01 


0.04 




-0.06 


-0.07 


-0.06 


-0.07 


-0.07 


-0.06 


-0.01 


0.02 


e (m 2 s~ 3 ) 


5.41 


4.94 


5.43 


6.00 


6.71 


6.98 


6.49 


6.45 


V (H 


163 


167 


163 


159 


155 


153 


156 


156 


r,, (ms) 


1.7 


1.8 


1.7 


1.6 


1.5 


1.5 


1.6 


1.6 


A (mm) 


7.4 


7.6 


7.5 


6.9 


6.7 


6.3 


6.4 


6.3 


Rx 


525 


536 


538 


485 


480 


443 


434 


423 



Table 4.2: The table shows the turbulence parameters, for loudspeaker RMS 
amplitude ratios (A) ranging from 0.33 to 2.5. K is the turbulent kinetic energy, 



u 



+ 2< 2 ),and|?7|is (f/ 2 



2U 2 

r-2 



1/2 ? 



T rir2 and T zr2 



are the average normalized 



Reynolds stresses measured along and normal to the axis of symmetry, respec- 
tively (see section 4~6| ). e is the energy dissipation rate estimated according to 
the procedure prescribed in section 5.1 77 and r v are the Kolmogorov length and 
time scales, respectively A is the Taylor length (see section 5.3) and R\ is the 
Taylor scale Reynolds number. Data were collected at the center of the soccer 
ball when the axis of symmetry of the forcing is horizontal. Similar results (not 
shown) were obtained when the axis of symmetry is vertical. 
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CHAPTER 5 

UNIVERSALITY IN THE STRUCTURE FUNCTIONS 



In this chapter, we examine how the inertial scales of the turbulence gener- 
ated by the apparatus introduced in chapter [2] are influenced by the anisotropy 
of the forcing, through the measurements of the second order moment of ve- 
locity increments. We measured the transverse structure functions, D r2r . 2 (z) = 
(«,((), 0,z) -< 2 (0,0,0)) 2 ) and D zz {n) = («(n, 0, 0) - <(0, 0, 0)) 2 ), with the 
separation vector lying nearly along the axis of symmetry or perpendicular to 
it, and show that, in an inertial range that spanned nearly a decade, the Kol- 
mogorov constant and the inertial range scaling exponents were scalar valued 
in unsheared anisotropic turbulence, within the experimental error of 10%. 



5.1 Isotropy of the structure functions 



To determine the inertial range of the turbulence for different anisotropics, we 
plot in figures |5.1[ a) and (b) the measure for energy per unit time and mass, 

(D r2r2 (z))V 2 (4C 2 /3)- 3 / 2 z' 1 and {D zz {r x )fl 2 (4C 2 /3)- 3 / 2 as functions of z 
and t\, respectively. No noise corrections were applied to the structure func- 
tions. The value of C 2 was taken to be 2.1 ( Sreenivasan| 1995| ). We found that 



the inertial range spanned approximately one decade, as indicated by the extent 
of flat region in figures |5?lj a) and (b). 



Figure 5.2 shows the scale-dependent measure of isotropy, namely 



D zz (ri) / D r2r2 (z), where r\ = z. In the limit of large separations, the ratio should 
approach u' 2 /u l2 2 , since the velocities w-(r) and u-(0) are uncorrected when r is 
large enough. Although we could not resolve this asymptotic limit due to ex- 
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Figure 5.1: The figure shows the (a) axial structure functions, D r21 . 2 (z), and (b) 
radial structure functions, D zz (ri). In both cases, the structure functions are di- 
vided by the inertial range scaling predicted by Kolmogorov. The error bars are 
smaller than the symbols. The values in the legends are the anisotropy mea- 
sured at the center of the soccer ball for various values of the anisotropy of the 
forcing. 
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Figure 5.2: The ratios between structure functions in different directions but in 
the same flow, D zz (r)/D r2r2 (z), approach the isotropic value of one when the 
separation distance, r x = z, decreases. Results for data taken at other values of 
the large-scale anisotropy are consistent with those shown here. The values in 
the legends are the anisotropy measured at the center of the soccer ball. 



perimental limitations, the values of this ratio in anisotropic cases did clearly 
separate from the isotropic value. Moving toward small scales, local isotropy 
requires that D zz (r)/D r2T2 (z) approaches one. The figure reveals that within the 
experimental error we observed local isotropy for separations smaller than the 
energy injection scale. 



5.2 Equality of transverse ESS scaling exponents 



Here, we investigate how the scaling exponent measured using the ESS method, 



as described in section 1A\ depends on the anisotropy of the forcing. We begin 
with an examination of the third order transverse structure function, ((5u(x)) 3 ), 



and its ESS variant, (\5u(x)\ ). Figures 5.3 'a) and (b) show ((5u r2 (z)) s ) and 
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((5u z (ri)) 3 } as functions of z and r\, respectively. No noise corrections were ap- 
plied to these structure functions. It can be seen that ((5u r2 (z)) 3 } and ((5u z (ri)) 3 } 
were nearly zero and varied unsystematic ally with separations, as all odd- 



order moments should if the turbulence is locally isotropic (Shen and Warhaft 



2002[ ). This contrasts the predictions in shear flows, where simple Kolmogorov- 
type argument assuming that ((5u(x)) 3 ) is proportional to the shear yields 
((Su(x)) 3 ) ~ x 4 / 3 (Lumley 1967|. The absence of a four-thirds scaling reflects 



the absence of shear in our flow. Figures 5.4 'a) and (b) show the scaling of 
{\5u r2 (z)\ 3 ) and (\5u z (ri)\ 3 ) with separations, z and r\. It can be seen that the 
scaling range spans approximately a decade from 2 mm to 20 mm. For reference, 
it might be noted that the empirical scaling for (\5u(x)\ 3 ) is (\Su(x)\ 3 ) ~ x 1,17 for 
the present measurement, for which the value 1.17 has been obtained by aver- 
aging the exponents obtained from straight-line fits, \og(\5u(x) | 3 |) ~ log x, to the 
individual curves. 

We now consider the ESS plots for D r2T2 (z) ~ (\5u r . 2 (z)\ 3 )^ ' and D zz (ri) ~ 



(n) . 



(\Su z (ri)\ 3 )^ 11 in figure \5.5\ We obtain two sets of exponents, ^ for the axial 



function, D r2r . 2 (z), and £2 for the radial function, D zz (r^), by calculating the 
slopes obtained from straight-line fits, log(D(r)) ~ log(|£-u(r)| 3 ), to the data. 
The exponents, Q 2 ^ an d Q / are shown in figure 5.6 The uncertainty in the 



exponents was estimated with the standard error of the slopes in the straight- 
line fits (see equation B.12 in appendix [B]). Table 5.1 lists the numerical values 
for and (2 • The scatter in the data sets for D zz (ri) with anisotropy ratios 
u' z /u' r = 1.71 and 1.97 could be due to a number of reasons. It may have been an 
artifact of the limitation in instrument spatial resolution, a residual shear in the 
turbulence, or an insufficiently developed turbulent region. For these reasons, 



we have plotted their exponents with gray symbols in figure 5.6 The average 
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7*i (mm) 

(b) 



Figure 5.3: The third-order structure functions with separations in the direction 
(a) along the axis of symmetry, and (b) perpendicular to the axis of symmetry 
The values in the color bar are the anisotropy measured at the center of the 
soccer ball. 
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Figure 5.4: The third-order moments of the absolute values of the velocity in- 
crements with separations in the direction (a) along the axis of symmetry, and 
(b) perpendicular to the axis of symmetry The dashed lines are proportional to 
a power law with exponent 1.17. The values in the color bar are the anisotropy 
measured at the center of the soccer ball. 
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0.59 ±0.03 
0.60 ±0.02 
0.63 ±0.03 
0.77 ±0.03 
0.94 ± 0.04 
1.16 ±0.05 
1.71 ±0.10 
1.98 ±0.10 



0.739 ±0.007 
0.733 ±0.007 
0.733 ±0.006 
0.730 ± 0.004 
0.719 ±0.004 
0.752 ±0.018 
0.847 ±0.072 
0.875 ±0.101 



0.721 ±0.005 
0.732 ±0.008 
0.726 ±0.011 
0.727 ±0.004 
0.724 ±0.004 
0.733 ±0.005 
0.730 ±0.007 
0.738 ±0.007 



1.024 ±0.017 
1.001 ±0.020 
1.009 ±0.024 
1.004 ±0.011 
0.993 ±0.012 
1.026 ±0.031 
1.161 ±0.111 
1.186 ±0.147 



Table 5.1: The numerical data for radial scaling exponents (Q )/ axial scaling 
exponents (Q )/ and the ratio between the two. 



c 2 



Source 



0.70 ±0.005 37- 82 

0.70 ±0.01 365 -605 

0.71 1000 

0.71 - 0.74 40 - 10 4 

0.73 ±0.02 480 



Camussi et al. 


<1996|> 




Noullez et al.l 


1997 


Zhou et al. 


d20 


oil 


Pearson and Antonia 


(2001) 



Table 5.2: The value of the scaling exponent of transverse second order structure 
functions, (2, measured using the ESS method for various laboratory flows. 



value for the scaling exponents, excluding the last two radial data sets, was 
0.73 ± 0.02, for which the error was the the standard deviation of the values 
of the exponents. The value for the mean exponent is largely consistent with 
values measured using similar transverse ESS method reported in the literature 



at comparable Reynolds numbers, see table 5.2 



We now consider the ratio formed by the two exponents, C2/C2 ' an d ob 



serve how it varies with anisotropy, u' z /u' r . As shown in figure 5.7 except for the 



last two values of anisotropy, the ratio is nearly 1. Again, noise or residual shear 
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(|^ r2 (z)| 3 ) (|^(n)| 3 ) 

Figure 5.5: The extended-self-similarity plots for the second-order transverse 
structure functions. The dashed lines are proportional to the power laws ob- 
tained from least squares fits to data. The values in the color bar are the 
anisotropy measured at the center of the soccer ball for various values of the 
anisotropy of the forcing. The error bars are smaller than the symbols. 
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Figure 5.6: The scaling exponents for the transverse structure functions mea- 
sured using the extended-self-similarity method. The O symbols are the scaling 
exponents for D r2r2 (z) and the □ symbols are for D zz (ri). The last two points in 
the radial data set have large uncertainties due to spatial resolution. They are 
marked in gray to show the general trend of the exponents. 



might have affected the last two radial structure functions, even though the ex- 

(z) 

ponents Q f or the corresponding axial structure functions were well within the 
mean of the estimated values. Thus, based on our measurements for anisotropy 
ratios 0.59 < u'Ju' r < 1.16, we conclude that and £2 are equal. 



5.3 Scalar-valuedness of the Kolmogorov constant 



Let us now examine the Kolmogorov constant, C 2 (see equation 1.19 and 1.201. 



Anisotropy of the large scales allows for a family of equations in the form of 
equation 1.20 



but with different constants, Ci , for each direction of r. A 
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Figure 5.7: The figure shows (2 /Q , the ratio between the scaling exponents of 
the radial structure functions and those of the axial structure functions. The last 
two points with large uncertainties due to spatial resolution are marked in gray. 



measure of inertial range anisotropy is /C2 , which we calculated in two 



ways. First, in figure 5.8 we plot D zz (ri) / D r . 2r2 (z) against the separations r\ and 



z. These are the same data sets shown in figure 5.2 For clarity, they are dis- 



played individually. We computed an average of the structure function ratio, 
(D zz (ri) / D r2r2 (z)), taken over the range of separations where the compensated 
structure functions displayed a reasonable plateau, or 7 mm < r x , z < 70 mm. 
Since the energy dissipation, e, is a scalar quantity, the ratio of structure func- 
tions in the inertial range of scales is equal to the ratio /C^ • To estimate 
the uncertainty in the measurements, we have used the standard error in the 
measurements of the ratios (see appendix|Bj equation B.3l. The second method 
was designed to mimic the one usually employed to measure C 2 (or e) when 
the data are collected in only one direction. That is, we estimated C 2 from the 



98 



E- 

J- 



1.5 
1 

0.5 
1.5 

1 

0.5 
1.5 

1 

0.5 
1.5 

1 F 

0.5 
1.5 

1 

0.5 
1.5 

1 

0.5 
1.5 

1 

0.5 
1.5 

1 



0.5 



i4 



-S— H— n— 



I I 

^-fl-n- *, 



- a- 5— n= s= jj j-§-u— a— * - 





0.59 



0.60 



0.63 



0.77 



0.94 



1.17 



1.71 



1.97 



10' 



10' 



ri , z (mm) 



Figure 5.8: The figure shows the ratio between structure functions in different 
directions in the same flow, D zz {r\)/D r2r2 (z), for all values of the anisotropy of 



the forcing. The data are the same as in figure 5.2 They are plotted separately 
to display individual trends. The two vertical dashed lines demarcate the range 
of scales between 7 and 70 mm where the compensated structure functions dis- 
played a reasonable plateau. The values in the color bar are the fluctuations 
anisotropy measured at the center of the soccer ball for various values of the 
anisotropy of the forcing. 
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maxima of the compensated structure functions, D(x) x 2//3 . In order to reduce 
the influence of noise, we fit a second-order polynomial function 

y = a + b log(x) + c (log(x)) 2 , (5.1) 



to each data set in the range of scales between 2 and 50 mm, see figure 5.9 Al- 
though this form is decidedly inappropriate since the structure function should 
approach a constant for large separations, it proves to be very instructive in ana- 
lyzing the data. We estimated the maximum value of e 2//3 for each curve by 
calculating the ordinate of the peak of the parabola, given analytically in terms 
of the fitting parameters a, b, and c 

Z/max = a~ h — . (5.2) 

4c 

The ratio between the maxima, max(D zz (ri) r^ 2 ^ 3 ) / ma.x(D r2r2 (z) z~ 2 ^ 3 ), was 
then calculated. The ratio so formed, the energy dissipation rate being a scalar 
quantity falls out of the ratio, gave a second measure of Ci r) /C { 2 z) . To evaluate 
the fit, we calculated the difference between experimental and predicted values, 
& = Vi ~ Vii where yi is the original value from experiment and %)i is the value 
predicted by the parabolic fit. The standard deviation, or variability, of & 

N -| 1/2 

A" 



Variability 



1 - 



(5.3) 



where N is the number of data points used in the fit, gave an estimate for the 
error in the measurements 



The two measures of C^/Cjj^ are shown in figure 



5.10 



and their numerical 



values are listed side by side in table 5.3 It can be seen that within measurement 
error, equals C^; the small downward trend is comparable to the scatter. 
Thus, we conclude that the Kolmogorov constant is scalar valued over the range 
of anisotropics studied here. 
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10° 10 1 10 2 10° 10 1 10 2 

z (mm) T\ (mm) 



Figure 5.9: The compensated axial and radial transverse structure functions for 
all values of anisotropy. Only the data in the range of scales between 2 and 
50 mm are used in the determination of the ratio of the Kolmogorov constant. 
Solid lines are the second-order, least-square, linear-log polynomial fits to the 
structure function data. Solid symbols are the maxima of the fitted curves. The 
error bars are smaller than the symbols. The values in the color bar are the 
anisotropy measured at the center of the soccer ball for various values of the 
anisotropy of the forcing. 
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Figure 5.10: The ratio between Kolmogorov constants measured in different di- 



rections Ci r) /Ci z) . The open circular symbols (C) are the ratio between the max- 
imum values of the compensated radial and axial structure functions. The open 
square symbols (□) are (D zz (ri)/D r2r2 (z)), the ratio of radial to axial structure 
functions averaged over a range of scales 7 < r < 70 mm. 



Let us conclude with the calculation of the flow parameters. We deter- 
mined e by taking the peak values of the compensated structure functions, as 



plotted in figure 5.1 assuming that the second order structure functions obey 
Kolmogorov-type scaling law. This assumption may not be valid in view of 



the anomaly of the scaling exponent discussed in section 5.2 but dimensional 
considerations compelled us to adhere to this definition. The dissipation rates 
depended on the anisotropy and reached a maximum value of 6.71 m 2 s~ 3 for 
the case of isotropic turbulence. The variation of the dissipation rate with vary- 
ing anisotropy of the agitation is most likely the result of our decision to fix the 



total kinetic energy of the turbulence, see equation 370 Fixing the turbulence 



kinetic energy does not guarantee a constant energy dissipation rate, since the 
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< 

u' r 


/ D zz (n) \ 


Vjariahilitv 

V CI 1 IClL'lllL V 


2/3\ 

max(D zz /r 1 ' ) 


Variabilifv 

v ci i i ci kj in l y 


\Dr 2 r 2 (z)/ 


max(D r2r2 / z 2 / 3 ) 


0.59 ± 0.03 


0.998 


±0.020 


1.023 


±0.014 


0.60 ±0.02 


0.989 


±0.024 


1.022 


±0.014 


0.63 ±0.03 


1.004 


±0.023 


1.037 


±0.014 


0.77 ±0.03 


1.036 


±0.016 


1.036 


±0.013 


0.94 ± 0.04 


1.012 


±0.007 


1.002 


±0.009 


1.16 ±0.05 


0.990 


±0.014 


0.979 


±0.012 


1.71 ±0.10 


1.028 


±0.022 


1.003 


±0.023 


1.98 ±0.10 


1.024 


±0.029 


1.028 


±0.027 



Table 5.3: The numerical data for the ratio of the Kolmogorov constant for trans- 
verse structure functions measured in nearly orthogonal directions, computed 
from averaging the ratio D zz {r\) / ' D r2T2 (z) over a range of scales between 7 and 
70 mm (second column), as well as from taking the ratio between the maxima of 
appropriately compensated structure functions, max(D zz /r^ 3 )/ m&x(D r2r2 / z 2 ^) 
(fourth column). The values in the third and the last columns are the corre- 
sponding estimated error. 



energy dissipation rate depends also on the large scales of the flow, which in 
turn depend on the anisotropy 

Under isotropic forcing, the Taylor scale was A = y/lhvvP/e = 6.7 mm, 
where v = 1.568 x 10~ 5 m 2 s _1 is the kinematic viscosity of air at room tem- 
perature, and the Taylor-microscale Reynolds number was R\ = Xu'/u = 480. 
The Reynolds numbers for the anisotropic cases were slightly lower. The corre- 
sponding Kolmogorov length and time scales were rj = (z/ 3 /e) 1 / 4 = 155 |im and 
T n = \fvjl = 1.5 ms. 



103 



5.4 Conclusions 



We investigated systematically the influence of the anisotropic agitation on the 
inertial scales of turbulence in flows for which the ratio of axial to radial RMS ve- 
locity fluctuations was between 0.6 and 2.3, with Taylor based Reynolds number 
R\ = 480. For scales smaller than the energy injection scale, the second-order 
transverse velocity structure functions were independent of anisotropy. Within 
the experimental uncertainty, the inertial range scaling exponent and the Kol- 
mogorov constant were independent of anisotropy. We expect anisotropic cor- 



rections, such as those introduced by the SO(3) decomposition (Biferale and Pro- 



caccia 



2005 ), to be smaller than the error in our measurements (about 5%). This 



will be true unless the corrections would cancel in the two directions we mea- 
sured; an outcome we consider unlikely. Our findings contrast with previous 
results by Shen and Warhaft (2002 ), who found that the scaling exponent for the 
transverse structure function in the direction of a mean shear was different from 
the scaling exponent in the direction of the mean wind by about 0.1; this indi- 
cates that the anisotropy produced by shear may be inherently different from 
that produced by turbulent fluctuations. Our result is of practical importance in 
providing a means of unambiguously estimating the parameters of anisotropic 
flows without shear, such as the energy dissipation rate and the Reynolds num- 
ber. 
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CHAPTER 6 
HIGHER ORDER STATISTICS 



In this chapter, we extend our investigation of the influence of anisotropy on 
the higher order statistics. We provide measurements of the structure functions, 



(\5u(x)\ n ), for n — 4, 5, 6, as described in section 1.4 and the nondimensional 



ratios, {\5u(x)\ 4 )/(\6u(x)\ 2 } 2 and (\5u(x)\ 6 } / (\5u(x)\ 2 ) 3 . 

It has now become common in studies of intermittency to look at higher 
order structure functions. Measurements of higher order structure functions 
are very difficult. This is because higher order moments are influenced by rare 
events in the tails of the probability distribution (PDF), thus requiring a longer 
data set. We can estimate how much time and storage space it requires to collect 
these rare events. Experiments (e.g. Anselmet et al. 1984[) suggest that the prob- 



ability distribution of the velocity increments, 5u, can be approximated by an 
exponential, P(x) ~ e~ x , where x = 5u/ {\8u\ 2 ) l l 2 . For good convergence of the 
nth order moment, (\5u\ n ), the rule of thumb is to resolve Su up to n (\5u\ 2 ) l l 2 , 



which corresponds to the location of the peak of x n P(x) (e.g. Staicu and Van 



de Water, 2003). Let us consider the sixth order moment. From figure 16 of 



Shen and Warhaft| ( |2000[ ), we estimate that P(x = 6) » 10" 4 . If 10 counts are 
required around x — 6, then the total number of samples we need to collect 
is 10/10 -4 = 10 5 . This is very close to the number of samples we actually col- 
lected for Su in the experiment. This amount of data occupied approximately 
50 megabytes of disc space in 20 minutes. The total amount of time for 28 dif- 
ferent separation distances and 16 anisotropics (8 for scanning along the sym- 
metry axis of the forcing, and another 8 perpendicular to it) was then 6.2 days 
and the total disc space was about 22 gigabytes. Now for the 12th order mo- 
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ment, P(x = 12) « 10" 9 ( |Anselmet et al.| ([T984) r figure 3). Taking 10 again as a 



reasonable number of counts, the total amount of time and disc space required 
would be about 1700 years and 2.1 petabytes. Thus, it can be seen that a precise 
measurement of the high-order moments of the velocity increments with LDV 
may take many generations of graduate students to accomplish. 

6.1 Higher order structure functions 



In figures 6.1 6.2 and 6.3 we show the ESS plots of transverse structure func- 
tions up to the sixth order. As can be seen, the quality of data for higher order 
moments of velocity increments quickly deteriorates when increasing the order. 
This is especially so for extreme prolate shape anisotropics and for separations 
in the dissipative range. Let us examine the important issue of convergence. 



Figure 6.4 shows, from top to bottom, the normalized probability density of 
5u r2 (z), denoted as P(x), where x = 5u r2 (z) / (\5u r2 (z)\ 2 ) 1 ^ 2 is the velocity incre- 
ments normalized by their standard deviation, and P(x) multiplied by x 2 , x 3 , 
x 4 , x 5 , and x 6 for a nearly spherically symmetric turbulent flow, u' z /u' r = 0.94, at 
R\ = 480. Here, the probe volume spacing was 20 mm, which corresponded to 
about 130 T]. We obtained the PDF from approximately 78000 data points and re- 
solved the PDF to approximately 8 standard deviations. There seemed to be no 
convergence problems for the second to fifth moments. The tails of these curves 
tended to flatten out at high x. The sixth moment was not as well converged 
as the lower orders; yet all moments showed strong asymmetry. This is clearly 



indicated in table 6.1 by the difference in the area integrated under the curves. 
The values of the third to the sixth moments were —0.107, 6.88, —0.689, and 141. 
The same features can be seen in the inertial range PDF obtained by |Shen and 
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Order of moment (n) 





2 


3 


4 


5 


6 


(-l) n f° oo x n P(x) dx 


0.512 


0.519 


1.14 


3.57 


14.6 


71.8 


J °° x n P(x) dx 


0.488 


0.468 


1.03 


3.31 


13.9 


69.6 


x n P{x) dx 


1.00 


0.987 


-0.107 


6.88 


-0.689 


141 



Table 6.1: The numerical integral of the PDF and its moments from the second 



to the sixth order of an isotropically forced turbulent flow, shown in figure 6.4 



The isotropy ratio was u' z /u' r = 0.94 and the Reynolds number was R\ = 480 



(see table 4.2 1. 



Warhaft ( 2000[ ) in a wind tunnel homogeneous nonsheared turbulent flow at 



R\ = 934 and a probe spacing of 153 77. All their third to sixth moments showed 
strong asymmetry; their values were —0.46, 5.14, —7.81, and 66.7. Notice that 
the asymmetry in Shen and Warhaft ( 2000[ ) is left-right inverted because their 



definition of the velocity increments is 5u(z) = u(0) — u(z) whereas we have 
defined it as 8u(z) = u(z) — u(0). The difference in signs between our odd order 
moments and those of Shen and Warhaft ( 2000[ ) has been corrected in the above 



comparison. 



Figure |63 shows the PDF of 5u r2 (z) for an oblate shape forcing, u' z /u' r = 0.60, 
at the same probe separation and R\ comparable to the spherically symmetric 
case. We obtained about 10 5 points. The improvement brought about by an 
increase in the sample size to the convergence of the higher order moments 
was marginal; yet the asymmetry in the PDF was evident, as indicated by an 
increase in the values of the odd order moments in table 16.21 The values for the 
third to the sixth moments were —0.19, 7.07, —2.11 and 154. The asymmetry in 
the PDF of 5u was highly reproducible from flow to flow with different values 
of anisotropy. 
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Figure 6.1: The extended-self-similarity plots for the fourth-order transverse 
structure functions. Dashed lines are proportional to power laws obtained from 
least squares fits to data. The values in the color bar are the anisotropy measured 
at the center of the soccer ball. 
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Figure 6.2: The extended-self-similarity plots for the fifth-order transverse struc- 
ture functions. Dashed lines are proportional to power laws obtained from least 



squares fits to data. For legend, see figure 6.1 
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Figure 6.3: The extended-self-similarity plots for the sixth-order transverse 
structure functions. Dashed lines are proportional to power laws obtained from 



least squares fits to data. For legend, see figure 6.1 



110 




Figure 6.4: The probability density of velocity increments, 5u T2 {z), and their 
moments for an inertial range separation z/rj 130 and spherically symmetric 
forcing, u'Ju' r = 0.94, atR x = 480 (see table gg). Herex = <5M r2 (2)/(|5M r2 (^)| 2 ) 1 / 2 . 
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Figure 6.5: The probability density of velocity increments, 5u r2 (z), and their 
moments for an inertial range sep arat ion z/rj ~ 130 and oblate shape forcing, 
u'Ju' r = 0.60, at R x = 536 (see table 5^. Here x = 5u r2 (z)/ (\5u r2 (z)\ 2 ) 1 ' 2 - 
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Order of moment (n) 





2 


3 


4 


5 


6 


(-l) n J° oo x n P(x)dx 


0.541 


0.530 


1.17 


3.79 


16.2 


83.4 


J °° x n P(x) dx 


0.459 


0.430 


0.98 


3.27 


14.1 


71.5 


J^oo xn P{x) dx 


1.00 


0.96 


-0.19 


7.07 


-2.11 


154 



Table 6.2: The numerical integral of the PDF and its moments from the second to 
the sixth order of figure 6.5 The flow anisotropy was u' z /u' r = 0.60 at R\ = 536 



(see table 4.2 1. 



6.2 Kurtosis and hyperkurtosis of structure functions 



We now turn to the normalized structure functions. When normalized by the 
second-order structure functions, the fourth and sixth order moments are given 
special names, namely the kurtosis 

(IM*)I 4 > 



KAx) 



and the hyper-kurtosis 



(\6u(x)\ 2 ) 2 

(IM*)I 6 > 



(6.1) 



(6.2) 



(|^(a;)| 2 ) 3 ' 

At the small scales, K 4 > 3 and K 6 > 15 because of intermittency (Frisch 1995[ >. 
At the large scales the PDF of 5u follows a normal distribution and K± = 3 and 

K 6 = 15. 

The measurements for the normalized axial and radial structure functions, 



K 4 (z), K A {r\), K e (z), and K 6 (rx), are shown in figures 6.6 and 6.7 It is evident 
that our flows were highly intermittent at the small scales. Both K 4 and K 6 
showed trends of decreasing to the gaussian values of 3 and 15, respectively. 
Shen and Warhaf t| ( [2000] ) and |Garg and Warhaft| ( |1998| ) have shown that K 4 and 
Kq in sheared and nonsheared turbulent flows asymptote to the same limits at 
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the integral scale, approximately at 2000 77. If we can extend our measurement 
that far out, fT 4 and K 6 in our experiment may continue to decrease and even- 
tually reach their respective gaussian values. This is not physically possible in 
our experiment because 2000 77 would correspond to a distance of about 0.3 m, 
which is very close to the boundary of the apparatus where the anisotropy of the 
flow is heavily influenced by the confining walls rather than by the anisotropy 
of the agitation. What we gain by generating a stationary turbulent flow with 
negligible mean within confining walls we lose by having a finite extend in our 
turbulent region. Nevertheless, it is remarkable that K 4 and K 6 in our flow 
show downward trends toward their gaussian values similar to those observed 
in the wind tunnel. Shen and Warhaft (2000) also noted the presence of a small 
bump at 200 rj and attributed this to the scatter in the data. This bump was 
also noticeable in our measurements of fT 4 and K 6 , albeit at 100 r\. The magni- 
tude of this bump was as large as the scatter and may have evaded detection 
in experiments. We are only able to find another instance of the appearance of 
this bump at approximately 150 rj in the measurement of Siebert et al. ( 2010[ ) in 
the turbulent region of stratocumulus clouds at R\ rs 5000. We note that the 
scale at which the bump occurs approximately coincides with the Taylor scales, 



A = 15 1//4 R 1 ^ 2 1] (see e.g. Pope (20001, page 200, for a derivation of A), which 
were 40 rj, 60 rj, and 140 77 for R x = 480, 934, and 5000, respectively. From the nu- 
merical experiments of Donzis and Sreenivasan ( 2010[ ), we infer that the second 
order structure functions at R\ = 140 — 1000 reach a plateau at around the Taylor 
scales A = 20 rj — 60 rj. Thus, we speculate that the reason for the occurrence of 
this bump may be due to transition range dynamics. 
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Figure 6.6: The kurtoses as a function of separations, (a) The axial kur- 
toses, K 4 (z) = (\5u r2 (z)\ 4 ) / (\5u r2 (z)\ 2 ) 2 , and (b) the radial kurtoses K±{rx) = 
{\8u z (ri)\ A ) / {\5u z (ri)\ 2 ) 2 . The dashed lines are at the gaussian value of 3. The 
values in the color bar are the anisotropy measured at the center of the soccer 
ball. Error bars (not shown) are smaller than the symbols. 
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Figure 6.7: The hyper-kurtoses as a function of separations, (a) shows K 6 (z) = 
(|<^ 2 (z)| 6 )/(|<^ 2 (z)| 2 ) 3 and (b) shows K«{n) = (\6u g { ri )\«) / (\6u t { ri )\y . The 
dashed lines are at the gaussian value of 15. The values in the color bar are the 
anisotropy measured at the center of the soccer ball. Error bars (not shown) are 
smaller than the symbols. 
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6.3 Higher order scaling exponents 



Here, we document the inertial range scaling exponents of higher order struc- 
ture functions. The method for obtaining the higher order scaling exponents is 



the same as in section 5.2 We calculated the slopes obtained from straight-line 



fits, \og{\Su(x)\ n ) ~ \og{\Su(x)\ 3 ) , to the data shown in figures 6.1 6.2 and 6.3 



Tables 6.3 6.4 and 6.5 list the values for the axial and radial scaling exponents, 
Q z ' and (p \ for p = 4, 5, and 6. The radial structure functions for u' z /u' r = 1.71 
and 1.98 were very noisy and the scatter in these two data sets might have been 
due to a residual shear in the turbulence or an insufficiently developed turbu- 
lent region. The errors in all the exponents were estimated with the standard 
error of the slopes in the straight-line fits (see appendix [B]). Within the experi- 
mental error, we found that the values of the exponents are in fair agreement 
with those found in the literature (see section [6~4] >. 

We now examine the ratio of radial to axial exponents for order 4 and 5. In 
we plot Q r " '/Q and C^VCs^ against anisotropy, u' z /u' r . Be- 



figures 



6.8 



and 



6.9 



cause the majority of the ratios between radial and axial exponents are close 
to one, we conjecture that (p and (p Z ^ are equal in anisotropic turbulence. We 
note that Ouellette ( 2006[ ) measured the Lagrangian scaling exponents up to the 
tenth order in a von Karman counter-rotating flow and found that the axial and 
radial exponents are the same to within experimental accuracy. Thus, our in- 
vestigation here concerning the equality of Eulerian scaling exponents forms an 
extension to the above work. 
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0.59 ±0.03 
0.60 ±0.02 
0.63 ±0.03 
0.77 ±0.03 
0.94 ± 0.04 
1.16 ±0.05 
1.71 ±0.10 
1.98 ±0.10 



1.186 ±0.021 
1.189 ±0.027 
1.223 ±0.023 
1.213 ±0.012 
1.227 ±0.011 
1.098 ±0.049 
0.923 ± 0.486 
0.989 ±0.629 



1.220 ±0.017 
1.184 ±0.020 
1.173 ±0.038 
1.216 ±0.013 
1.215 ±0.015 
1.199 ±0.013 
1.231 ±0.022 
1.199 ±0.016 



0.972 ±0.031 
1.005 ±0.040 
1.043 ±0.053 
0.998 ±0.021 
1.010 ±0.021 
0.916 ±0.051 
0.750 ± 0.400 
0.824 ± 0.535 



Table 6.3: The values of the fourth order radial and axial scaling exponents, Q 
and ([ z \ measured in our experiment using ESS. The last column shows the 
ratio between the two. 



u'm cf d 2) ti ] /ti z 



0.59 ±0.03 
0.60 ±0.02 
0.63 ±0.03 
0.77 ±0.03 
0.94 ± 0.04 
1.16 ±0.05 
1.71 ±0.10 
1.98 ±0.10 



1.338 ±0.052 
1.327 ±0.072 
1.441 ±0.061 
1.394 ±0.032 
1.429 ±0.029 
1.118 ±0.099 
0.923 ±0.852 
0.989 ± 1.102 



1.407 ±0.043 
1.312 ±0.049 
1.290 ±0.082 
1.403 ±0.036 
1.396 ±0.039 
1.355 ±0.036 
1.458 ± 0.054 
1.373 ±0.038 



0.951 ±0.066 
1.011 ±0.092 
1.117 ± 0.118 
0.993 ± 0.048 
1.023 ±0.050 
0.825 ±0.095 
0.633 ±0.608 
0.720 ± 0.822 



Table 6.4: The numerical data for the fifth order radial and axial scaling expo- 
nents, (s \ £5, and the ratio between the two. 



6.4 Anomalous scaling 



We now turn to the issue of how the transverse scaling exponents, ( p , vary with 
order, p. In the absence of intermittency, Kolmogorov scaling predicts a linear 



relationship between ( p and p, see equation 1.49 Figure 6.10 and table 6.6 show 



the mean of C P in tables 5.1 6.3 6.4 and 6.5 as a function of p. The error bars were 
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0.59 ±0.03 
0.60 ±0.02 
0.63 ±0.03 
0.77 ±0.03 
0.94 ± 0.04 
1.16 ±0.05 
1.71 ±0.10 
1.98 ±0.10 



1.490 ±0.085 
1.447 ±0.120 
1.668 ±0.104 
1.562 ±0.055 
1.630 ±0.051 
1.131 ±0.144 
0.583 ±0.170 
0.426 ±0.175 



1.584 ±0.076 
1.424 ±0.082 
1.395 ±0.127 
1.584 ±0.067 
1.570 ±0.069 
1.499 ±0.063 
1.696 ±0.090 
1.544 ±0.062 



0.941 ±0.099 
1.017 ±0.143 
1.196 ±0.183 
0.986 ±0.076 
1.038 ±0.078 
0.754 ±0.128 
0.344 ±0.119 
0.276 ±0.124 



Table 6.5: The numerical data for the sixth order radial and axial scaling expo- 
nents, Q r \ £g , and the ratio between the two. 




Figure 6.8: The figure shows C4 /Ci / the ratio between the fourth-order scaling 
exponents of the radial structure functions and those of the axial structure func- 
tions. The last two points with large uncertainties due to spatial resolution are 
marked in gray. 
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Figure 6.9: The figure shows Cs^/Cs^/ the ratio between the fifth-order scaling 
exponents of the radial structure functions and those of the axial structure func- 
tions. The last two points with large uncertainties due to spatial resolution are 
marked in gray. 



the standard deviation of the values of ( p in each table. Our observation is that, 
within this error, the exponents show an anomalous scaling (departure from 
Kolmogorov prediction), similar to what is observed for the scaling exponents 
measured in a nonsheared wind tunnel homogeneous and isotropic turbulent 
flow at an R\ of 134 ( Shen and Warhaftj 2002[ ). The values of both sets of expo- 
nents are below the values obtained by |Shen and Warhaft ( 2002[ ) at an R\ of 863. 
The difference is less than 6% at the fourth order and increases monotonically 
up to about 13% at the sixth order. We suspect that the reason for the differ- 
ence is a combination of Reynolds number and large scale anisotropy effects. 
We believe that if we could achieve higher Reynolds numbers, the values of the 
scaling exponents would increase. 
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0.5^ ' ' ' ' ' — 
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V 

Figure 6.10: The circular symbols (O) are our measurements of the inertial range 
scaling exponents, C, p , of the transverse structure functions using ESS method. 
The triangular symbols (A) are the exponents measured by IShen and Warhaft 
(2002) in an unsheared wind tunnel homogeneous and isotropic turbulent flow 
at an R x of 134, and the square symbols (□) are those measured by them at an 
R\ of 863. The dashed line is the Kolmogorov prediction (equation 1.491. The 



solid line is the She-Leveque prediction (equation 1.54). 



6.5 Conclusions 



We studied the influence of the large scale anisotropy on the structure func- 
tions of order 4, 5, and 6. We found indications that anisotropy in the velocity 
field intensifies the asymmetry of the probability density of the velocity incre- 
ments. We have shown that the small scales of the flows are highly intermittent, 
as indicated by the measurements of the kurtosis and hyper-kurtosis, and that 
both kurtosis and hyper-kurtosis asymptote to their gaussian limits at the large 
scales. Within the experimental uncertainty, we found evidence that the scal- 
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Order p 




R\ 






480 


134 


863 




Present work 


S&W (2002) 


S&W (2002) 


2 


0.731 ±0.017 


0.68 ±0.05 


0.69 ±0.05 


3 


1 


1.00 ±0.05 


0.98 ±0.05 


4 


1.198 ±0.088 


1.17 ±0.05 


1.24 ±0.05 


5 


1.360 ±0.167 


1.37 ±0.05 


1.49 ±0.05 


6 


1.516 ±0.096 


1.53 ±0.05 


1.71 ±0.05 



Table 6.6: The inertial range transverse scaling exponents, ( p , measured in the 



current experiment using the ESS method and those measured by Shen and 



Warhaft (2002), denoted as S&W (2002) above, in an unsheared wind tunnel 



homogeneous and isotropic turbulent flow at two different R\ of 134 and 863. 



ing exponents measured using the ESS method are independent of anisotropy. 
The exponents showed a departure from the Kolmogorov prediction and their 
values were below those measured in a nonsheared wind tunnel homogeneous 
turbulent flow at an R x of 863 ( |Shen and Warhaftl |2002| ). This is a subtle dif- 
ference and we believe the difference would become smaller if we increase the 
Reynolds number. 
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CHAPTER 7 
THE INTEGRAL SCALES OF TURBULENCE 



Recently, there is increasing evidence to suggest that macroscopic properties 
of a turbulent flow are linked to the turbulence energy spectrum. Some of the 
macroscopic properties considered include the cross-sectional mean velocity of 
gravity-driven open channel flows (Gioia and Bombardelli 20011 ), the depth of 
turbulent cauldron (Gioia and Bombardelli, 2005), the friction factor ( |Gioia and 



Chakrabort^[2TO)[Guttenberg and Goldenfeldj [2009) |Tran et al.}|20l0) , and the 



mean- velocity profile of pipe flows (Gioia et al. 2010[ ). The authors established 
their findings based on the assumption that Kolmogorov's theory for homoge- 
neous and isotropic turbulence is also applicable to anisotropic and inhomoge- 
neous flows ( |Knight and Sirovich| |1990f |Moser| |1994f |Lundgren| |2003| ). They 
explain their findings in terms of the varying habits of momentum transfer with 
varying sizes of turbulent eddies. Here, we study the physical space equiva- 
lent of the spectrum - the correlation of velocity fluctuations separated in space. 
Our results highlight the interesting feature that correlation functions measured 
in different flow directions collapse onto one self-similar curve when appropri- 
ately scaled with the anisotropy that drives the turbulence at the large scale. 
The scaling form allows us to express the integral length, which is represen- 
tative of the length scale of the most energetic eddies (Batchelor ( 1956[ ), chap- 
ter 6), as a power law of its characteristic velocity, with an exponent related to 
the inertial range scaling exponent. Our work complements the existing stud- 
ies and suggests that self-similarity and Kolmogorov theory are relevant to the 
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phenomenological description of the large scale of anisotropic turbulent flows. 



7.1 The integral scale and correlation function 



We begin with a brief introduction to the integral length of the correlation func- 
tion in order to facilitate the discussion that follows. 

As described in the introductory chapter, an important tenet in the Richard- 
son picture of the turbulence energy cascade is that the energy dissipation rate 



per unit mass in a turbulent flow is finite (see e.g. Frisch (1995), 5.2). This is 



summarized in the following empirical law (Taylor] 1935[ > 

Au' 3 



(7.1) 



Here, A is a number of order unity, u' is a velocity scale characteristic of the 
most energetic eddies, and L their characteristic length scales. This order 
of magnitude assumption has found support from a theoretical consideration 
by Lohse (1994). Starting from an energy balance equation derived from the 
Navier-Stokes equation with a mean field approach, Lohse (1994]) showed that 
A m (a/C 2 ) 3 / 2 , where C 2 is the Kolmogorov constant for the second-order struc- 



ture function (see equations 1.19| and 1.20) and a ~ 1.25. Recent updates on 



the empirical law of energy dissipation rate can be found in Sreenivasan ( 1984[ ) 
and Pearson et al. ( 2002[ ). A simple interpretation of the above empirical law is 
that the decay time of the turbulence is a few characteristic periods of the most 
energetic eddies (Batchelor (1956), 6.1). Identifying u' with the RMS velocity 
fluctuation, u' 2 /e can be regarded as the time scale of the decay of energy, and 
L/u' can be regarded as the characteristic period of the most energetic eddies. 

The above empirical relation has found wide use in estimating the energy 
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dissipation rate of turbulent flows. Measuring vl is relatively simple. Single- 
point measurements of the turbulent velocity fluctuation suffice to yield u' . In 
contrast to u' , measuring L is not so straightforward because it is not directly 
accessible (see e.g. Batchelor ( 19561 ), appendix to 6.1). If we consider the three- 
dimensional energy spectrum, &(k), it is permissible to think of the part of 
the function around the maximum, at wavenumber k , as responsible for sup- 
plying energy into the flow (the most energetic eddies). The reciprocal of the 
wavenumber at which this maximum occurs, Le = 1/ko, then characterizes 
the principal energy bearing part of the spectrum (the most energetic eddies). 
On the other hand, if the contribution from higher wavenumbers is negligible, 



we may approximate L with the transverse integral length (see e.g. Batchelor 

< [T956| )) 

_ 3 7T /^K-^fiQdK 

8 J Q <S>{k)&k 

which, for isotropic turbulence, is related to the one-dimentional transverse cor- 
relation function in the following way 



L g = / g(r)dr. (7.3) 
Jo 

Here, g(r) is the diagonal element of the much more general correlation tensor 

9ij{r) = , = , (7.4) 

(uf{r)) «(0)) 



which is the velocity correlation tensor introduced in section 1.1 normalized by 
the fluctuations. The careful reader may note that, in general, the velocity fluctu- 
ations are functions of space and time. The time dependence of the fluctuations, 
and consequently gij{r), has been suppressed for the simplicity of discussion. 
If the velocity field satisfies the homogeneity and isotropy conditions, then the 



correlation tensor can be expressed in terms of two scalar functions (Batchelor 
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([1956), 3.4). With reference to experimental work, we usually make the choice 
r = re\ and introduce the longitudinal correlation function 



f(r) = 9u(r) 



(7.5) 



and the transverse correlation functions 



g(r) = 922{r) = gs 3 (r) . 



(7.6) 



In this way, we obtain the form for the correlation tensor first derived by 
Karman and Howarth ( 1938[ ) and systematically formulated by Robertson ( 1940[ > 



9a(r) = g(r)6ij + [f(r)-g(r)] 



(7.7) 



f(r) and g(r) are related by the continuity equation in much the same way as the 
longitudinal and transverse structure functions are (e.g. Pope[ 2000[ ). Relevant 
to our experiments are the following two transverse correlation functions 

« 2 (0,0,zK 2 (0,0,0)> 



g r2r2 (0,0,z) 

9zz (n,o,o) 



v/«(o,o,z))«(o,o,o)) 

«(r 1; 0,0) <(0,0,0)) 



(7.8) 
(7.9) 



V(«?(ri,0,0)>( U «(0,0,0)>' 
They are related simply to the structure functions by expanding the product in 
equation 1.30 

D(x) 



(7.10) 



(u'*(0)) \/ <u' 2 (*)> y/(u»(x))(u*(0))\ ' 
which, under the assumption of homogeneity (u' 2 (x)} = (w /2 (0)) = u' 2 , reduces 



to 



9{x) 



D(x) 



2u 



12 



(7.11) 



Here, x is understood to represent z or r\. Given that L g is an important control 
parameter in aerodynamic stability analyses of bridges with sharp edges (e.g. 
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Fransos and Bruno 2010 ), as well as being a key parameter in the design of gas 



turbine engines (e.g. Van Fossen and Ching| 1994; Van Fossen et al. 1995 Barrett 



and H ollin gsworth[ [20011 |Carullo et a!} |20TT| ), a more fundamental understand- 



ing of L g , and thus of g(r), is needed. 



7.2 Scaling and self -similarity 



There is no rigorous theory that yields the shape for the longitudinal and trans- 
verse correlation functions in fully developed turbulence, but we have a rough 



picture of how these curves should look like from experiments (e.g. Comte 



Bellot and Corrsin 1971[ ) and numerical simulations (e.g. Kim et al. 1987[ >. 
Batchelor and Townsend (1948]) studied the decay of turbulence and showed 
that in the final period of decay of the turbulence, when the effects of inertial 
forces are negligible, the longitudinal correlation function asymptotes to a self- 
similar function 



f(r,t) = exp - 



8ut 



(7.12) 



corroborating the hypothesis made by |Karman and Howarth ( 1938 ) (for more 
details on the self-preservation hypothesis, see e.g. Monin and Yaglom ( 1975) , 
16). Temporally self-similar correlation functions were also derived analytically 



in decay problems of Burgers turbulence (Gurbatov et al. 1997) and passive 
scalar in Kraichnan's model ( Eyink and Xin| 2000[ ). |Ewing et al. ( 2007[ ) exam- 
ined the correlation function of streamwise velocity component in the far field 
of an axisymmetric jet and showed that the similarity solutions of the govern- 
ing equations for the correlation functions depend only on the separation dis- 
tance between the points in the streamwise similarity coordinate, namely ip' — ip, 
where the similarity coordinate ip = In x is measured from a virtual origin in 
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space x. The spatial self-similarity of correlation functions in anisotropic turbu- 
lence considered here complements the existing works. 

The similarity argument to be presented here is a special case of Widom's 
scaling (see e.g. Huang ( 19871 ), 16.5) with the dimension of homogeneity d 
equals zero (E yinkJ 2011), and is in the same vein as Goldenfeld (2006). Two 
notable features of this approach are the universal scaling function that gov- 
erns the correlation function and the reduction of the number of variables in the 
function. The scaling form we propose for the transverse correlation function is 



g(x,u') =y{xu' p ) 



(7.13) 



where £f(£) is a universal scaling function of a single variable £, and is an 
exponent to be determined. To determine 0, we connect the scaling function 
to that in the inertial range. We note that g(x, u') scales as 



g(x,u) oc 



u 



(7.14) 



in the inertial range. This can be seen by substituting the inertial range scaling of 
the structure function, D(x) oc x^, into equation 7.11 This requires that oc 



in the inertial range, and therefore = — 2/£. Thus, the scaling form for the 
correlation function is 

>'- 2 K) . (7.15) 



g[x,u 



[xu 



The scaling form predicts that correlation functions will collapse onto a single 
curve when plotted as g against x u'~ 2 ^. 



The scaling form in 7.15 implies a relationship between the integral length 
and the large scale velocity fluctuation. By changing the variable of integration 

to £ = x u'^ 2 ^, we obtain 



/2/C 



(7.16) 



128 



where J? = f °° d£ is a constant which may depend on some other parame- 
ters of the flow, e.g. R\ or (. This dependence can be eliminated if we form the 
ratio of two integral scales obtained from two different correlation functions 



i , , • ( 7 - 17 ) 

L 2 \u 2 J 

We cannot test these predictions in isotropic turbulence because the correlation 
functions and RMS velocity fluctuation are isotropic by definition. If we can 
drive the flow away from isotropy, as we can in our flow generator, we can test 
these predictions. 



7.3 Integral length scaling 



Figure 7.1 shows 8 pairs of correlation functions, g r2r2 (z) and g zz (ri), obtained 
in 8 different flows with different large scale anisotropics, but with approxi- 
mately the same energy dissipation rate and Reynolds number. We found that 
correlation functions measured in different directions in anisotropic flows were 
different when plotted against separation. Note that we normalized the corre- 
lation functions by the velocity fluctuations at both points, and r\ (or z), al- 
though homogeneity would allow us to set them equal, because it compensated 
for the small inhomogeneity of our flow. In addition, we divided the correlation 
functions by their values at zero separation, that is, when the distance between 
probe volumes was zero. The values of the correlations at zero separation were 
not equal to one because of noise in the signals, but were approximately equal 
to 0.98, for all correlation functions. 

In order to test equation 7.17 we need to estimate the integral scale. |Batch- 



elor ( 1948[ ) observed that transverse correlation functions cross zero before they 
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Figure 7.1: The figure shows the 8 pairs of transverse correlation functions ob- 
tained in 8 different flows with varying degree of large scale anisotropy. Circles 
(O) are for g r2 r 2 (z) and crosses (x) are g Z z{ri)- The error bars are smaller than the 
symbols. The values in the color bar are the fluctuations anisotropy measured 
at the center of the soccer ball for various values of the anisotropy of the forcing. 
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return to zero asymptotically, but it is not unusual that experimental measure- 
ments do not resolve this limit because of the finite size of experiments (e.g. 
Hwang and Eaton} |2004} |de Jong et aTj [2009| |Siebert et alj |2010) . In fact, the 



first zero crossing would occur outside the region of our flow that is approxi- 
mately homogeneous. On the other hand, Comte-Bellot and Corrsin ( 1971[ ) ar- 
gued that integral scales in finite flows should strictly be zero, if data could be 
collected that extended far enough. This can be seen by considering the case 
k = in the Wiener-Khinchin relation between the one-dimensional transverse 



spectrum, E 22 (k), and its correlation function, #22(7*) (see equations 1.27 and 

vr£ 22 (0) 



1.28), which yields 



2 (u' 2 ) 



(7.18) 



If we consider, however, the spectrum given by E 2 2(k) = u(k) u*(k), where the 
Fourier coefficients of the velocity field averaged over a distance L in physical 
space are u(k) = \ f Q L u'(r) e~ lkr dr, given that the mean velocity f Q L u'(r) dr is 
zero, it follows that the coefficient of the zeroth Fourier mode is zero, u(0) = 0. 
It can be immediately seen that £^22(0) = 0, and thus L g = 0. Only hypotheti- 



cal infinite flows could have nonzero integral scales. Comte-Bellot and Corrsin 



( |1971 ) note that the only way to construct a nonzero integral scale from real data 
is to make 'simple extrapolations' of partially measured correlation functions in 
a reasonable way. We chose to extrapolate the correlation functions with ex- 
ponential functions (e.g. Van Fossen and Ching} 1994} Van Fossen et all} 1994} 
de Jong et al} [2009) . We are aware that physically plausible extrapolations of 



the transverse correlation functions must include zero crossing, but Lenschow 



and Stankov (1986) found that, following the suggestion of Comte-Bellot and 



Corrsin ( 1971[ ), approximating the one-dimensional longitudinal spectrum (see 



equation 1.29) with a downward parabola with zero slope at zero wavenum- 
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ber gives a correlation function with oscillatory tail. Integrating this correlation 
function up to its first zero-crossing yields a value that overestimates the inte- 
gral length by no more than 5%. We found empirically that using more elaborate 
functions, such as the modified Bessel function, for the extrapolations did not 
alter our results for the integral length scaling exponent. The difference between 
the exponents calculated using different fitting functions is less than 10%. We 
suggest that this is because we only consider ratios of integral scales, not the in- 
tegral scales themselves, and the influence of the two neglected negative parts 



of the correlation functions cancel each other. Figure 7.2 shows how a typical 
extrapolation is done. For each correlation function, we fit an exponential using 
a least squares algorithm (e.g. Draper and Smith 1998[ ) to the last 15 data points 
collected between < x < x max , where the maximum separation, x max , ranged 
from 8 mm to 120 mm. We then chose the value of the maximum separation, 
x max/ that minimized x 2 > the sum of squares of the vertical differences between 
the experimental, g(xi), and the fitted values of the correlation function g(xi) 

N 

^ = X>fo) - S(*i)) 2 , (7.19) 
i=i 

where iV is the total number of samples used in the fit. Figure |7.3| shows the 
variation of x 2 with x max . It can be inferred that x max = 13 mm yields the fit with 
the smallest sum of squares of residuals. Data for separations x > x max were 
extrapolated with exponential tails, whose integral we added to the numerical 
integration of the data, to yield the integral length. The values of the integral 
scales did depend on a; max , leading to a variation in the value of the exponent 
(see figure [7~4] >. We then used the standard error (see appendix [B]) of value of the 
integral scale calculated with different x max to estimate the error of our integral 
scale measurement. 
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Figure 7.2: The extrapolation of the partially measured correlation function. The 
blue circular symbols (O) are the experimental data. The dashed line is the 
exponential extrapolated from the data. The area in the region shaded in blue is 
numerically integrated and the area in the region shaded in pink is analytically 
calculated. The sum of the two yields the integral length. 



We note that, since u' { is a transverse velocity the subscript in u\ indicates 
the direction of the velocity and the subscript in L t indicates the direction of the 
separation over which the integral length is calculated. Therefore, the expected 
scaling is a power law with a negative exponent 



XL 



u: 



-2/C 



(7.20) 



in conformity with our definition of axial and radial directions. We test this 



relationship in figure 7.5 Note that we have used the definition given in equa- 
tion 



4.2 for u' z /u' r , in compliance with the definition of v! in 7.17 The integral 
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10 2 



x max (mm) 

Figure 7.3: The sum of squares of residuals in the determination of x max that 
yields the exponential from a least squares fit of an exponential tail to the corre- 
lation function data. 



length scaling exponent obtained from a least squares straight-line fit 

(7.21) 



L, W 



L r V it 



to the data is (5 = —2.56 ± 0.07, for which the uncertainty was estimated with 
the standard error of the slope in the straight-line fit (see appendix |B|). To check 
for consistency, we calculated the value for (3 with the mean value of the iner- 



tial range scaling exponents measured using ESS method, £ 2 (see section 5.2 1, 
for which we find (3 2 = —2/(2 = —2.74 ± 0.06, assuming that the percentage 
error remains the same (A/3 2 //3 2 = A£ 2 /£ 2 ). The value for (3 2 is no more than 
7% greater than (3. Minor differences between the two scaling exponents may 
arise from the difference between the exponential and the real correlation curve, 
which is at present unknown. A function that takes into account the transition 
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-2.5 




Smax (mm) 

Figure 7.4: The variation in the value of the exponent with the maximum sepa- 
ration in the least squares fit of exponential tails to the data. 

range between the inertial scales and the large scales will be able to account for 
the difference. Nevertheless, even with these simple methods, we find excellent 
agreement between the two exponents. 

In addition, we note that the value of the exponent showed a departure from 
the Kolmogorov prediction, which has a value of —3 for (3 if the value of ( is 2/3. 
The difference is very likely due to a Reynolds number effect. We suspect that 
the difference would become smaller if we increase the Reynolds number. 
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Figure 7.5: The main figure shows the ratio of integral scales, L z / L r , as a func- 
tion of the anisotropy in the fluctuating velocities, u' z /u' r . The solid line is the 
least squares fit to the data and the dashed line is the power law with exponent 
—3 predicted by the Kolmogorov inertial range scaling argument. The inset 
shows the same data on a logarithmic scale. The straight line in the inset has a 
slope of —2.56. 



u' z /u' r L z (mm) L r (mm) L z /L r 



0.67 ±0.03 
0.69 ±0.03 
0.71 ±0.03 
0.84 ±0.03 
1.00 ±0.04 
1.16±0.03 
1.58 ±0.09 
1.76 ±0.09 



100 ±4.3 
95.6±4.8 

93.0 ±4.1 
76.9 ±3.3 

67.5 ±2.5 
61.0±1.9 

45.1 ±1.0 

42.6 ±0.5 



33.5 ±0.6 
34.2 ±0.9 
36.1 ±1.0 
50.2±1.9 
68.0±2.1 
95.4 ±2.6 
144 ±3.9 
155±4.5 



2.99 ±0.18 
2.80 ±0.21 
2.58 ±0.18 
1.53±0.12 
0.99 ±0.07 
0.64 ±0.04 
0.31 ±0.02 
0.28 ±0.01 



Table 7.1: The numerical data for axial integral lengths (L z ), radial integral 
lengths (L r ), and the ratio between the two. 
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7.4 Correlation functions collapse 



The test of data collapse is shown in figure 7.6 The value for the inertial range 



scaling exponent is £ = 0.73, obtained using the ESS method (see section 5.2 1. It 
can be seen that correlation functions stretching out in different directions on the 
plane in figure [7T| collapse onto one single curve. Small deviations from the data 
collapse are visible, especially at large separations, which may be a reflection of 
the inhomogeneity occurring outside the turbulent boundary, where the effects 
of individual jets are felt. 



7.5 Conclusions 



We have presented, in the foregoing pages, a similarity argument applicable to 
the correlation functions whose main predictions is the power-law scaling of 
the ratio of integral lengths with the large scale anisotropy ratio, with an ex- 
ponent governed by the inertial range scaling exponent. Within experimental 
accuracy, the scaling exponent determined by direct measurement of the ratio 
of integral lengths and the large scale anisotropy ratio was found to be in excel- 
lent agreement with the value derived from the inertial range scaling exponent 
determined with ESS method. Further, we presented evidence to support the 
claim that correlation functions can be made to collapse when separations were 
rescaled with a scaling factor derived from the similarity hypothesis. 
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Figure 7.6: The figure shows the same correlation functions from figure [7T| plot- 
ted against the reduced variables, T\ u' z ~ 2 ^ and zu^ 2 ^. Here, ( = 0.73. Circles 
(O) are for g r2T2 (z) and crosses (x ) are for g zz {ri). The values in the color bar are 
the fluctuations anisotropy measured at the center of the soccer ball. 



138 



CHAPTER 8 
SUMMARY AND OUTLOOK 



In this dissertation, we have succeeded in isolating the effect of shear from 
the large scale anisotropy. We saw that anisotropy in the velocity field has dif- 
ference effects on the inertial and integral scales of turbulence. 

In chapter|4]we demonstrated that asymmetry of forcing leads to asymmetry 
of turbulence in a controllable way We saw that the flows are approximately 
homogeneous and axisymmetric, have negligible mean and shear. 

In chapter |5j we investigated second-order moments of velocity increments 
measured in different directions in anisotropic turbulent flows. We ruled out 
the possibility that the scaling exponent and the Kolmogorov constant might 
depend on the direction in which the structure function is measured. This is 
in contrast with previously published results and suggests that the anisotropy 
produced by shear may be inherently different from that produced only in the 
fluctuations. 

In chapter [6j we examined moments of the velocity increments up to the 
sixth order. We found evidence that the asymmetry in the probability density of 
the velocity increments is enhanced when the large scale anisotropy is present. 
In addition, we found indications that the inertial range scaling exponents are 
very likely to be independent of anisotropy. 

In chapter |7j we presented a scaling argument to describe the anisotropy 
observed at the large scale. We saw that correlation functions measured in dif- 
ferent directions of the flow collapse onto a single curve when appropriately 
scaled. The scaling implies a power-law relationship between the ratio of in- 
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tegral lengths and velocity fluctuation anisotropy, whose exponent is closely 
linked to the inertial range scaling exponent. We found that the exponent mea- 
sured from the power-law relationship is consistent with the value estimated 
from the inertial range scaling exponent with extended-self-similarity method. 
This power-law relationship suggests that self-similarity and Kolmogorov the- 
ory are relevant to the large-scale phenomenology of turbulence. 

It is important that the reader is aware of some of the difficulties and short- 
comings of the scaling argument presented. A point of vagueness is the com- 
plete neglect of the dissipation range; correlation function in this range has 
been shown to follow an analytic scaling law. The second point is the observed 
anomaly and intermittency in the inertial range scaling exponent. These ques- 
tions will not be settled by the present measurements, but these measurements 
will provide more information on which to base a future study. 

Further insight into these questions may be provided with tools like the 



particle image velocimetry (e.g. Raff el et al. 2007) and the Lagrangian parti 



cle tracking (e.g. |Ouellette] 2006). These measurement systems would allow 



us to probe the small scales with even higher precision at higher data rate. If 
further experiments were to be conducted, then we should investigate the an- 
gular dependence of the longitudinal structure functions and correlation func- 
tions. Do they follow the same scaling law? Or are they intrinsically differ- 
ent from transverse functions, as measurements in shear flows suggest? In ad- 
dition, we should examine the angular dependence of higher-order structure 
functions with longer measurement time, derivative structure functions, and 
the Lagrangian structure functions in our system. 
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To conclude, let us discuss a fascinating extension of the present work to tur- 
bulence in dimensions other than three. A well-founded scientific proof must 
explore and exhaust all the parameter space in the problem. A parameter we 
have not been able to vary is the inertial range scaling exponent, ( 2/ because 
our system is inherently three-dimensional and it obeys the three-dimensional 



forced Navier-Stokes equation. Two-dimensional turbulence (e.g. Kraichnan 



and Montgomery! 1980[ ) has captured the fascination of the researchers because 



many geophysical and magneto-hydrodynamical phenomena can be formu- 
lated as a two-dimensional problem. In the inertial range of two-dimensional 
turbulence, Kolmogorov-type dimensional reasoning gave r 2 scaling for the 



second-order velocity structure functions (e.g. Kellay and Goldburg 2002). The 
question would then be, does the ratio of integral lengths scale with the fluc- 
tuations anisotropy ratio following a power-law with a different inertial range 



scaling exponent? Measurements in thin layers of conducting fluid (e.g. Xia 



et al. 2011 ) and in soap films (e.g. Kellay and Goldburg, 2002[ ) will provide the 



answer to this question. 

We have presented in this thesis a flow apparatus that permits a system- 
atic exploration of turbulence away from isotropy and described some unique 
features of this system through measurements of structure functions and corre- 
lation functions. We believe that understanding the large scale phenomenology 
of turbulence may elucidate the underlying mechanism of turbulent motions at 
the small scales. Any useful theory of turbulence must incorporate some degree 
of anisotropy, and the purpose of our system is to provide a test bed. 
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APPENDIX A 

EVALUATING THE DIMENSIONLESS INTEGRAL A, 



The integral 

/oo 
(l-0|x|- (7+1) dx , (0< 7 <2), (A.l) 
-oo 

may be written as a sum of two integrals 

/0 ;>oo 
(1 -e")kr (7+1) dx+ / (1 -e")|x|- (7+1) dx. (A.2) 
-oo JO 

By making a substitution of variables y = —x in the first integral and employ- 
ing the trigonometric identity cos(x) = ^(e lx + e~ tx ), we can combine the two 
integrals to obtain 

/■oo 

A 1 = 2 (1 -cos(x))x~ (7+1) dx, (A.3) 
Jo 

Integrating the above integral by parts, the resulting boundary term can be 
shown to vanish at the boundary and we are left with 

2 f°° 

Ay = - x" 7 sm(x) dx . (A.4) 



7 Jo 

This integral can be evaluated by considering a complex integral of the form 

3 = j> 2~ 7 e iz dz , (A.5) 

and by applying Cauchy's theorem to the closed curve consisting of two line 
segments along positive Re z and Im z axes and two quadrants in the upper half 
plane; one very large with radius R and one very small with radius S, as shown 



in figure A.l The point z = has to be avoided; we do this by following the 
small circular quadrant of radius 5. The integral over the big circular quadrant 
C*2 tends to as R — > oo. The integral over the small circular quadrant near the 
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Re z 



Figure A.l: Contour for the integral A.5. The pole at zero has been avoided by 
following the small circular quadrant near the origin. The integrals over the two 
circular quadrants tend to zero as R — > oo and 8 — > 0. Only the integrals over 



the line segments contribute to the integral A.5 



origin C 4 can also be seen to approach as 5 — > 0. Only the integrals over the 
line segments remain and the complex integral reduces to 



a 



oo 



x~ 7 e ix dx - e^ 1 "^/ 2 / e~ y dy . 

i) Jo 



oo 



(A.6) 



Because we have avoided the pole at the origin, the integral S3 equals by 
Cauchy's theorem. After a substitution of T(l — 7) = f °° y~ 7 e~ y dy, the imagi- 
nary part of S3 yields 



J x 7 sin(x) dx = T(l — 7) cos^ 



7 IT 
2" 



(A.7) 



From Euler's reflection formula, T(l — 7) T("y) = n/ sin(7 7r), and the double- 
angle formula, sin(2 6) = 2 sin 9 cos 9, the above integral can be written in the 
form 

(A.8) 



x 7 sin(x) dx = -. , 

1 ' 2T( 7 ) sin( 7 7r/2) 
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Finally, with a substitution of the above in equation A.4 and an explicit use of 
the recurrence relation T(7 + 1) = 7 1X7), we find 



A, 



r( 7 + 1) sin(77r/2) 



(A.9) 



Gradshteyn and Ryzhik ( |2007 ) (3.761, equation 4) give 

[a>0; < |Re/i| < 1] . (A.10) 



x M 1 sinfax 



ax = sin — 

af* V 2 



The integral in equation A.4 is a special case of that found in Gradshteyn and 
Ryzhik| ( [2007| ) with a = 1 and y. = 1 - 7. 
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APPENDIX B 
STATISTICAL TOOLS 



In this thesis, unless otherwise stated, the following statistical estimators 
are used. For a sample (data set) with sample size iV and observed values 
xi, x 2 , ■ ■ ■ , xn, the sample mean, (x), is the sum of the observations divided by 
the sample size 

1 N 

i=i 

The measure of spread, or variability, in the sample is the sample standard de- 
viation 



<7~ = 



1 N 

— - (x)f 



N 

i=i 



1/2 



(B.2) 



The measure of spread, or variability, in the sample mean is the standard error 
(the central limit theorem) 

SE = 4^ • (B.3) 

y/N 



For a sample with N pairs of independent (xj) and dependent variables (j/j), 
if a straight line, jji = a + bx it is appropriate for the range of values studied, the 
linear regression coefficients, a and b, are those for which the sum of squared 
residuals, J^fiVi ~ Vif' ^ s minimized. Defining the following sums of squares 

N 



s xx = Y,( x i-( x )) 2 > ( B - 4 ) 

1=1 

N 

i=\ 
N 

s x V = ^2( x i - ( x ))(vi - (y)) , ( B - 6 ) 



1=1 
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the slope, b, and intercept, a, are given by 

, _ ^xy 



a=(y)-b (x) . 



(B.7) 
(B.8) 



Let \ji = a+bxi be the value predicted by the least-squares fit, then the difference 
between the original value and the predicted value is given by & = yi — y~i. The 
measure for the spread in £ s is the sum of squares of residuals 



1 N 



N-2 



i=l 



1/2 



or it can be rewritten in terms of the sums of squares 



q q c2 

°xx °yy ° X y 



1/2 



(N - 2) S xx 

The standard errors for the intercept, a, and the slope, b, are 



SE(a) = S£ 
SE{b) = 



1 (x) 2 
h 1 - L - 



1/2 



5*.,: 



(B.9) 



(B.10) 



(B.ll) 
(B.12) 



For a reference on linear regression, see e.g. Acton (1966) 
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APPENDIX C 

GENERATING EXPONENTIALLY CORRELATED COLORED NOISE 



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% This Matlab script calculates the exponentially % 
% correlated colored noise using the algorithm by % 
% Fox et al. (PRA, 38, 5938-5940, 1988) . % 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 



% Setting parameters for the algorithm 



a = 0; 
b = 0; 

lambda = 1 . / (0 . 1) ; 

D = 1 . / lambda ; 

N = 2. "11; 

f = 50; 

Fs = 3000; 

dt = l./Fs; 

t = dt . * [0 : 1 : (N-l) 



% Initializing random number 

% Initializing random number 

% 1/Correlation time (1/seconds) 

% Correlation time (seconds), D*lambda 

% Number of samples 

% Carrier frequency (Hz) 

% Sampling frequency 

% Sample time 

; % Time array 



exp (-lambda . *dt) ; 



% Here begins the algorithm 

e = zeros (N, 1); % Initializing noise array 

for i=l : (N-l) 

a = rand ( 1 ) ; 

b = rand ( 1 ) ; 

h = sqrt (-2 . *D . * lambda . * ( 1-E . " 2 ) . *log (a) ) .*cos(2.*pi.*b); 
e (i + 1) = e (i) . *E + h; 

end 
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28 % Modulate the amplitude of the carrier wave 

29 amp = e . *sin (2 . *pi . * f . *t ) ; 

30 

31 % Display the signal 

32 figure; 

33 plot (t, amp, 1 k- ' ) ; 

34 xlabel ( ' Time (s) ' ) ; 

35 ylabel ( 'Voltage (V) ' ) ; 

36 

37 % Take the Fourier transform of the signal 

38 Nfft = 2"nextpow2 (N) ; 

39 famp = f ft (amp, Nfft) /N; 

40 freq = (Fs . 12 ) . *linspace ( , 1 , Nf f t /2 + 1 ) ; 

41 y = abs (famp (1 :Nf ft/2 + 1) ). "2; 

42 

43 % Display the spectrum 

44 figure; 

45 loglog(freq, y, ' ko ' ) ; 

46 xlabel (' Frequency (Hz)'); 

47 ylabel (' [Amplitude r 2 (A.U.)'); 

48 

49 % Get rid of noise in the spectrum 

so windowsize = 30; 

51 yf = filter (ones ( 1 , windowsize) /windowsize, 1 , y) ; 

52 

53 % Display the filtered spectrum 

54 figure; 

55 loglog(freq, yf, ' ko ' ) ; 

56 xlabel (' Frequency (Hz)'); 

57 ylabel (' | Amplitude r 2 (A.U.)'); 
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APPENDIX D 

CALCULATING STATISTICS WITH INTER-ARRIVAL TIME WEIGHTING 



1 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

2 % This Matlab function calculates statistics 

3 % according to inter-arrival time weighting 

4 % 

5 % [out, sigma, num] = InterArrival (t, v) ; 

6 % 

7 % Each vector is weighted by the time past since previous 

8 % measurement . 

9 % Measurements made after 5 times of the mean arrival time 

10 % are rejected. 

n %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

12 

13 function [out, varargout] = InterArrival (t , v) 

14 

is % Check for error in input arguments 

16 error (nargchk (1, 2, nargin) ) ; 

17 if (nargin < 2) | | isempty (v) 

is error (' Input argument must contain two vectors'); 

19 end 

20 

21 if length (t) js length (v) 

22 error (' Vectors must have the same length'); 

23 end 

24 

25 % Specify the criterion for rejection 

26 limit = 5; 

27 
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28 % Calculate the inter-arrival time 

29 dt = dif f (t) ; 

30 v = v (2 : length (v) ) ; 

31 

32 % Reject measurements with dt > 5*mean(dt) 

33 accept = logical (dt < limit*mean (dt, 1) ) ; 

34 

35 dt = dt (accept); 

36 v = v (accept ) ; 

37 

38 % Calculate the weighted sum 

39 out = sum(v.*dt) / sum(dt); 

40 

« % Provide the standard deviation and the number of statistics 

42 if (nargout > 1) 

43 varargout(l) = {std(v)}; 

44 end 

45 if (nargout == 3) 

46 varargout(2) = { length (v) } ; 

47 end 
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APPENDIX E 
RESAMPLING LDV VELOCITY SIGNALS 



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% This Matlab function resamples LDV velocity signal 

% tin will be rounded up to the nearest multiples 

% of tsep. 

% If there are more than one velocity within a particular 

% tin, an average of the velocities is calculated. 

% syntax: 

% [vout, tout, nout] = LDVResample (vin, tin, tsep); 

% inputs: 

% vin velocity vector 

% tin time vector 

% tsep separation time 

% outputs: 

% vout resampled velocity vector 

% tout resampled time vector in multiples of tsep 

% nout number of elements within a bin 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

function [vout, tout, nout] = LDVResample (vin, tin, tsep) 

% Check for error in input arguments 
error (nargchk (1,3, nargin) ) ; 
if nargin < 2 | | isempty(tin) 

error ('Time array is missing!'); 

end 
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if nargin < 3 | | isempty (tsep) 

error (' Specif y separation time.'); 

end 



if ischar(vin) | | ischar(tin) | | ischar(tsep) 
error (' Input arguments must be numeric.') 

end 



% Initialize vectors 

vout = zeros ( 1+round (max (tin) /tsep) , 1); 
nout = vout; 



% Take the sum of velocities within a bin 
for i = 1 : length (tin) 

j = 1 + round (tin (i) /tsep) ; 

vout ( j ) = vout ( j ) + vin(i); 

nout(j) = nout(j) + 1; 

end 



% Calculate average 
nout = nout + (nout==0); 
vout = vout ./ nout; 

tout = tsep* [ : round (max (tin) /tsep) ]' ; 
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APPENDIX F 

CALCULATING VELOCITY AUTOCORRELATION 



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% This Matlab function calculates velocity correlation 
% for equally spaced data 

% syntax: 

% [fout, stat] = LDVCorr(vl, v2, tl, t2); 

% inputs: 

% vl, v2 — velocity vectors 

% tl, t2 — time vectors, equally spaced 

% outputs: 

% fout — velocity correlation function 

% stat — [vl_RMS, v2_RMS, vl_mean, v2_mean, 

% nfout, nvl_mean, nv2_mean] 

% RMS — RMS of coincident velocities 

% mean — mean of nonzero velocities. 

% n's — number of statistics 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

function [fout, varargout] = LDVCorr(vl, v2, tl, t2) 

% Check for error in input arguments 
error (nargchk (1,4, nargin) ) ; 
if nargin < 2 | | isempty(v2) 
error ('v2 is missing!'); 

end 

if nargin < 3 | | isempty(tl) 

error('tl for vl is missing!'); 
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end 

if nargin < 4 | | isempty(t2) 

error ('t2 for v2 is missing!'); 

end 

% Locate nonzero elements 

cutoff = min ([ length (t 1 ) length (t2 )]) ; 

nonzero = find (vl ( 1 : cutoff) . *v2 ( 1 : cutoff) ) ; 

% Calculate the correlation 

foo = (vl (nonzero) -mean (vl (find(vl) ) ) ) . * (v2 (nonzero) -mean (v2 ( f ind (v2 ) ) ) ) 

( std (vl (nonzero) -mean (vl (find(vl) ) ) ) *std(v2 (nonzero) -mean (v2 (find ( v2 ) 
fout = sum (foo) /length (nonzero) ; 

% Provide the RMS, the mean and the number of statistics 
if (nargout > 1) 

stat = zeros (1,7) ; 

RMS1 = std(vl (nonzero) -mean (vl (find(vl) ))) ; 

RMS 2 = std (v2 (nonzero) -mean (v2 (find (v2) ))) ; 

MEAN1 = mean (vl (find (vl) )) ; 

MEAN 2 = mean (v2 (find (v2) )) ; 

NFOUT = length (nonzero) ; 

NMEAN1 = length (find (vl) ) ; 

NMEAN2 = length (find (v2) ) ; 

stat ( : ) = [RMS1 RMS 2 MEAN1 MEAN2 NFOUT NMEAN1 NMEAN2 ] ; 
varargout(l) = {stat}; 

end 
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APPENDIX G 

CALCULATING LEAST-SQUARES FIT COEFFICIENTS 

For an introduction to the theory of linear regression, the reader could con- 
sult IDrapirandSmith] fl998l . 



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% This Matlab function calculates the least 

% square fit coefficients 

% Usage: 

% [B] = LinearMultiVarFit (Y, XI, X2, XN, flag) ; 

% Input : 

% Y - fit function 

% XI, X2, ... XN - variables (2D arrays) 

% flag - 'zero' for regression through the origin. 

% default is no crossing at origin. 

% Output : 

% B - an array of fit coefficients 
% 

% LinearMultiVarFit calculates the fit coefficients of 

% the multivariate function 

% Y = BO + B1*X1 + B2*X2 + . . . + BN*XN 

% with matrix left division in the least square sense. 

% If X's and Y are M by N matrices, B(:,j) is a 

% column vector [BO Bl B2 ... BN] of fit coefficients 

% for each set of Y(:,j), Xl(:,j), X2(:,j), ... XN ( : , j ) 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

function B = LinearMultiVarFit (varargin) 
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% Check for error in input argument 
if (nargin < 2) 

error (' Input argument must contain at least two arrays' 

end 

% Check for the number of input arrays 
if isstr (varargin { end} ) 

nin = nargin-1; 

zeroflag = varargin { end} ; 

else 

nin = nargin; 

zeroflag = 'boobs'; 

end 

% Form an ND array 
for i=l : nin 

X(:,:,i) = varargin{i}; 

end 

[M, N, P] = size (X) ; 

% Initialize array 
if (M==l) 

X = permute (X, [2 13]); 

[NX, NY, NZ] = size (X) ; 

B = zeros (NZ, 1) ; 

else 

if (N==l) 

[NX, NY, NZ] = size (X) ; 
B = zeros (NZ, 1) ; 

else 

[NX, NY, NZ] = size (X) ; 
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57 B = zeros (NZ, NY) ; 

58 end 

59 end 

60 

61 % Compute fit coefficients 

62 if ( strcmpi ( zerof lag, ' zero ') ) 



63 B ( 1 , : ) = ; 

64 Mat = zeros (NX, NZ-1) ; 

65 for i=l:NY 

66 Mat ( : , : ) = X ( : , i , 2 : end) ; 

67 out = Mat\X ( : , i, 1) ; 

68 B(2:end,i) = out; 

69 clear out 

70 end 

71 else 

72 Mat = zeros (NX, NZ) ; 

73 Mat(:,l) = ones (NX, 1); 

74 for i=l:NY 

75 Mat(:,2:end) = X ( : , i , 2 : end) ; 

76 out = Mat\X ( : , i, 1) ; 

77 B ( : , i ) = out ; 

78 clear out 

79 end 
so end 
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